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SUR DE NOUVELLES SOLUTIONS OSCILLATOIRES DES 
EQUATIONS DE LA GRAVITATION 


ALPHONSE MATTE 


1. Introduction. Nous savons quedansle champ gravitationnel de la relativité 
(i.e. défini par les équations R,, = 0) les ondes dites gravitationnelles se propa- 
gept avec la vitesse de la lumiére. Nous allons maintenant voir, dans le présent 
travail que si, dans le méme continuum, on impose aux potentiels g,, une oscil- 
lation de fréquence assez élevée, il en résulte une onde du type maxwellien. En 
d’autres termes les éléments de courbure R,,,, vérifient certaines équations, a 
savoir, R*,s,,. = 0, que nous allons déduire directement de R,, = 0 et qui se ré- 
duisent aux équations de Maxwell suivant une approximation qui est précisément 
basée sur la valeur élevée de la fréquence. II se trouve alors que |'onde est 
maxwellienne en premiére approximation, méme si la courbure du milieu de 
propagation prend une valeur finie, et pour des oscillations dont |’amplitude 
peut étre aussi petite qu’on veut. Nous soulignons ce point parce que dans le 
cas du résultat classique concernant les ondes gravitationnelles la partie princi- 
pale de la courbure est précisément due 4 l’onde et celle-ci se propage dans un 
milieu qui est presqu’euclidien. 

En vue d’obtenir le résultat annoncé il apparaitra essentiel de séparer l’espace 
et le temps de sorte que nous adopterons, pour |’élément d’arc, la forme classique: 


(1.1) ds* = c*dt* — do’, 
avec 
do* = a,x'dx" (r,s = 1, 2, 3) 


dans laquelle les fonctions a,, et c dépendent des quatre variables x’ et ¢ (¢ pour 
x®). Nous verrons alors comment il est possible de remplacer les équations 
R,, = 0 par un systéme invariant par rapport aux variables d’espace x’ et dans 
lequel le temps figure comme un paramétre n’ayant plus le caractére tensoriei. 
Ce systéme fournira ensuite de nouvelles équations dans lesquelles les éléments 
de courbure R.g,, seront groupés en deux tenseurs E,, et H,, tridimensionnels 
(r, s = 1, 2, 3) et de second ordre qui, en premiére approximation, et du point 
de vue de la structure des équations, joueront le méme réle que les vecteurs E 
et H dans les équations de Maxwell. 


2. Notations. Nous rencontrerons des équations invariantes dans l’espace- 
temps ainsi que d'autres invariantes seulement dans les sections spatiales. Pour 
faire la distinction nous conviendrous que les lettres grecques employées comme 
indices doivent prendre les quatre valeurs 0, 1, 2, 3, alors que les caractéres 
latins seront limités aux seules valeurs 1, 2, 3. 


Recu le 13 juin 1951; revisé le 12 mars 1952. 
l 
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Par ailleurs tout symbole tensoriel surmonté d’une barre comme par exemple: 
Rn Regur, Pus 


représentera un systéme 4-dimensionnel de l’espace-temps, alors que pour les 
systémes correspondant des sections spatiales la barre sera omise. De cette facon 
il sera possible de distinguer entre un systéme de l’espace 3-dimensionnel tel 
que R,,,, et un autre systéme 3-dimensionnel comme  - représentant une 
partie des composantes du tenseur Rag,,. 

Cependant nous emploierons deux lettres différentes pour représenter les 
tenseurs métriques: a savoir g,, pour l’espace temps et a,, pour les sections 
spatiales. I] est entendu qu’on peut toujours choisir les coordonnées pour que 
go; = 0 de sorte qu’on a |’élément d ‘arc (1.1) dans lequel: 


(2.1) Zoo = e’, Gry = — Brs- 


Alors le systéme a,, détermine la métrique de toute une famille de sections 
spatiales paramétrées suivant le temps ¢, et l’espace 3-dimensionnel fainsi 
déterminé sera souvent représenté par le symbole E(3). 

Un calcul facile donne: 


g=-— ca, a = ||a,,||, 
(2.1a) 2” a“ c, g’ on 0, 2" an a” 
Nous écrirons 
1 da,, 
(2.2) Ays= a ot 


Le systéme ainsi défini est un tenseur de E(3) puisque d’une part, les éléments 
dx‘/dx’* d’une transformation ne dépendant pas du temps et |’opérateur 0/dt 
d’autre part sont permutables. Cependant il est important de noter que la 
relation correspondante écrite au moyen des composantes contravariantes 
définirait le tenseur de signe contraire. Effectivement si on dérive les deux 
membres de l’identité 





a” ‘nn aut" i™ 
on obtient 
da™ mr ns OA mn $ da" r da™ 
FF ay Tom a Toe age 
d’ou 
ox da" = g™"q™ OAmn 
ot ot’ 
c’est-a-dire 
re _ 1da™ 
_ “= 2c Ot 


Ce tenseur est évidemment symétrique. 








arr = & 


——— > el 
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Dans goo = c’, la fonction c(x’,t) est un invariant de E(3) et son gradient 
peut se présenter dans nos équations invariants 3-dimensionnelles comme un 


vecteur. Nous adopterons la notation: 


a. aw 
ax" " at 


Pour ce qui est des symboles de Christoffel nous adopterons la notation: 


_ 1( afer ote: _ 2a) __» 
Pas = 5 (285 + ax ax’)? Pet = oT pr 
Alors avec l’aide de g,, = — a,, dans (2.1) nous avons immédiatement: 


= —T,,,. Par ailleurs 


et de la méme fagon il est aisé de vérifier les résultats qui suivent: 





1 0a,, 
Too = Clo, Ts = 2 a’ 
1 da,, 
Poo, =~ Ga Dos => > oH’ 
T',00 ns CCr, P vos = — Dye, 
o_ & o_ 1 
Poo _ Cc ’ Fes —_ c y 
(2.3) Too’ = a’*cl,, &r..’ = cA “a 
r 0 _ Cr r oo. nr ¢ 
r0 c ’ rs are: 


Dyes 


3. Les équations R,, = 0 sous forme invariante dans E(3). Considérons le 
systéme Ryrs, qui représente seulement une partie des composantes du tenseur 
R.s,»; nous allons voir que nous avons 1a un tenseur de E(3) aussi longtemps que 
rien n'est changé dans la séparation de l’espace et du temps. Cela veut dire 


que nous considérons des transformations de coordonnées du type 


Alors la relation 


5 . ax* ax” ax” ax’ 
urae = Heater art Oe” ag” ax 


donne immédiatement (puisqu’aucun terme significatif n’est obtenu lorsque les 
] ] 


indices muets prennent la valeur 0) 


RP... = R. ax” Ax" Ox" 
— come ox” ax” ax” ax'*’ 


ce qui prouve notre avancé. 
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Les autres composantes de A,,,, donnent encore deux de ces tenseurs 3- 
dimensionnels dont l’ordre est égal au nombre des indices qui n’ont pas été rem- 
placés par la valeur 0. La vérification du caractére tensoriel s’effectue comme 
précédemment et nous avons en premier lieu 


ax* ax’ ax” ax’ 


= _ Ox Ox Ox 
R 700; = Rasyr ax” ax” ax” ax” ’ 


ot les sommes donnent des termes non nuls seulement pour 8 = un = 0 de 
sorte que 
= = @x* dx" 
R’ 002 = Raoon ax” ax’ 
Enfin on montre, toujours de la méme facon que le systéme R,om, est un tenseur 
de troisiéme ordre de E(3). 

Notons qu’on peut toujours faire usage des propriétés de symétrie de Rags,» 
qui donneront par example 


Roms == Rormn = Ranro- 


En conséquence tous les éléments du tenseur Ryg,, se trouvent ainsi regroupés 
en trois tenseurs 3-dimensionnels de l’espace E(3). Le passage au composantes 
contravariantes sera illustré dans des calculs subséquents. 

Nous emploierons la méme technique pour séparer les dérivées Rag,»,¢ en un 
certain nombre de tenseurs 3-dimensionnels comme, par exemple, Ryrss,2. Ici le 
dernier indice est surmonté d’une barre pour spécifier qu’il s’agit de la dérivation 
covariante par rapport a l’espace-temps. Le lecteur n’aura aucune peine a 
s’assurer que notre raisonnement sur les transformations de coordonnées dans 
E(3) vaut encore pour établir le caractére tensoriel des cinq indices. Bien 
entendu ce nouveau tenseur est distinct de Ryrs»,, obtenu par la dérivation dans 
E(3) et nous verrons plus loin comment on peut développer le premier de fagon a 
tout ramener 4 la dérivation dans ce dernier espace. 

Occupons nous d’abord de la séparation des équations R,, = 0 tel qu’annoncé. 
Une partie du systéme s’écrit 
Ry = o” Rese = 0 
ou, avec l’aide de (2.1a), 

(3.1) 7" O"Bares + c*Rorso = 0. 


Pour le reste on a les équations 


> 
$ 
ll 
va, 
wD 
= 
& 
ll 
S 


qui donnent immédiatement 








ee 








c 
é 
I 
é 





rent w 
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(3.2) a**Ruone = 0, 
(3.3) a” Rete = 0. 


Ces équations représentent une forme invariante dans E(3) du systéme 
R,, = 0. On pourrait les écrire en fonction des tenseurs a,, et A,, mais cela ne 
sera pas nécessaire pour atteindre le but que nous poursuivons. 


4. Symétrie des nouveaux tenseurs. Nous utiliserons le tenseur 
wate, eres = ake res, 
ol ¢,,, = e”** prennent les valeurs suivantes: 
0 si deux quelconques des indices rst ont la méme valeur. 
1 si rst représente une permutation paire des nombres 1, 2, 3. 
—1 si rst représente une permutation impaire des mémes nombres. 

Nous tenons a souligner que nous avons 1a les deux types de coordonnées 
d’un méme tenseur puisque les deux systémes coincident en coordonnées norm- 
ales (a = 1). Il est donc entendu que !’on peut passer d’un A I'autre indifférem- 
ment, lorsque dans nos équations tensorielles certains indices sont élevés ou 
abaissés. 

Nous avons d’ailleurs: 

(4.1) Evert = 26," 
ou 5,* représente la matrice unitaire. La vérification est aisée et il en résulte que 
ce dernier systéme est un tenseur mixte. 

Un autre tenseur mixte qui nous sera trés utile s’introduit comme suit: 


(4.2) Sue” = Cum€ 
On trouve aisément pour la valeur des composantes: 
1 siusv, r=u, s=v, 
—-1 siuxv, r=%7, s=uUu, 
0 dans les autres cas. 
Posons alors: 


(4.3) 4E” = Rooms 'e™, 
(4.4) 2cH,’ = Rowse”. 


Puisque Reema = Rmnwe nous avons E*’ = E" et, bien entendu, cela entraine la 
symétrie des composantes covariantes: 


(4.5) Ey, = Es. 


Les équations (4.3) peuvent étre résolues pour les composantes Ryyema: multi- 
pliant les deux membres par €,ap€s2y et faisant usage de (4.2) on a 


AE canrecys = Russ”, 


et puisque les seules valeurs utiles pour les indices (uv) sont (ab) et (ba) avec les 
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valeurs correspondantes 6,°°=1 et 6°*= —1 onobtient 


AE” cay rrys = (Rosas ~  _ +. 


2 R, sadly ! 


Finalement on a 
(4.6) Rares = EF esp stsys: 
Par ailleurs (3.1) et (4.6) donnent successivement 


—25 urp uo mn 
c Rorso =a Rarse =aeE Eurm€son 


— ope £66 em 


= RS eeten = En tng "Qtr } 

-_ Ey"dn: = Eas: | 
= Ev. ~ Brelte « ' 
f 
Alors par contraction de la derniére égalité et avec l'aide de (3.3) on a ] 
¢*a"*Rorw = E,’ — 3E," = — 2E," = 0, | 

d’ot les résultats suivants: | 
d 

(4.7) E;" = 0, 
: P ( 
(4.8) Er, = C Rye, = 0" Rares ; 
On peut noter que Il’invariant EZ," s’annule par suite des équations de la 

relativité alors que le tenseur E,, est symétrique identiquement. Nous allons ( 


retrouver des résultats identiques dans le cas de H,, sauf que l’ordre est renversé. | 
De la relation (4.4) on tire successivement | ( 


WH" = aR ue”, 


( 

cH” ern = A Riousdnr” = 2a"*R oar, 
qui s’annule a cause de (3.2). Mais H’*e,,, peut s’annuler seulement si , 
(4.9) H” = H". | “ 
D’un autre coté . , d 
2cH,’ = Ryo” = — Rome’, d 

qui s’annule identiquement par suite de l’identité classique 

Resur + Rages + Rava, = 9. c 


On a donc 
(4.10) H,’ = 0. 








nl 








EQUATIONS DE LA GRAVITATION 7 


Enfin une autre relation qui sera utile par la suite, se déduit de (4.4) en multi- 

pliant les deux membres par €snn: 
2H," €smn = Rrowdmn” = 2R oma 

d’od 
(4.11) Rios = CH" énse- 

Les deux tenseurs que nous venons d’introduire sont a la base d’un nouveau 
formalisme qui va nous permettre de mettre en lumiére un aspect jusque 1a 
ignoré, de la structure des équations de la gravitation. A vrai dire nous nous 


intéresserons moins 4 ces équations qu’a certaines de leurs conséquences dif- 
férentielles avec lesquelles nous allons tout de suite faire connaissance. 


5. Un systéme d’équations fondamentales. Partons de I’identité de Bianchi, 
| + Betas + Beene = 0. 


Aprés contraction avec g“ les deux premiers termes s’annulent a cause de 
R,, = 0 et il reste 


(5.1) o Revers = 0. 


Occupons nous maintenant de séparer ces équations suivant la technique 
déja employée. Une partie du systéme peut s’écrire g“R,,,. = 0 ou 


(5.2) — a Rares + C Rovers = 0. 
Et pour la série au complet nous obtenons encore: 
(S.3) — a" Runso,5 + CRonso.5 = 0, 
(5.4) a” "Rows. = 0, 

(5.5) a” Rromo,5 = 0. 


Disons maintenant que c’est aprés avoir transcrit le systéme (5.2) a (5.5) en 
fonction des tenseurs E,, et H,, qu’on y retrouvera la structure des équations 
de Maxwell. Bien entendu nous avons ici affaire aux dérivées dans l’espace-temps 
et il nous faut d’abord développer les termes de facon 4 tout ramener aux 
dérivées dans E(3). Nous avons donc recours au développement classique de la 
dérivée covariante et nous avons par exemple 


) R ree «cp ep *«p ep 
Revesa _ oe < a Rarse = Pas Ruase = es Rauras — Pes | 
Comparons cette égalité avec 


a = ohhaca as a a Ter Retos ca Faso Mkiets 7 a 
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On en tire, avec l’aide de f,,' = T,,', 
(5.6) Reset ” | — + = a, coe Re _ 3, a 


Par ailleurs on a directement 
(5.7) Rosset = OR eon _ Poo’ Rares — To,"Resss — Po,"Rores hes Po."Rorse- 


Nous suspendons momentanément les calculs en vue d’introduire une méthode 
d’approximation qui va nous permettre de simplifier en négligeant un grand 
nombre de termes. Ces considérations feront l’objet de la section suivante. 


6. Ondes a haute fréquence. Nous considérerons d’abord le cas d’un onde 
dont la fréquence prend une valeur m trés élevée et c’est sur ce nombre m que 
notre méthode d’approximation sera basée. Comme la vitesse de propagation 
de nos ondes va se révéler égale a l’invariant c, nous pouvons d’avance profiter 
de ce fait et choisir l’unité de temps de telle sorte que les potentiels c* et a,, 
soient tous du méme ordre de grandeur. Nous adopterons donc le centimétre 
comme unité de longueur et nous choisirons comme unité de temps I’intervalle 
nécessaire 4 la lumiére pour parcourir cette distance, ce qui entraine évidemment 
que la valeur de la vitesse de la lumiére dans le vide est égale 4 un. Par ailleurs 
une fréquence de 10‘ signifie: 10‘ ondes au centimétre aussi bien qu’un nombre 
égal de périodes par unité de temps. Alors, en vue des considérations qui vont 
suivre il est important de noter que notre choix d’unités entraine le résultat 
suivant: non seulement les potentiels c? et a,, sont du méme ordre de grandeur 
mais encore leur taux de variation par rapport au temps (par suite de |’oscil- 
lation) est du méme ordre de grandeur que le taux de variation par rapport a 
la distance. En d’autres termes les dérivées du premier ordre par rapport au 
temps et celles par rapport aux variables d’espace prennent toutes des vaieurs 
qui sont du méme ordre de grandeur (sauf évidemment au voisinage des singu- 
larités; nous considérons le cas d’un domaine dans lequel la courbure reste finie). 
Pour fixer les idées remarquons qu’une fréquence m = 10‘ correspond a une 
onde dans |’infra-rouge, 4 peu prés a la limite du spectre visible. D’un autre 
coté, et pour éviter toute confusion avec ces données expérimentales concernant 
la lumiére, rappelons encore une fois que nous voulons simplement imposer aux 
potentiels g,, une oscillation de fréquence suffisament élevée et montrer qu’alors 
il en résulte une onde caractérisée en premiére approximation par des équations 
ayant méme structure que celles du systéme de Maxwell. 

Pour commencer, il sera commode de nous restreindre au cas ov I’intensité 
d’oscillation des potentiels est de l’ordre de n~' et nous prendrons n > 10*. 
Alors dans un domaine de l’espace-temps s’étendant 4 un nombre raisonnable 
d’ondes et de périodes les oscillations de nos potentiels peuvent étre représentées 
en premiére approximation par des fonctions sinusoidales. C’est-a-dire qu’il 
suffit de s’en rapporter 4 la relation 


(1 in 2 nz) = cos 24nx 
dx an” - ” 
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pour s’assurer que, si l’intensité d’oscillation des potentiels est du méme ordre 
de grandeur que n~', les dérivées du premier ordre de ces mémes potentiels 
prennent des valeurs du méme ordre que n°, c’est-a-dire qu’elles sont des grand- 
eurs “‘finies’’ du point de vue de notre échelle suivant les puissances de n. Par 
ailleurs les dérivées du second ordre sont des grandeurs du méme ordre que n et 
comme les dites dérivées se présentent linéairement dans les équations R,, = 0 
nous aurons pour une premiére approximation des termes du méme ordre que n. 

Cependant nous sommes plutét intéressés au systéme (5.2) a (5.5) dans lequel 
les potentiels apparaissent avec leurs dérivées jusqu’au troisiéme ordre, ces 
derniéres linéairement; ce qui veut dire que nous aurons des termes du méme 
ordre de grandeur que m?. Nous allons tout de suite nous assurer que ces termes 
sont les plus grands et que nous les avons tous en ne conservant que les dérivées 
du troisiéme ordre. Tout d’abord les coefficients a,, et c~? dans (5.2) sont au 
plus du méme ordre de grandeur que n° ou “‘finis.”” D’un autre coté un simple 
examen de (5.6), (5.7) et des calculs intermédiaires nous montre qu’il nous reste 
a décider au sujet de termes tels que 


a, ae To,"Roase- 


Or ces termes résultent du produit d’un facteur de grandeur “‘finie” comme 
Tnu* avec un facteur d’ordre un comme 


ef _ sas 


(6.1) R rao = ax" ax” 


+ Ses Bes fas Pee Bees 
Nous en concluons que les dérivées covariantes par rapport aux deux métriques, 
celle de l’espace-temps d’une part et celle de E(3) d’autre part ne différent que 
par des termes négligeables en premiére approximation puis-qu’ils sont d’ordre 
un au plus. De méme le second membre de (5.7) se réduit A son premier terme. 
Nous nous servirons de ces résultats pour transcrire le systéme (5.2) a (5.5) en 
fonction des tenseurs E,, et H,, et en cours de route, d'autres termes comme 
CoRorse Ct CoRowum Seront négligés toujours par suite des considérations précédentes. 
Nous emploierons le symbole — pour indiquer une égalité qui n’est vraie qu’en 
premiére approximation. 

Partant d’abord de l’équation (5.2) on passe aux dérivées dans E(3) et on 
or ensuite les deux membres par e”**; il en résulte (en tenant compte de 

alias Rorses 


Db so —2 mso 0 R, se 
a" Rurse.n é" +c * Fre ~ 0. 


Le premier terme peut étre tranformé avec l'aide de (4.6); on obtient (puisque 
€aur joue le r6le d’une constante par rapport a la dérivation covariante) 


x ne? = "E”. a 


25,"a"E*’ €aur 


ar at 
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Par ailleurs, utilisant (4.11) et négligeant des termes tels que coH," et 








n D€nse 
Br 
on a 
mse OR os m oH,” 0H,” 
2c5, at 2c a 
Alors notre équation peut s’écrire 
m _. 0H,” 
OE 5 tour + Co he ~0, 


ou encore 


m anr —l ir oH” 
Bent +6¢€ 6 at 


~ 0; 
et finalement par suite de relations telles que 


a” 0H,” we aH™ 
ot ot ’ 
on a 


Transcrivons encore une fois ce résultat en y ajoutant les équations (5.3) a 
(5.5) transformées de la méme facon, quoique beaucoup plus aisément avec 
l'aide de (4.8) et (4.11). On obtient pour le systéme au complet 


1 dH,’ 
c @é 
ues __ 1 OE," 
(6.2) {rus € can" 


a ~w 0, 
H”™ ,~ 0. 


Rest + ~~ 0, 





Il est pour ainsi dire superflu de souligner ici qu’on retrouve les opérateurs vec- 
toriels de la divergence et du rotationnel, dans les termes E"* , et E,y ,e“** de nos 
équations (6.2). En conséquence on reconnait aisément la structure des équa- 
tions de Maxwell étendue aux tenseurs E,, et H,,. 

Il importe maintenant de noter que dans les dérivées de second ordre E,, us 
les indices (uv) sont permutables en premiére approximation puisque, ainsi que 
nous allons le montrer 4 I’instant, la différence entre dérivées covariantes et 
dérivées ordinaires provient de termes qui tombent dans la partie négligeable. 
Effectivement nous avons vu que les éléments £,,,, sont du méme ordre de 
grandeur que 7? alors que les I,” sont finis comme n° de sorte que dans 


 - 
ox” 
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le premier terme au second membre prédomine avec une valeur de |'ordre de 
n®. On a donc 


OE’, OE’, 


E sue ™ ax" _ ax’ ax" ’ 





d’ou 
r —7 
E id E 2,ue- 


En conséquence on peut s’en rapporter aux calculs classiques pour résoudre 
le systéme (6.2) et caractériser les ondes au moyen de |’équation 


1 a°¢ 
(6.3) Ad — 2 oF ~ 0. 

Il est entendu ici que la fonction ¢ intervenant dans les coefficients de nos 
équations (6.2) et (6.3) est considérée comme constante en premiére approxi- 
mation et dans un voisinage qui peut contenir un certain nombre d’ondes (d’une 
part les potentiels oscillent avec des intensités qui sont de l’ordre de n= et 
d’autre part les variations de potentiel dues a la gravitation sont aussi trés 
petites sauf au voisinage des singularités). Cette restriction a un voisinage est 
évidemment rendue nécessaire par la courbure du milieu de propagation mais 
elle n’en premet pas moins de reconnaitre les caractéristiques essentielles des 
ondes électro-magnétiques dans ce nouvel aspect de la structure des équations 
de la gravitation. 


7. Extension des résultats. Jusqu’a présent nos résultats sont limités au cas 
ot les potentiels oscillent avec une amplitude de l’ordre de n~'. Pour une ampli- 
tude plus grande, par exemple du méme ordre que les g,,, notre procédé d’ap- 
proximation n'est plus valide car dans (6.1) les derniers termes sont du méme 
ordre de grandeur que les dérivées secondes des premiers termes. En réalité, il 
n’y a pas lieu de considérer ce cas extréme pour lequel les potentiels changeraient 
périodiquement de signe. 

Par contre nous allons montrer que l’onde reste maxwellienne pour des arupli- 
tudes aussi petites qu’on veut et indépendamment de la courbure du milieu de 
propagation. Considérons donc le cas od l’amplitude d’oscillation des potentiels 
est de l’ordre de n~’ avec r > 1 et posons gy» = Yu» + h,, de telle fagon que les 
h,, représentent en premiére approximation la partie oscillatoire de ces mémes 
potentiels. Effectivement on peut toujours s’arranger pour que les h,, soient du 
méme ordre de grandeur que l’amplitude d’oscillation, 4 savoir de l’ordre de 
n~', et que par ailleurs les y,, oscillent avec des amplitudes < n~"—"'. 

Notons alors les résultats suivants: 


7.1. Les dérivées d’ordre k des h,, oscillent avec des amplitudes de l’ordre de 
n~'+* et l’ordre de grandeur de ces quantités ne dépasse pas celui de leur ampli- 
tude d’oscillation. 


7.2. Les Yu» et leurs dérivées jusqu’au deuxiéme ordre sont des grandeurs 
“finies’”’ au plus. 
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Disons maintenant que sans nous inquiéter pour le moment de savoir 4 quoi 
cela conduit nous allons séparer les équations (5.1) de fagon a ne garder que les 
termes oscillant avec la plus grande intensité. Nous allons voir tout de suite que 
nous avons tous ces termes en ne gardant que ceux qui contiennent les dérivées 
du troisiéme ordre des h,, et qu’ainsi l’intensité d’oscillation sera de l’ordre de 
n~'+*. Effectivement un simple examen de (5.6), (5.7) et (6.1) nous montre que 
dans les termes ne contenant pas de dérivées du troisiéme ordre seulement un 
des facteurs peut étre une dérivée seconde des h,, et que tous les autres facteurs 
sont des grandeurs fimies au plus de sorte que les termes en question oscillent a 
ordre — r + 2 au plus. 

Par ailleurs on peut effectuer le développement de (5.1) en introduisant 
d’abord les E,,, et Hy,» comme dans la section précédente et, puisque les 
dérivées du troisiéme ordre interviennent linéairement, on peut dédoubler ces 
tenseurs d’une fagon qu’on peut représenter symboliquement par 


Eye = Em.o(y) + En. o(h) 


pour ne garder que la partie oscillant avec la plus grande intensité, soit a l’ordre 
—r+3. Les En (hk) et Hw,.(%) sont alors des tenseurs pour lesquels la mé- 
trique est déterminée en premiére approximation par les valeurs moyennes ‘up. 

Bien entendu les termes que nous voulons maintenant laisser de coté ne sont 
pas nécessairement plus petits que les autres (c’est plut6t le contraire pour 
r > 2), cependant si les équations (5.2) a (5.5) doivent étre vérifiées il faut que 
les termes de la partie oscillatoire se compensent entre eux en premiére ap- 
proximation et que leur somme (dans chaque équation) se réduise 4 une quantité 
n’oscillant plus qu’a l’ordre — r + 2. Enfin l’ordre de grandeur de ces derniéres 
quantités (ou sommes) ne dépassera pas — r + 2 si les h,, sont convenablement 
choisis. En effet, il n’y a lA qu’une question de partie constante en premiére 
approximation qui peut étre renvoyée dans les y,,. Bref on est de nouveau 
conduit aux équations (6.2) pour la partie oscillatoire et les ondes de petite 
amplitude possédent toujours les caractéristiques que comporte la structure 
maxwellienne. 

Il est encore possible d’étendre nos résultats au cas oi la fréquence n’est pas 
nécessairement trés élevée mais cette fois nous serons limités 4 un milieu de 
propagation presqu’euclidien. Posons donc 


c a 1 ao hoo, 867 — 6,’ _— Boss 


ot les h,, et hoo sont des quantités élémentaires du premier ordre. Nous avons 
donc recours ici 4 un procédé classique d’approximation qui a permis, en par- 
ticulier, de caractériser les ondes dites gravitationnelles au moyen de |’équation 
(6.3). Cependant I’introduction de nos tenseurs E,, et H,, va nous permettre 
d’apporter un élément nouveau 4 ce résultat. 

Dans une premiére approximation suivant cette méthode on ne conserve que 
les termes linéaires. Alors le lecteur n’aura aucune peine a s’assurer que les 
équations (5.1) se réduisent aux dérivées de troiséme ordre des h,, et qu’ainsi 
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on retrouve encore le systéme (6.2) avec la différence que les dérivées covariantes 
sont remplacées par les dérivées ordinaires. Nous en déduisons que dans le 
domaine des fréquences de la radio nos ondes sont encore maxwelliennes pourvu 
que la courbure du milieu de propagation soit assez petite. 

Ce n’est pas sans raison que nous faisons ici la distinction pour exclure de nos 
déductions les ondes gravitationnelles proprement dites. En effet dans ce cas la 
fréquence prend des valeurs extrémement petites (de l’ordre de 10-"* dans le cas 
du mouvement de la terre autour du soleil et avec l’unité de temps comme plus 
haut) ce qui introduit dans le probléme un aspect critique que nous ne voulons 
pas discuter ici. 


8. Conclusion. Ainsi donc la structure des équations de Maxwell se retrouve 
en premiére approximation et comme conséquence directe des équations R,, = 0 
pour caractériser dans presque tous les cas les ondes qui pourraient exister dans 
le continuum espace-temps de la relativité. En résumé si la fréquence prend des 
valeurs correspondant 4 |l’infra-rouge ou plus grandes nos ondes sont caractér- 
isées par les équations (6.2) et (6.3) indépendamment de la courbure du milieu; 
par contre si la fréquence est moins élevée nos équations restent valides seule- 
ment si la courbure est petite. Notons encore une fois qu’en derniére analyse ce 
sont les éléments de courbure R.g,, qui vérifient de telles équations et que cela 
est rendu possible simplement par la fagon, 4 tout le moins remarquable, suivant 
laquelle les dits éléments se laissent regrouper en deux tenseurs Z,, at H,, jouant 
le r6le des vecteurs E et H des équations de Maxwell. 

En conséquence Il’invariance du phénoméne ne fait pas de doute car les 
oscillations qu’on imposerait arbitrairement au réseau de coordonneés ne pour- 
raient @tre caractérisées par des équations invariantes. Nous avons donc bien 
une onde et nous faisons tout-de-suite le rapprochement avec les ondes électro- 
magnétiques. Or ces derniéres se laissent représenter au moyen de vecteurs, du 
moins dans les aspects que nous leur connaissons, alors que pour les nouvelles 
ondes les tenseurs E,, et H,, présentent tout de méme une structure plus com- 
plexe et laissent entrevoir de grandes possibilités. 

Pouvons nous dire maintenant que, dans le continuum R,, = 0, nous avons 
identifié les ondes électro-magnétiques sous une forme nouvelle dont nous ne 
connaissions que quelques aspects? I faut avouer que pour le moment tout essai 
d’interprétation dans ce sens pose des problémes dont la solution ne semble pas 
encore 4 notre portée. Tout d’abord les ondes que nous venons de caractériser 
existent-elles? Ou, plus précisément, existe-t-il des solutions de ce genre qui 
soient réguliéres a l’infini? Voila une question qui nous laisse perplexes devant 
la complexité des équations R,, = 0. Et pourtant il serait du plus haut intérét 
de savoir déterminer les possibilités de ces équations en ce qui concerne des 
singularités accompagnées d’un champ oscillatoire stationnaire ou irradiant de 
l’énergie, et qui entrainerait sans doute un phénoméne de quantification. 
Quoiqu’il en soit nous terminons ce travail en exposant deux autres résultats qui 
prolongent encore l’analogie avec les équations de |’électro-magnétisme. 
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Note A. Dans cette section nous nous occuperons de deux invariants de 
l’espace-temps et, bien entendu toujours dans le cas de go, = 0, nous trouverons 
moyen de les exprimer en fonction des tenseurs E,, et H,,. 

Tout d’abord on vérifie sans peine: 


Risk” -< Rank” ry ) ae 


| + l—_ + Rw ht + 4R,o,k” 
(A.1) 


Rrak™™ + 4B, 0k” + 4Rok””. 
D’un autre coté, avec l’aide de (2.1a) et (4.8) ona 
Re = gg oo” Regus = € aa Ryo = CE", 
d’od 
(A.2) Rok” = E,,E”. 
Puis étant donné que Rous est un tenseur dans E(3) la relation (4.3): 
gt yO ae’ 
entraine (compte tenu de (2.1a) comme précédemment 
4E,, = RR enastuees 
d’od 
(A.3) 16E,,E" = R™ Rereyinn bus” = 4Renuh™’. 
Rappelons maintenant (4.4): 
2cH,” = Rio’. 
Nous avons aussi comme pour les relations précédentes et tenant compte 
de (2.1a) 
cH", = — OR ™ exes. 
Enfin par contraction il vient 
40°H,"H', = — Rough bmn” 
ou 
(A.4) 4H,,H" = — 2Riwk™™. 
Alors les relations (A.2) 4 (A.4) permettent d’écrire (A.1) sous la forme 
(A.5) Reak” = 8(E,,E" — H,H"). 


Nous ferons de méme avec un autre invariant dont l’expression fait intervenir 
le tenseur universel suivant 


— 
€afys = 2 Copy, 


~.——_ 
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avec €asy3 = 0 ou + 1 d’aprés la régle: 


0 si deux quelconques des indices afyé ont la méme valeur 
1 si aSyé est une permutation paire nes nombres 0, 1, 2, 3. 

—1 si a@yé est une permutation impaire des mémes nombres. 
L’ensemble des composantes covariantes non nulles de ce tenseur peuvent 
aisément étre représentées au moyen d’un tenseur de E(3) comme il a été 
expliqué dans §3. Le nouveau tenseur sera de troisiéme ordre puisqu’un des 
quatre indices aSyé doit prendre la valeur zéro. On notera d’abord que les 
permutations Orst et rst ont la méme parité (c’est-a-dire le méme signe) par 
rapport aux nombres 0, 1, 2, 3, et 1, 2, 3 respectivement. On a d’ailleurs gt = 
icat avec i = (— 1)! de sorte que 


€ors: = ICE ps4. 
Alors on a successivement 


Haier pe 
= Mel Goes 


RB’ Raan” enree + 2’ Rao” eure 
= 2B" Ran” tusoe + 2R Ran” Coors 

+ 4 Rao” eusos + 4B Rao” coors 
= SiR e” come + 4h’ Rn” toms 
4ic(2R™ Ryo” €ue + BR Ran’ Ess): 


Les deux termes de la parenthése se retrouvent 4 un facteur prés par la con- 
traction de E,, avec H"* de deux facons différentes A partir de (4.3), (4.4) et 
(4.8). Les calculs sont analogues 4 ceux qui ont donné (A.2) a (A.4) et entrainent 
finalement 


(A.6) RY’ Ras” eure = 32iE,,H". 


Il faut remarquer que les relations (A.5) et (A.6) sont possibles en 
autant que la forme (1.1) de notre élément d’arc est conservée. Cependant les 
invariants au premier membre nous permettent de conclure que les fonctions 


E,,E” = } A la E,,H™ 


sont elles-mémes invariantes par rapport a la transformation de Lorentz. On 
retrouve ainsi une propriété que la forme relativiste des équations électro- 
magnétiques confére aux invariants 


E’—H’, E-H. 


Note B. Partant de (6.1) nous ferons r = 0 et nous séparerons les sommes 
relatives a l’indice a. 
Alors 


R,,, = 2Eoet _ Poe pp np 4 By,°Py," — Po,°Ts,”. 
Ox Ox 
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Puis avec l'aide de (2.3), 
Run. = O(cA’ +) _ 9(cA",) 





+ c(T yA ‘s — ty ‘a + = - af”. 


Ox ax" 
aA", 0A", nat nat 
= 4 ax _ A") + (Tis A> Tt A s)- 


Ajoutons a ce dernier résultat le terme identiquement nul 
c(— r,, A", + I,,'A",), 


et nous reconnaitrons immédiatement le développement de deux dérivées covari- 
antes de sorte que 


(B.1) Rosy = c(A"s.s — A" s.9)- 
Pour v = n on retrouve les trois équations Ry, = 0 sous la forme 
(B.2) A*s.s — A*en = 0. 
D’ailleurs en tenant compte de (2.1) l’équation (B.1) peut s’écrire 

Raose = C(Ans,e — Anse); 
d’oi, aprés multiplication des deux membres par ¢*” et substitution au moyen 
de (4.4) 

2cH,” = c(Aas.e — Anes)€” = 2Ans,e €’”. 

On a donc finalement 
(B.3) H,! = Au.” 
et on retrouve la forme de la relation classique H = rot A étendue au tenseur 
Ars. 
Université Laval 
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THE COORDINATE CONDITIONS AND THE EQUATIONS 
OF MOTION 


L. INFELD 


1. Introduction. The problem of the field equations and the equations of 
motion in general relativity theory is now sufficiently clarified. The equations 
of motion can be deduced from pure field equations by treating matter as 
singularities, [2; 3], or from field equations with the energy momentum tensor 
[4]. Recently two papers appeared in which the problem of the coordinate 
system was considered [5; 8]. The two papers are in general agreement as far 
as the role of the coordinate system is concerned. Yet there are some differences 
which require clarification. 

For this purpose it will be useful to analyse a very simple and well-known 
example that will emphasize some of the essential ideas by which field and 
motion are connected. The problem of motion is complicated in the case of two 
bodies but much simpler in the case of one body, when we may use the geodetic 
principle, which, by the way, can also be deduced from the field equations 
{1; 6]. Therefore, we shall begin by returning to this rather elementary problem 
and analyse some of its features, which appear in a changed form in the two-body 
problem. An interesting by-product of this work is a method of deducing 
perihelion motion, which seems to me simpler and more instructive than any 
others I know. 


2. The Schwarzschild solution and the change of the coordinate system. 


It is usual to write the quadratic form of the Schwarzschild solution in the 
following way: 


r 


—1 
(2.1) ds = (1 _ 2m) _ (1 _ =) dr® — r°(d0” + sin’é d¢’). 


What are the characteristic features of this coordinate system? For large r 
the only difference between the Galilean and the Schwarzschild field consists of 
additional terms that are essentially the Newtonian potentials. For r+ @ the 
field becomes Galilean. As a rule it is in such a coordinate system that we 
calculate the perihelian motion. We ask: Is it possible to find a coordinate 
system in which the equations of the path are the ordinary equations of a 
conic without any perihelion motion? This does not mean that the perihelion 
motion can be wiped out. But it does mean that we can introduce an artificial 
observer moving with the angular velocity of the perihelion motion. 

To show this, at least for the mean path, it is sufficient to introduce a change 
in one variable, that is ¢. We write: 


Received March 28, 1952. 
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(2.2) ¢=V¥+ wy). 


We shall assume not that w(yW) is small, but that its derivatives are small, that 
is we shall consistently omit all expressions of higher order than dw/dy. This is 
our only assumption. 

Thus from (2.2) it follows that: 


- 2 fe) 
de = av(1 +27 , 


Therefore, in our new coordinate system, the metrical form (2.1) takes the 

following form: 
-1 

(2.3) ds* = (1 - 2m a - (1 ~ a) dr*® — r°de* — r’sin’o av(1 az), 

The next step is to calculate the differential equations of the geodetic line. 
This calculation is absolutely straightforward and little difficulty is added by 
taking (2.3) instead of (2.1) We begin with the procedure described extensively 
in any book on general relativity theory and take @ = $2. We then obtain three 
equations. The first is simply (2.1) or, in our coordinate system, (2.3) for @ = 4a 
and d@ = 0. Then we have two integrals. One of them which represents the 
conservation of momenta is, because of (2.2), 


(2.4) 


The third equation, also integrated, is 


(2.5) . (1 - —) 


or precisely the same as in the original coordinate system represented by (2.1). 
Thus (2.3), (2.4), (2.5) represent our equations of motion. 

We can write (2.4): 
(2.4a) r 


therefore: 
dy _h dw dy 


° ow ae dy ds 


Thus, neglecting higher powers of dw/dy, we can write (2.4): 


dy _ af . $e) 
(2.6) ds r i dy/° 


Besides this, we have equations (2.5) and (2.3) for d@ = 0, 6 = 42, which we 
shall rewrite now: 
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oan (1-8) G-2)"- A) 038) = 


Here we replace: 


-i1 

a by (1 - an) ; according to (2.5), 

dr dr dy 

a ™ deas’ 

“ by ( - de), according to (2.6). 
We obtain 

i 2)( r doy | 2) —— 

(2.7) c (* oe oe ef) 8 eS 


We could have done this still more simply by taking the usual equation in 
the coordinate system (2.1), that is the equation: 


_ $)\% - Hf _2=). _ 2m 
© (4 Pee ess 


and introducing into this equation the transformation (2.2). Equation (2.7) 
follows immediately. 

Now we write (2.7) in the following form, introducing as usual u = r~: 

7) 2 (#)(2) 
2 2 2.2 2 2 3 

—h\—) — hu -1 n\ =) (= = 0. 
(2.8) C (ae u + 2mu + 2h dy) \av + 2mh*u 

Our question is whether by a proper choice of the function w we can rigorously 
solve the equation (2.8) by 
(2.9) “= =r enn. 


This is indeed possible and to show it we split the equation (2.8) into two parts 
which we shall satisfy separately: 


2 
(2.10) Cc” — (2) — h*u® — 1 + 2mu = 0, 
2 
(2.11) C-c'+ 2w'(#) (2) + 2mh'u® = 0. 


Here (2.10) is the Newtonian equation, in which, provisionally, the constant 
C’ may be regarded as unknown. Its solution therefore is (2.9), where e and p 
are functions of C’ and h. Now as to the equation (2.11), we introduce into it 
(2.9), regarding it as an equation for w and C’. We obtain from (2.11): 


(2.11a) (42) sinty + rae + 3e cos ¥ + 3e’cos*y + e’cos* y) 
(C* — C*)p* | 
* = 
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or, assuming sin? y = 0: 


dw , 3mcosy 3m _, mecos'y 
GI®) Gs* oe ant b * bane 


1 (m , 3m oS") - 
+ sag f+ es" 


p 
Thus we can solve this equation, by putting: 
2mh* . 6mh* 
2.12 C®8 =a C+ SS + — 
(2.12) pre oP 
and 


dw , 3mcosy _ 3m , me cos'y _ 

(2.13) dy + pe sin” y p + psiny 6. 

Equation (2.12) says that the ellipse is slightly changed from the one which 
corresponds to the constant C. Remembering that 


h? is of the order r4w? 


(r being the “distance”’ of the planet from the sun) and that 


p is of the order of r, 


we see that 
C? — C” is of the order mv?/r 


(v being the “‘velocity’’). Therefore (2.10) in which C? was replaced by C” 
gives the “relativistic” correction to the “energy.” 


Equation (2.13) is not sufficient to determine w in the neighbourhood of 2rn. 
But it is sufficient to determine the perihelion motion. We have: 


$= ¥+ o(y) = f(y). 
Let us write 
on =f(2an) = 2x0 + w, 
bata = f(2an + 2x) = 2x(m + 1) + waa, 
then the definition of the perihelion motion is: 
nti — On — 20 = Wart — Op. 
Now integrating (2.13) from 2rn to 2x(m + 1) we have 
Wn+1 — Wy = Omx/p, 


which formula expresses the famous perihelion motion. 

In order that the perihelion motion describe reality, it must be, in a certain 
sense at least, independent of the coordinate system. Indeed, the perihelion 
motion is independent of the coordinate system, as long as the coordinate 
system is Galilean at infinity. 

The description of the perihelion motion of, say, Mercury is based on the 
assumption that the observer is so far away and in such a coordinate system 








ch 


“2 
, 
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that his field is Galilean. Otherwise, nothing else is required from the coordinate 
system. The one represented by (2.1) would do just as well as any other as long 
as it is Galilean at infinity. 

Imagine an observer (in a Galilean local coordinate system) very far from the 
“sun,” pointing with one rigid rod towards the sun and with another towards 
Mercury, at a moment when the angle between the two rods is smallest. In 
principle, this can be done, since, in a Galilean coordinate system, all the 
concepts—trigid rods, angles, etc., are the familiar Euclidean concepts. Then the 
rigid rod pointing toward Mercury will describe an ellipse. But after one rota- 
tion, the perihelion position (that is, the minimum angle between the two rods) 
will be slightly shifted. After very many rotations, the original perihelion position 
from which the idealized observer started will be reached again. We can ask: 
After how many rotations of Mercury will the perihelion point complete one 
rotation? Such a formulation of the question has nothing to do with any special 
coordinate system but it does assume that the chosen coordinate system is 
Galilean at infinity, because only in such a coordinate system can we use freely 
the concepts of rigid body, angles, etc. Thus we can have either Newtonian 
motion and a non-Galilean field at infinity, or perihelion motion and a Galilean 
field at infinity. 


3. The two-body problem. The above discussion should help to clarify the 
problem of the coordinate system for a two-body problem. It was solved in a 
series of papers, not by assuming the principle of a geodetic line but by deducing 
the equations of motion from the field equations. One of the essential ideas used 
then is the “new approximation method.”’ By its use we find the field in the 
second and third approximations. Then from this knowledge of the field we 
deduce the equations of motion in the fourth approximation. These are the 
Newtonian equations of motion. The next step is to find the field in the fourth 
and fifth approximations. Finally, we deduce the equations of motion up to the 
sixth approximation. These equations give us the perihelion motion of the two 
(or many) body problem. It would be extremely difficult, technically, to go 
beyond the sixth approximation. But it is not necessary to do so. All contributions 
to the motion in higher than the sixth approximation can be wiped out by a 
proper choice of a coordinate system [5]. 

Thus having described the procedure in general terms we shall write down the 
field equations in the second and third approximations [3, p. 219]. These are: 


(3.1) Jou =9 
= Yom, es + Yoome = JO. Om m,s = 1, 2, 3, 
where 
Yur = hye — Anon he, zp, »,o,p = 0,1,2,3 


Zur = Nur + hy» 


_— Sus; jo = 1; tim = 0. 
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Here, and later, Greek indices take the values 0, 1, 2, 3 and Latin indices the 
values 1, 2, 3. One more remark: the stroke means differentiation: 


Ogu» 
Yur: = ro = Zur.s 


Yur = ey chery = ALur,0 (ax” = T). 


Thus the differentiation with respect to 7 raises the order, since \ is the para- 
meter of the development and its order is indicated by a number below the y’s. 

The essential point in the use of the “new approximation method” consisted 
in the fact that we began the development of ‘oo, Yom, aNd Ym_ with Yoo Yom: 
and Yon respectively : 


Yo = d*yo0 + r*yo0 +... 
(3.2) Yom = dy0m + rym +... 
Yun = X*ymn + d"Ymn E+ + 


If we look at (3.1) we see that yoo must be a harmonic function. Let us choose 
an harmonic function that represents two moving singularities. But then yon 
is not uniquely determined. Indeed let us replace 


Yom by Yom + 4o,m 
3 3 3 
ao being an arbitrary function. Then yom + @o,m 7s a solution. Thus we could ask 
3 3 


whether the arbitrary function ao,,, can have any influence upon the equations 
3 


of motion. It certainly does not have any influence on the equations of motion 
of the fourth order (which are the Newtonian equations of motion). This follows 
from a theorem proved before [3, p. 233]. But a straightforward, though trouble- 
some, calculation [9] shows that it has no influence, either, upon the equations 
of motion up to the sixth order. Thus it is not true [8] as Papapetrou claims 
that we could have changed the equations of motion by changing the coordinate 
condition involving Yoo and ‘Yom- Thus nothing can change the equations of 


motion up to the sixth order, as long as we stick to the following procedures: 


1. To introduce solutions of Laplace’s equations for Yoo that represent the 
Newtonian fields. 


2. To use the “new approximation method”, starting with yoo, Yom, Ymn 
respectively. ’ : ’ 

Once these procedures are assumed and no singularities later arbitrarily 
introduced, the equations of motion up to the sixth order are uniquely deter- 
mined, giving us, in consequence, the perihelion motion of a double star. 





ly 
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But we have seen in the first part, that we can wipe out the perihelion motion 
by a proper choice of the coordinate system. There seems to be a contradiction 
between our statement here and the result of the first part of this paper. The 
solution is simple enough: to introduce a coordinate system that would wipe out 
the perihelion motion we would have to violate the second of our assumptions 
just quoted. 

Although the corresponding calculation is elementary, we shall show it in 
detail in order to clarify the problem 

Let us rewrite (2.1) in the form: 


| e 
(3.3) ds* = (1 - 2m) a _ ( _ 2m) dp’ — p'(dé” + sin’ d¢’). 
We introduce the transformation 


(3.4) p=r+m=r+m, 


since m is of the second order. Again we rewrite (3.3) keeping only expressions 
up to the second order: 


(3.5) ds* = (1 — 2mr™)dt? — dr*(1 + 2mr™)((dx')*? + (dx.)" + (dx*)*) 


9 


(Here the connection between 7, 6, ¢@ and x', x?, x* is obvious.) 


Thus we have in (3.5) 


hoo = = mr: = = Sinn2mr* 

2 2 2 2 
Therefore: 

oo = thoo + $h,, = — 4mr™, 

2 2 

Yom = 0, 


Yun = hun — b5mnlss + b3nnhoo = 0. 
2 


Thus this solution satisfies our conditions, because ym, = 0 and because Yoo 
represents a Newtonian field belonging to a particle. 

To wipe out the perihelion motiog we would have to introduce a rotation in 
the “plane” of the motion with an angular velocity w chosen so as to wipe out 
such a motion. Thus we introduce the transformation: 


1 1 . 
x = £cosw— e’sin w 


x” = t'sin w + tcos w 


In this coordinate system we have the following quadratic form: 
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(3.6) ds* = [1 — 2mr* + a (()* + (€)*) Jae + 2ehadt de? — 2¢° ade ae’ 
2 
— dr*(1 + 2mr~")((de")* + (dt*)* + (ae’)’). 
2 
If @ = @ is of order one, we see that our quadratic form violates both the 
1 
conditions on which our approximation procedure is built. First, yoo does not 


vanish for r— ©, and does not represent a Newtonian field of a particle. 
Secondly, Yo, does not start with yo, but with yo,,. Thirdly, 
3 1 


Yan = 40°((E)” + EY") Sms 
that is, Y¥nn does not start with ymax. 
4 


Thus, returning to the two-body problem, we see that the conditions of our 
method are sufficiently stringent to insure uniquely the equations of motion up 
to the sixth order without any additional conditions concerning the coordinate 
system. 

Therefore the choice of a coordinate system has no influence upon the equa- 
tions of motion up to the sixth order. A straightforward calculation shows that 
this is true for Yom that the addition of an arbitrary function do,, does not 


change the surface integrals that determine the equations of motion. But as 
has been shown before [3, p. 13] the choice of the coordinate system in Yma, 
4 


Yoo, ‘Ymo, has no influence upon the equations of motion of the sixth order. 
4 5 


But it will have influence upon the equations of motion of the eighth order. 
Therefore it is possible to use such a coordinate system for Yom, Ymn» Y00. Yom: 
3 + 4 5 


so that the contributions to the surface integrals coming from A are zero. We 
8 


know [5], that such a coordinate system exists, though it would not be an 
easy task to find it explicitly. Similarly we can introduce such coordinate 


conditions in ‘Yma, Yoo, Ymo that the contributions to the surface integral A 
6 6 7 10 


will be zero. Thus it is possible to regard the equations of motion that we ob- 
tained before [2] as the exact equations of motion in a properly chosen coordinate 
system. We can get rid of all contributions to the surface integral beginning 
with A, but we can not get rid of A and A; these expressions, and also the 


equations of motion up to the sixth order, are determined not by the choice of 
our coordinate system, but essentially by the procedure concerning the use of 
our new approximation method. 
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MULTIPLE QUANTIZATION 
A. E. SCHEIDEGGER 


1. Introduction. The efforts of most theoretical physicists of this century 
have been directed towards that branch of the physical science which is 
commonly called “Quantum Theory.” Physically, Quantum Theory was 
postulated because of a vast amount of physical evidence which led to the 
postulates of states, observables, superposition, and commutation relations. 
From these four postulates, all quantum mechanics follows. 

It is well known that the differential equations which are obtained in Quantum 
Theory can be quantized for a second time. There is a strong physical motive 
for doing this, namely, the necessity for a field description, required both on 
general theoretical grounds and to explain, for example, pair production. 
Mathematically the “‘second”’ quantization is expressed in exactly the same 
way as the first one. In this fashion a certain mathematical scheme has been 
set up for “quantizing” a differential equation. This scheme consists of several 
stages whereby the original dynamical variables are finally replaced by operators. 
This scheme has been applied to many examples. Because of the fact that it is 
mathematically rather involved, this was not an easy task. 

Once the mathematical scheme of Quantum Theory has been obtained, it is 
of considerable interest to investigate the latter even without any further 
reference to its physical background. It is possible that the mathematical form 
suggests interesting generalizations which, in the end, might even have a physical 
application. This is what we propose to do here. 

A careful investigation of the mathematical formalism of Quantum Theory 
shows that it can be regarded as a process which can be applied to almost any 
differential equation. Thus one might say that quantization is a corrective 
mathematical process which is applicable to unsatisfactory physical differential 
equations. This is an entirely new aspect of quantum mechanics which is quite 
apart from any physical arguments for introducing it. It is suggested by the 
mathematical form of the equations and it remains to be seen whether it has any 
physical applications. 

If we accept the interpretation of quantization as a corrective mathematical 
process, the possibility of repeated or multiple quantization presents itself, as 
it seems to us, quite naturally. There are, of course, no direct physical reasons 
for such a repetition of the quantization process, except for the fact that it 
yields physical results when applied the first two times. From the mathematical 
standpoint, however, a further repetition seems a most natural generalization. 
In addition, it might be of interest to know, even for a physicist, what the 
quantization process does to his equations when repeatedly applied. Perhaps one 
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might speculate that each quantization makes the equations describing a physical 
process more and more accurate. 

Thus, we propose to give first a detailed account of the mathematical induce- 
ments for the notion of multiple quantization. The idea thus established will 
be discussed in the rest of the paper. 

The results of the present investigations are as follows. 

As a first example the Schroedinger-Gordon field is quantized for the third 
time. To achieve this, it is necessary to obtain a Schroedinger-representation 
for the Schroedinger-Gordon field, since the customary Heisenberg matrix- 
representation does not permit a further quantization. Then, the energy eigen- 
values of the third quantization of the Schroedinger-Gordon field are obtained; 
they are found to be 4 >> fw, Nz. 

In order to quantize an equation for an infinite number of times, one has to 
know what the quantization of a Schroedinger equation with an arbitrary 
Hamiltonian yields (with every equation after a few steps one seems to arrive 
at a Schroedinger equation). It is found that in this case the (m + 1)th quanti- 
zation describes a Bose ensemble of systems as described by the mth quantization. 

In the case where the quantization is performed using anticommutation, 
instead of commuation rules, one obtains a Fermi, instead of a Bose ensemble. 
In this instance it turns out that the repetition of the quantization procedure is, 
physically, simply an easy way of doing the quantum statistics of fields. 


2. Inducements for multiple quantization: notation. First, we should like 
to show how those ideas which we intend to discuss in this paper suggest them- 
selves quite naturally. In order to do so, we have to summarize the principles 
of quantization in their mathematical form, emphasizing especially their 
algebraic features. At the same time, this will serve to explain the notation of 
the present paper. The physical reasons which require the mathematical steps 
will not be discussed here, since we are only concerned with the mathematical 
formalism. 

In classical theory one starts with a series of real generalized co-ordinates 
dx (k = 1, . .. , f) which describe the state of the mechanical system concerned. 
The co-ordinates g, are functions of a parameter ¢ which, in applications, may 
be identified with the physical time. Dots denote derivatives with respect to 
this parameter. The equations of motion follow from a variational principle 


(2.1) 6 fz dt = 0, 
where L = L(q,, Gy, ¢) denotes the Lagrangian function of the system. The 
variational principle is equivalent to the Euler-Lagrangian differential equations 


d aL aL 
(2.2) or ~e 


Instead of using this representation one may pass to the canonical formalism 
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in the well-known manner. The algebraic nature of the latter is seen if one 


introduces Poisson brackets between arbitrary functions F, G, .. . of gq and 
Px = OL/0G,. Then, the equations of motion may be written as follows 
(2.3) F = [F, Hi), 


H denoting the Hamiltonian function of the system. 

The Poisson brackets have some remarkable algebraic properties (see, e.g., 
Dirac’s book [1, p. 85-86]), which permit one to reduce any Poisson brackets of 
arbitrary functions of p,, g, to elementary Poisson brackets between p, and q,. 
Thus, if one adds to the algebraic rules of Poisson brackets 


(2.4) [Pe qi) = — db; 


(6 = Kronecker symbol), then the problem of motion can be solved by algebra 
alone. 

This fact makes the quantization of a classical theory possible. The dynamical 
variables are re-interpreted as algebraic entities for which only laws of multiplica- 
tion need to be defined. If the meaning of ‘‘Poisson bracket” is newly fixed in 
terms of the commutator of the two entities concerned, then the motion of the 
quantum mechanical variables is entirely determined by the Poisson bracket 
rules and (2.3). 

One can obtain a representation of the quantum mechanical algebraic entities 
by identifying them with linear Hermitian operators in a Hilbert space. Actually, 
a proper Hilbert space is not quite sufficient because one has to admit vectors 
of infinite length. Nevertheless, we shall refer in this paper to such a general 
space simply by calling it “Hilbert” space. This draws one’s attention to the 
elements of such a generalized Hilbert space. Mathematically, they are a kind 
of vectors denoted by ¢, to each of which a dual exists (¢*) so that a scalar 
product (y*¢) may be formed. Physically, they are interpreted as representatives 
of the mechanical system such that 


(2.5) A = (¢*(Ag)) 


is the average value for the dynamical observable A in a statistical ensemble of 
measurements. 

The representation of quantum mechanics in this form corresponds to the 
Heisenberg-representation. By throwing the time-dependence from the opera- 
tors into the Hilbert vectors g, one can pass over to the Schroedinger representa- 
tion. Thus, one arrives at a differential equation for the time-dependence of the 
Hilbert vector ¢: 


(2.6) e=+Hy. 


hh 


In short, we note that starting from a “vector” 


(2.7) q = (qi. -- 4), 


ee 
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which satisfies a differential equation (2.2), one finally arrives through the 
process of quantization at a new differential equation of a similar type. The 
vector g originally characterizing the state of the system is replaced by the 
vector ¢ serving the same purpose. 

Thus, quantization is mathematically a correcting process for “classical” 
equations. But since the results are equations similar to those from which one 
had started, it must be asked why this correction should be applied only once 
to the classical theory. 


3. General theory. We note that in order to make a quantization procedure 
according to the canonical rules established in the last section, one needs a 
Hilbert vector g(k, ¢) which depends on a time-like parameter. The Hilbert 
vector must be subject to a differential equation of the type (2.2) regarding its 
time-dependence: 


(3.1) F(q(k t), G(R, t) G(R, t)) = 0, 


so that the canonical formalism can be set up. After this has been done, one can 


proceed to the Heisenberg equations of motion and hence to the Schroedinger 
equations of motion 


(3.2) * U(x) = Hv(x), 


which complete the circle from a Hilbert vector to a Hilbert vector. 

If (3.1) is of the type of a Schroedinger equation like (3.2), which will be 
usually the case if one wants to quantize for the second or third time, there are 
several difficulties to reckon with. First we note that x in (3.2) is a continuous 
variable with an infinite range, whereas k in g, assumes only a finite and discrete 
set of values. This difficulty, however, has been overcome successfully in the 
quantum theory of fields. 

Furthermore, the y’s in (3.2) are complex, and the differential equation (3.2) 
contains only first-order time-derivatives. The latter fact is insignificant, as 
(3.1) can be replaced by first-order equations; but the reality of the components 
of the y's is essential if one wants to maintain the analogy to the ordinary 
(first) quantization. We shall see, first, how this reality can be obtained by 
transforming (3.2), although there exists the possibility of abandoning the 
correspondence between y and a canonical variable, which is indeed sometimes 
done in the quantum theory of fields. 

In order to obtain real equations of the form of (3.1) from (3.2), one has to 
split ¥ into its real and imaginary parts 


(3.3) V=Wt we 


where both ¥ and yw are real. Then one has to try to set up equations for 
¥1 and y¥2 alone. Splitting (3.2) into its real and imaginary parts yields, if H is 
assumed to be real (which will be ordinarily the case), 
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a mo oh 
From the first of these equations one obtains 
(3.5) ¥2 = — AH, 
and from the second 

(3.6) vi = hAH-Y, 


which, inserted into either equation (3.4), yields 
(3.7) —-fh=-HHY, 


and an identical equation for ¥2. It maybe noted that the Schroedinger-Gordon 
equation is of this form (3.7). 

However, as it was said above, it is net necessary to transform (3.2) into real 
equations. If one keeps equation (3.2), one has to abandon the direct analogy 
between y and ordinary dynamical variables. Thus it may be possible to derive 
(3.2) from a Lagrangian principle wherein the Lagrangian function will depend 
on y and y* as well as on their derivatives. Then one can formally define the 
canonical conjugates * and x* and write down the Hamiltonian. Doing this it 
should be kept in mind that x and y are no longer canonically conjugate in the 
classical sense; it would be difficult to define Poisson brackets between func- 
tionals depending on them. One may, however, assume the same commutation 
relations as if x and y were conjugate. 

This treatment is very useful if one is only interested in the matrix formulation 
of the quantized theory. Then, x and y are simply represented by non-Hermitian 
matrices. This is sometimes so convenient that even in the case where one starts 
with equation (3.7), one transforms the latter into an equation where the 
variables are no longer real; this is well known from the transition to momentum 
space in the Schroedinger-Gordon case. Using complex variables, one is, however, 
generally at a loss if one wants to set up the Schroedinger formulation of the 
theory in order to proceed with further quantizations. In this case, one has to 
transform back to variables which are canonical. 

In the next section, we shall demonstrate these ideas in a model example. 


4. Example: Third quantization’ of the Schroedinger-Gordon equation. 
Assuming the Hamiltonian 


(4.1) H=cV Pp +m 


leads, after the first quantization, to the Schroedinger-Gordon equation (see 


(6]): 
1 a ; 
(4.2) la~§ t-s)eo0 


1The name “third quantization” has been used by Nambu [4] in quite a different context. 
Our paper has nothing to do with the work of Nambu. 


———$—S ~~ 
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which can be written in the following Hamiltonian form: 
(4.3) H =J$d*x; $ = 3{x* + clgrad y|? + cy}. 


If x, ¥ are taken as Hermitian operators, (4.3) represents the second quantiza- 
tion. A Fourier transform of the variables leads to 


(4.4) H = p> {Pa* Pe + on" g2* ge}. 


This equation is usually solved by the matrix method. If we want to carry on 
our quantization procedure for the third time, we must establish first a Schroe- 
dinger equation for the second quantization, instead of (4.4). The transition 
from a Heisenberg to a Schroedinger formulation is usually done in such a way 
that one of the canonical variables is taken as the argument of the new Schroe- 
dinger function and the conjugate one represented as a differential operator. 
Unfortunately, this procedure cannot be applied to (4.4) because p and g are 
not canonical since they are not real. They satisfy, however, the following 
relations (see [6]): 


(4.5) pi* = pari a = G-- 


The ordinary transition from y, x to g, p is thus not a canonical transformation. 
Therefore, one has to introduce for each “dimension” k a variable & which is 
canonical. If one finds a representation of p;, g, by differential operators acting 
upon the &’s such that the commutation relations and the conditions (4.5) are 
satisfied, then one has obtained the Schroedinger representation of the second 
quantization. It is easy to verify that this can be done by setting 





+ {( ) 
— = i(4) if tte setas 
$;: | ts] 0 \ 
(4.6b) De = 4 (hux)’ i \b -t.+ ab, + ae_S 


A straightforward check shows that the commutation relations for p, g are 
satisfied as well as (4.5). For the latter condition one has to remember that the 
star indicates the Hermitian adjoint of an operator and that the Hermitian 
adjoint of 0/0 is —0/dt. The procedure used here has been investigated by 
Infeld and Hull [2, §6] to which the reader is referred for the details. The total 
Hamiltonian of the field may be represented by 


(4.7) H = >> Mi, (é,? — 9°/0#,”). 


Since we can split up the eigenfunctions for every “dimension” k, it is possible 
to make the following ‘‘ansatz”’: 


(4.8) ¢ =[] exp (i Ey ) x (Ex). 
k 


Then, the Schroedinger equation 
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h ° 
(4.9) H¢= 7? 


yields the following condition for the ¢,'s: 

(4.10) Hic (& — 9°/O&:) ex = Ex or. 
As is well known, equation (4.10) is soluble only for 
(4.11) E, = }ho, (2N; + 1), 


which is the same result as obtained by the matrix method [2, p. 40]. A general 
state vector is obtained by a linear superposition of the particular ones in (4.8). 

The third quantization of the Schroedinger-Gordon equation is performed if 
we are able to quantize (4.9), or (4.10). If we keep the time dependence in the 
g's, (4.10) may be written 


(4.12) hia, (& — 3° /d&') oe = 


So. 
qm 


The task is to quantize (4.12). We note that this equation is of the form of a 
non-relativistic Schroedinger equation, if we set 
(4.13) m = fi/w,; U = 4&7 fay. 


The Hamiltonian formulation of the non-relativistic Schroedinger field has been 
set up long ago [5, p. 338]. Using those earlier results, we can write down the 
Hamiltonian which refers to the “dimension” k: 


(4.14) H, = jiu f {grad m, grad gy + &? mm Gx} dé. 


Thus, the total Hamiltonian is simply the sum of the expressions (4.14) taken 
over all the “dimensions” k: 


(4.15) H = hid we (grad Te grad Cr + hy Te ¢x) dé,. 


The eigenvalues of the total Hamiltonian are of particular interest. We note 
that the eigenvalues of each H;, have been calculated; they are 


(4.16) (H,)...= > N,ME™, 


where the N,’s range from 0 to ~, and E, is eigenvalue of the equation 


2 
(4.17) { . bw, ie? + thu, e) wa(é) = EY, (é). 
Thus we have 
(4.18) E, = th,(2n + 1). 


Therefore, we find the following eigenvalues of H,: 


(4.19) (Hy)... = rod) N, (2n + 1), 
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and of H: 
(4.20) (H)... = >> Hw, N,” (2n + 1). 


We may observe that the N,’s are arbitrary integer numbers ranging from 
0 to ~. But then also 
(4.21) N, = >> (2n + 1)N,” 
are arbitrary integer numbers ranging from 0 to ~, and we might as well 
represent the eigenvalues of H in the following from 


(4.22) (H)... =X Hie, Ne. 


These eigenvalues are very similar to those which were obtained after the 
second quantization (cf. (4.11)). Speaking in terms of particles we can say that 
our equations describe an ensemble of 9, particles of the energy }fiu,, k ranging 
over all the lattice points in momentum space. 

However, it is interesting to observe that the “‘particles’’ have only half the 
energy of those of the Schroedinger-Gordon equation. Furthermore, the infinite 
term in the energy has disappeared, since the lowest state is 2, = 0 for every k. 
Therefore, it seems better to interpret (4.22) physically in a different way so 
that the usual! physical interpretation of the second quantization can be retained. 
It is easily seen that the present scheme can be taken as describing a Bose 
ensemble of fields (instead of an ensemble of particles) each of which satisfies the 
Schroedinger-Gordon equation. In this instance the new quantization is equiva- 
lent to considering Bose statistics of the meson fields. 

It will be seen in the next section that quantizing a Schroedinger equation can 
always be interpreted as taking a Bose ensemble of systems as present before 
the new quantization, and doing statistics with it. Thus our last remark is in 
agreement with more general principles. 


5. Infinite quantization. After having seen that it should always be 
possible to perform a quantization of a first- or second-order differential equation, 
one may wonder how often this process can be repeated. 

We may note that after a very few steps, one always arrives at a Schroedinger 
type of equation; i.e., 

i 


(5.1) =; H ¢. 


The arguments of the ¢-function are the parameter ¢ and the canonical variables 
before the quantization. This is due to the fact that one has to set up the Hamil- 
tonian formalism for the equation from which one wants to start, and the 
Hamiltonian formalism immediately leads to (5.1). One may check that with the 
Schroedinger-Gordon case where we arrived at a Schroedinger equation after the 
third step, as was seen in the last section. 
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Therefore, if we are able to find a recursion formula which connects the 
eigenvalues of "H in (5.1) with the eigenvalues of the Hamiltonian "*'H of 
the next quantization, then we are able to proceed with the quantizations as 
often as we like. 

In order to obtain that (5.1) is the Euler differential equation of a Lagrangian 
principle, the Lagrangian density must be of the form 


(S.2) L = th o*> + Bd, ¢*) 
so that the canonically conjugate z of ¢ is 


oy . 
(5.3) «= —: = th ¢*. 

00 
It would be rather difficult to find all the terms in the Lagrangian (5.2) explicitly 
so that such a general equation as (5.1) would follow from it. Nevertheless, we 
can find the Hamiltonian **'H in matrix representation, in terms of the matrix 
elements of "H, provided that we assume the existence of such a Hamiltonian 
formalism and the corresponding quantization. 

For, if such a Hamiltonian formalism and the corresponding quantization 

exists, then the consistency of (5.1) with the quantum formalism requires 


“ - i n+l = in 
or 
(5.5) [¢, """H] = "H ¢. 


Expanding ¢ into the complete orthogonal set of eigenfunctions of "H yields 
(5.6) = Do an Un (x), 
(5.7) x = th @ = th>. an* Un* (x). 


The canonical commutation relations for x and ¢ are satisfied if 


II 


(5.8) [an, Om’ | _ len’, Om’*] = 0; [an, Om’*| = bam’ 
holds. Inserting the expansions for and ¢ into (5.5) yields 
(5.9) _ Om Um (x) "HH — we 9 Am Um(x) = > "Em Um(X)Om, 


where "E,, denotes the mth eigenvalue of "H. We multiply (5.9) by u,* and 
integrate over x; one then obtains 


(5.10) a," —*'H a, = "E; a, 
which is a matrix equation for **'H. The solution is obviously 


(5.11) “'H =), "E, a;*a;, 
P| 


because of the following relation which holds in virtue of (5.8): 





-_—~ ee 
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(5.12) [ay, aad & f| = ym "Ej(a a;* —— a;* a; ay) = "E, lax, a,* ja, = "EZ. a,. 
i 


A diagonal representation of "*'/7 is obtained if we choose 


(5.13) 0 0 0 
0 vi V1 0 1 
aq = 0 J/2 Py a,* = V2 > N, = a,*a, = 2 
i : 0 
which obviously conforms with (5.8). (The matrices (5.13) refer to the index k 


only, a, is the unit matrix with respect to all other indices.) Thus, the diagonal 
elements of **'H become 


(5.14) ate) = > N, °F; m= (No N, o° ae 


This means that the (n + 1)th quantization simply describes a statistical en- 
semble of systems as described by the nth quantization, wherein the counting prescrip- 
tions of Einstein and Bose hold. (Each combination of N,'s is enumerated just 
once.”) The lowest eigenvalue is always zero, no matter what the "£;,'s are; it 
never diverges. 

If we perform one further quantization, the eigenvalues will be 


(5.15) HE, = Nad, Ni*Ex =D ME 


where Jt, again may assume any integer values from 0 to ~. 

Therefore we see that the eigenvalues of the energy do not change any more 
after one has had a Schroedinger equation and quantized it once more. Thus, 
even if one repeats quantization an infinite number of times, the energy eigenval- 
ues remain final after a few steps. Therefore the eigenvalues of the Schroedinger- 
Gordon equation will always be as given in (4.22), notwithstanding how many 
quantizations are performed afterwards. 


6. The anticanonical formalism. It seems desirable to find out what 
becomes of the formalism of the last section in the case of anticanonical quantiza- 
tion. 

The anticanonical quantization is distinguished from the canonical procedure 
by the fact that for x, @ anticommutation (rather than commutation) relations 
are assumed. The equation of motion, however, holds in its canonical form. 

In order to find the anticanonical formulae that correspond to what we have 
done in the last section, one can keep practically ail the equations that were 
developed there, except that one has to find a new representation for the matrices 
a, and a,*. These representations have been proposed by Jordan and Wigner 
[3] in connection with the quantization of the Dirac electron field. The same 
matrices may be used in the general formalism as required in the present investi- 
gations. Thus we choose 


*It has been known before that a Schroedinger equation with an arbitrary potential leads to a 
Bose ensemble of the non-quantized system [5]; but the proof had never been given for an 
arbitrary Hamiltonian, as far as I know. 
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(6.1) a =a(°}), 


(6.2) a,* = f° °), 


The index k labels the state of the system as before. The factors ™ are equal 
to +1 or — 1 according to whether the number of occupied states with labels 
n < k is even or odd. 

Using this representation it is easy to see that the Hamiltonian is again 
given by 


(6.3) ott = » "E,; a;* a;. 
| 


This satisfies indeed the matrix equation (5.10). 

This means that the (m + 1)th anticanonical quantization of a Schroedinger 
equation describes a statistical ensemble of such systems as are described by 
the mth Schroedinger equation. But now, since we are using the anticanonical 
formalism, we see that this is a Fermi-Dirac ensemble. 

This result holds under the same premises and assumptions as before, that is, 
the consistency of the formalism has to be presumed in each case. 
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THE PROPAGATION OF A PLANE SHOCK INTO 
A QUIET ATMOSPHERE 


M. H. MARTIN 


Introduction. In this note the problem of determining the propagation of a 
plane shock moving through a polytropic gas into an undisturbed body of the 
gas, leaving a non-isentropic disturbance behind it, is reduced to the solution 
of the problem of Cauchy for a Monge-Ampére partial differential equation of 
special type. The application of this result to specific examples is being studied 
with a view to later publication. 


1. First principles. In one-dimensional unsteady flow of a gas the velocity 
u, density p, and pressure p must satisfy [3, p. 508] the underdetermined system 


(1) p(u,u + u;) + p, = 0, (pu), +p, = 9, 


of partial differential equations implied by the physical hypotheses of conserva- 
tion of momentum and of mass. By multiplying the second equation by u and 
adding the result to the first, we arrive at the equivalent system 


(p + pu), + (pu), = 0, (pu)e + pr = 0. 


If u, p, p are functions of x, t which satisfy this system, we can infer the existence 
of two functions £, y of x, t defined by 


d= = pudx — (p+ pu’) dt, dy = pdx — pudt. 


Along a trajectory x = x(t) of a gas particle in the (x,t)-plane we have dx/di = 
u and consequently y is constant along such a trajectory. The curves w(x, /) is 
constant in the (x, ¢)-plane are accordingly the trajectories of gas particles and 
the function y may therefore be termed the trajectory function. 

One readily sees that the above system can be given the form 


d= = — pdt + udy, dy = pdx — pudt, 
which, in turn can be replaced by 
(2) dt = udy + tdp, dy = pdx — pudt, 


provided we set — = E + pt. 

The form of this system suggests that ¥, p be taken for independent variables. 
This amounts to introducing the trajectories and isobars in the (x, #)-plane 
as a system of curvilinear coordinates in this plane. Obviously this cannot be 
done if the trajectories and isobars coincide and we accordingly exclude from 
consideration those flows in which each gas particle retains a constant pressure 
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during its motion.’ The introduction of ¥, p as independent variables into the 
system (2) yields the relations 


u= &, t = &, xy = & bop +1, Xp = by fpp, 


where tr = p-' = r(y, p), is the specific volume. When x is eliminated by 
partial differentiation from the last two equations above we find that ¢ = 
t(¥, p) is a solution of the Monge-Ampére partial differential equation.? 


(3) tw ton — ten = Tp. 


Here +r = r(¥, ~) is an arbitrary function whose appearance reflects the 
underdeterminedness of the system (1). Once r(y, ») has been specified and a 
solution § = £(W, ) of this equation has been found an unsteady flow is presented 
by 


(4) fies x(¥,p) _ fic bun + r)dy + by top dp}, t= £(¥,0), is fy (¥,p). 


For a given y the first two equations present the trajectory of a gas particle 
in the (x, #)-plane parametrically in terms of the pressure p; the remaining 
equation provides the velocity u at each point of the trajectory, and the density 
may be obtained from +r = r(y, p). 

The arbitrariness in the function 7(y, ») entering into (3) may be limited by 
further hypotheses of a physical nature. 

To begin with, let us suppose that the gas is polytropic with the caloric 
equation of state [1, p. 10] 


S—S.)/es ~ 
p on e lee . ¢ v 


which we write in the form 
(S—S.)/¢, —n 


rTr=e -p', n = 1/7. 
Secondly, from the hypothesis of conservation of energy it follows that the 
specific entropy S is constant along a trajectory [1, pp. 15-16], ie., S = S(y), 
the function S(y~) being termed the entropy distribution function. Thus the 
arbitrary function r(y, ») must take the special form 


(5) t= 3(¥)p", n= 1/y, 


where 4(y) is an arbitary function, determined by the choice of the entropy 
distribution function. 


2. Progressive condensation shock. For a condensation shock, carrying 
in back of it the values u, 7, p of the velocity, specific volume, and pressure, 
traveling into a quiet atmosphere in which these quantities have the constant 
values to, To, Po, the shock conditions take the form [3, pp. 512-513] 


‘Although I have not investigated the point, I suspect that flows with this property form a 
very restricted class. 

*This equation was obtained earlier by the author as a minor formal modification of a 
method developed for the treatment of steady plane flows. See [2, pp. 149-150]. 
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(i) «=u +V (b — po)(t — 7), 
” (i) #= w+ mV @- Pom) 
Gi) 2 - Y+1re— y— Ir 
bo (y+1)r — (y — 1)r0’ 


where denotes the velocity of propagation of the shock with respect to a fixed 
plane. 


We shall prove the following theorem for a polytropic gas. 


THEOREM. Once the entropy distribution function S(W) is specified, the deter- 
mination of a condensation shock moving into a quiet atmosphere (and the flow 
immediately in back of it) reduces to the solution of a problem of Cauchy for the 
Monge-Ampeére partial differential equation 

bw to — bvp + nd(y)p"” = 0. 

If we substitute r from (5) into (6 iii), we find 

(7) ute p"** — 8(4) > + pote p* — u'8(v)po = 0, uw = (y — 1)/(y + 1)- 


This equation defines p as a function of y and the curve p = p(y) in the 
(Wy, p)-plane is termed the carrier. 
From (4) and (6 i) the values of & are prescribed along the carrier by 


ty = uo tv {ro — 5(v)p "(W)} {>(Y) — pol. 


To obtain the values prescribed for &, along the carrier, we employ (6 ii). First 


we observe, from (4), that 
7 d 
£=&+ /te 


and consequently, from (6 ii), 








= =V {ro — 6(¥)P “(Y)}/{o(Y) — Pol, 


to determine £, along the carrier up to an arbitrarily additive constant. Of the 
three shock conditions (6), the last serves to determine the carrier in the (y¥, p)- 
plane, and the first two provide the required Cauchy data upon the carrier. 

If = &(y, p) is a solution of this problem of Cauchy, the shock curve in the 
(x, t)-plane and the flow behind it is obtained, at least locally, by mapping one 
side of the neighbourhood of the carrier in the (Wy, »)-plane on to the (x, ¢)-plane 
by the first two formulae in (4). 
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CORRESPONDING POLYHEDRA IN THE THREE SPACES 
OF CONSTANT CURVATURE 


ERNST ROESER 


1. Introduction. The five Platonic solids can be drawn in elliptic or hyper- 
bolic space just as well as in Euclidean space. Their numerical properties are, 
of course, the same in all three. So are the various angles subtended at the centre. 
But the face-angles and dihedral angles are greater in elliptic space, smaller in 
hyperbolic. It is a special feature of the non-Euclidean spaces that we cannot 
change the size of a solid without changing its shape. The edge-length 2/, circum- 
radius R, and in-radius r are conveniently expressed as functions of the mid- 
radius p: the distance from the centre to the mid-point of an edge. With each 
non-Euclidean polyhedron we shall associate a corresponding Euclidean poly- 
hedron of a definite size, whose mid-radius serves as a parameter for expressing 
all the metrical properties of the non-Euclidean polyhedron. When an elliptic 
polyhedron and a hyperbolic polyhedron have the same parameter, the edge- 
length of the former is the complement of the angle of parallelism corresponding 
to the edge-length of the latter. 

For these purposes it is immaterial whether we take the space of positive 
curvature to be elliptic (single-elliptic) or spherical (double-elliptic). To simplify 
the formulae we shall employ the natural unit of length. In hyperbolic space this 
means that a horocyclic arc of length 1 has the tangent at either end parallel 
to the diameter through the other end. 


2. Corresponding polygons. Let us begin by recapitulating the results of an 
earlier paper [5]. 

In Euclidean space, spherical triangles are usually drawn on a sphere of 
radius 1, because then the geodesic lengths of their sides are equal to the angles 
subtended at the centre of the sphere. In hyperbolic space the same result is 
achieved by using a sphere of radius 


arg sinh 1 = log (/2 + 1). 


On such a sphere, consider a circle of angular radius ¢, < $2 or straight radius 
op, SO that 


(1) sinh o, = sin o, < 1. 


Through this circle draw a plane and a horosphere (say, for definiteness, the 
horosphere curved in the same direction as the sphere). Draw p tangents to the 
circle in its plane, so as to form a regular p-gon, {p}, of in-radius ¢, and side 
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2l,, say. With the same points of contact, draw  great-circular tangents on the 
sphere, and p horocyclic tangents on the horosphere, so as to form a spherical 
{p} of in-radius ¢, and side 2/,, and a horospherical {p} of in-radius ¢, and side 





2 arg sinh 1 








2.8 
P 








— 
Fic. 1 Fic. 2 


2/,. The relation between each in-radius and side can be computed from a right- 
angled triangle (Fig. 1). Using the appropriate trigonometry (hyperbolic, 
spherical, or Euclidean), we find 


° ° T T 
(2) tanh J, = sinh o, tan > tan /, = sin o, tan —, l, = o, tan rs 


p 


Since ¢, is the semi-chord of a circular arc o, and of a horocyclic arc o, (Fig. 2), 
we have 


(3) sinh o, = sin o, = a, 
[7, p. 62, with S = 1], and therefore 
(4) tanh /, = tan /, = /,. 


This last result can also be seen directly, since /, and /, are a circular arc and 
its tangent in a diametral plane of the sphere, while /, and /, are a horocyclic 
arc and its tangent in a diametral plane of the horosphere. In particular, when 
the hyperbolic {p} is “asymptotic” (with its vertices at infinity), 


lL = @, l, = 1. 


By (1) and (2), tanh J, < tan 2/p; therefore a hyperbolic {p} has no corres- 
ponding spherical {p} if 


rT 
tanh /, > tan —. 
p 
(This failure cannot occur when p = 3 or 4, but may occur for a pentagon or 


higher polygon.) By (4), 
tan /, = tanh J, < 1; 
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therefore a spherical {p} has no corresponding hyperbolic {p} if 1, > 4, that 
is, if the {p} is a spherical triangle whose sides are obtuse. Thus the only 
regular polygon for which the correspondence holds both ways, for every possible 
length of side, is the square. 

The equation tanh /, = tan /, may be expressed geometrically [4, p. 3] by 
saying that 2/, is the complement of Lobatschefsky’s angle of parallelism for the 
distance 21,. 


3. Regular polyhedra. Analogous considerations in hyperbolic 4-space yield 
a definition for “‘corresponding” polyhedra. Instead of a common in-circle, these 
have a common mid-sphere (sphere touching all the edges) with its radius 
measured three different ways. 

Let {p, g}, where (p — 2) (g — 2) < 4, denote a regular polyhedron whose 
faces are p-gons, g at each vertex, in Euclidean or non-Euclidean space, so that 
the number of edges is [2, p. 13] 


E = 2pq/ [4 — (p — 2) q — 2)). 


Let ¢ and y denote the angles subtended at the centre by a half-edge and by 
the in-radius of a face. Then 


cos + 3 sin r/h sin x/h 

cos ¢ = oes 3/P mm ¢ 2 3", tan 2 . 

sin 2/q sin 4/q cos 4/p 

(5) cos «/q sin 2/h sin «/h 
on f a =. sin yy = ——_, tan ¥ = ——, 

sin x/p sin r/p cos m/qg 


where h = +/(4E + 1) — 1 [2, pp. 19, 21]. 
Let R denote the circum-radius, p the mid-radius, 7 the in-radius, ¢ the in- 
radius of a face, 2a the dihedral angle, and 26 the face-angle. These can be 


























p 
Fic. 3 


computed by solving the three right-angled triangles shown in Figure 3. The 


results are collected in Table I. 
The p and r of [4] have been interchanged for the sake of agreement with 
Sommerville [7, p. 123]. Note, however, that his 
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are our 


2i, 2¢, 26, 2a, p, q; C. 


For the extension to m dimensions see [8, p. 22]. 

In terms of Schlafli functions [1, pp. 15, 23, 28], the volume of the spherical 
\p, g} is ES(4e — x/p, x/q, $4 — a), while that of the hyperbolic {p, g} is 
EiS(}ax — 2/p, x/q, 4x — a). 

TABLE I 
THE GENERAL REGULAR POLYHEDRON 

















Elliptic or spherical space | Euclidean space Hyperbolic space 
tan / = sin ptan@ l= ptan@ tanh / = sinh p tan ¢ 
tan R = tan psec ¢ R = psec? tanh R = tanh psec @ 
tan r = tan pcos y r= pcosy tanh r = tanh pcos 
sin ¢ = sin psiny o = psiny sinh o = sinh psin ¥ 
cot a = cos p tan y a=tr-y cot a = cosh ptany 
sin @ = sec! cos x/p 6 = x — x/p sin @ = sech! cos x/p | 











4. A digression on rectangular polyhedra. It is interesting to see which 
non-Euclidean polyhedra {p, g} share with the Euclidean cube {4, 3} the 
property that the edges at each vertex form an orthogonal trihedron. This 
requires 

q = 3, a= 6 = dr. 


Since @ is greaier or less than $4 — x/p according as the space is elliptic or 
hyperbolic, we have p < 4 in the former case and p > 4 in the latter. Thus we 
find an elliptic rectangular {3, 3} with 

cos p = cot y = V}, 
a hyperbolic rectangular {5, 3} with 

cosh p = coty = 7, r= 4(/5 +1), 
and also, if we allow p to be infinite, a hyperbolic rectangular {6, 3} inscribed 
in a horosphere. In each case the edge 2/ is given by 

cos/ or cosh! = +/2 cos x/p. 


Table II lists the chief properties of these polyhedra, with the cube of mid- 
radius 1 for comparison. (Note that the mid-radius of the rectangular dode- 
cahedron is equal to its edge: p = 2/.) 

By repeatedly reflecting such a solid {p, 3} in its faces, we obtain a regular 
honeycomb { >, 3, 4}, filling the whole space. 

{3, 3, 4}, when drawn on a hypersphere in Euclidean 4-space, has the same 
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vertices as the regular cross polytope §,, analogous to the octahedron £; [2, 
p. 121]. There are sixteen cells {3, 3}. By identifying antipodes we derive a 
honeycomb of eight rectangular tetrahedra filling elliptic space [7, p. 124]. 

{4, 3, 4} is the ordinary space-filling of cubes. 

{5, 3, 4} was described by Schlegel and Sommerville [6, p. 444; 8, p. 17]. 

{6, 3, 4} was described by Coxeter and Whitrow [3, pp. 426, 427]. 


TABLE II 
THE RECTANGULAR POLYHEDRA 








ye: | 
Elliptic or spherical {3,3} lEuclidean {4,3}| Hyperbolic {5,3} |Hyperbolic {6,3} 





R=1}n,  sin’R sinh*R=3 | R= 








p = in, sin? p =} p= 1 sinh? p=rT p= 
r=jr, sin?’r =} | P=} sinh? r = }r | r= © 
l=, cos2/=0 | 21 = 2 cosh 2} = 7 | cosh 2/ = 2 





5. Corresponding polyhedra. The formulae in Table I apply to polyhedra 
{p, q} in the three spaces separately. In each case we can assign a value to p and 
deduce all the other properties. Corresponding polyhedra are given by identify- 
ing the Euclidean p with the elliptic sin p and the hyperbolic sinh p, say 

sin Ps = Pe = sinh Ph» 
whence 
tan /, = ], = tanh I, 21, = 4a — II(2),), 


and sin ¢, = o, = sinh a. 


The remaining relations are not quite so simple, but we can express all the 
properties in terms of p, and /, ( = p, tan @¢), as follows: 


Pe Pe 
tan R, = —-.; sec 4, R, = p-sec¢, tanh R, = —;—> sec 9; 
*— V(1 — pe’) . = " * V(1 + pe) ° 
a ee ” a Se Lae , 
tan?r, = Wi = 52) cosy, 7’ =p,.cosy, tanhr, = Va + p2) cos ¥; 


cota, = Vi-p. tany, a, =4nr-—y, cota, = V 1 + per tan y; 
° = -. te T T ° ce ) a T 
sin@, = V 1+/, cos—, 06 =%3tr-—-, sinh, = V 1 —l1, cos-. 
p p p 
It is interesting to observe that 

cot? a, + cot? a, = 2 tan? y = 2 cot? a, 

and 
- 2 . 9 27 - 2 

sin’ 6, + sin’ 6, = 2 cos rs = 2sin 8@,, 

whence 
cos 26, + cos 20, = 2 cos 28,. 
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6. Limitations of the correspondence. In elliptic or spherical space, the 
largest possible {, g} has 
R, =p, 77,2 br, l, = ¢, CG, @ ¥, a, = }r, 6, = «/q. 


Referring to a table of the five Platonic solids [2, p. 293], we see that ¢ > }4 
for the tetrahedron {3, 3}, @ = 49 for the octahedron {3, 4}, and @ < } for 
the remaining three. Since 


tan /, = tanh J, < 1, 


an elliptic {p, g} has no corresponding hyperbolic {p, g} if 1, > 41, that is, if 
the elliptic {p, g} is a tetrahedron of edge greater than 41. 
In hyperbolic space, the largest possible {p, g} has 


R, =], = ©, tanh p, = cos ¢, sinh p, = cot ¢, 6 = 0, 


T rT 
tanh r, = cos ¢cos ¥ = cot — cot —, 
p q 
° ° T 
sinh ¢, = cot ¢sin ¥ = cot ? 
Tr 
cota, = csc ¢tan yy = tan q’ 


[1, p. 28]. Since tanh J, = sinh p, tan ¢ = sin p, tan @ < tan ¢, a hyperbolic 
{p, g} has no corresponding elliptic {p, q} if 


tanh /, > tan ¢ 


(which implies ¢ < 3x). Thus, for the correspondence to hold, a hyperbolic 
hexahedron or dodecahedron or icosahedron must not be too large. 

The only regular polyhedron for which the correspondence holds both ways, 
for every possible size, is the octahedron. 
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A SIMPLE BOUNDING FORMULA FOR INTEGRALS 
J. L. SYNGE 


1. Introduction and summary. In this paper I establish the following 
bounding formula for the integral of a function of m variables: 


A 
(1.1) Jra-—— > m, Fi + Ei + Es, 


: 


the “error terms” £, and E, being bounded by 
(1.2) IE,| < F(D—T), |E. < RMT. 


Here F is the function, D the domain of integration in Euclidean n-space (and 
also the n-volume of that domain), and F an upper bound of |F| in D. D is 
divided up into congruent cells with m + 1 vertices (segments on a line, triangles 
in a plane, tetrahedra in space), the total m-volume occupied by these cells being 
T, so that D—T is the n-volume of that part of D not covered by the 
cells. A is the m-volume of a cell. The summation runs over all the vertices of the 
cells, F, being the value of F at the vertex i, and m, the number of cells meeting 
at this vertex. M is an upper bound to |d?F/ds*\, for differentiation in all directions 
and at all positions in D, ds being an element of length in the direction of differen- 
tiation. The quantity & is a parameter with the dimensions [length}’, depending 
only on the shape and size of the cell, and given by 


(1.3) k=3(r, —P/A), P= fr dD, 


where 7» is the radius of the circumscribed sphere of a cell and P the polar 
moment of inertia of a cell relative to its circumcentre. In particular 


( forn=1, k= poh? (h = length of segment-cell) ; 
form =2, k= ge(s2 + 522+ 537) (51, 52, 53 = sides of triangular 
cell); 





14) ) torn m3, k= 1G¢e — RB) 


(ro = radius of circumscribed sphere, 
R = distance between circumcentre and centroid of the cell). 


For n = 1 we would naturally divide D into equal segments, these segments 
being the cells. Then D = 7, E, = 0, and (1.1), (1.2) may be written 


b 
(1.5) f F dx = th(Fo + 2Fi +... + 2F,1+ F,) + E:, 
|E,| <h*? M(b— a); 
” Received July 16, 1952. 
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here h is the length of a segment and Fo, F;,.. . the values of F at the end points 
and points of division. This simple formula is very easy to prove, being in fact 
a by-product in the proof of the Euler-Maclaurin formula [3, p. 142 (6)]. Its 
significance here is as the prototype of the general multi-dimensional formula 
(1.1). 

If we like we can write (1.1) in a more general form applicable to the case 
where the cells are not all the same. Or we can improve the bounds by using 
different upper bounds M in different parts of D. But these considerations 
complicate the formula and we shall confine ourselves to establishing (1.1). 

Geometrical intuition is a great help and the plan we shall adopt is to prove 
(1.1) for » = 2 and then recognize that the same method is available for n = 3 
and indeed for all values of m. 


2. Bounds for a function and its integral when the function vanishes 
at the vertices of a triangle. Generalization to » dimensions. As a first 
step towards the proof of (1.1), we consider a triangle ABC and a function ¢, 
continuous and with continuous first and second derivatives. We impose on @ 
the condition 


(2.1) @=0 atA, Band C, 
and we write as bound for the second derivative 
(2.2) \d*o/ds*| < M, 





this inequality holding for differentiation in all directions and at all positions in 
the triangle. 


Let G be the circumcentre of ABC and ro the radius of its circumscribed 
circle. Consider the two functions 


(2.3) &, =} M(r? — 7), o, = 4 M(r’ — r,?), 


where r is distance from G. Note that in ABC, 9%, is positive and %, negative, 
both functions vanishing at A, B, and C. If we draw three-dimensional graphs 
of these two functions, they enclose a paraboloidal lens, with A, B, and C lying 
on its edge. We note that for all directions and positions 


(2.4) @4,/ds* = — M, d@’@,/ds* = M. 
The essential theorem we require is this: By virtue of (2.1) and (2.2), the three- 


dimensional graph of } within the triangle ABC cannot pass outside the paraboloidal 
lens formed from ®, and %.; equivalently, 


(2.5) %<¢< 9; 


in the triangle ABC. 
To prove this, consider first the side AB. On it 


(2.6) *, (d, — ¢) <0, &, — @ = OatA and B. 
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Hence , — ¢ is positive (or zero) on AB. Similarly ¢ — ®: is positive (or zero). 
Thus on AB, and similarly on BC and CA, the inequalities (2.5) hold. 

Now take any point E on AB and join it to C. On EC we have the inequality 
(2.6) and the boundary conditions 


(2.7) ®,—¢=0atC, ,—¢>OatE. 


Hence ©; — ¢ is positive (or zero) on CE. Similarly ¢ — %, is positive (or zero). 
Since lines such as CE traverse the whole triangle, the inequalities (2.5) are 
established. 


It follows then that, for integration over the triangle, 


| 
(2.8) | f.eap i< fia: dD = }M(r,°A — P), 


p= f{rap, 
A 


where dD is an element of area, A the area of ABC and P the polar moment of 
inertia with respect to the circumcentre. We may also write this as 


|¢ 
(2.9) | J,eap | <kMa, k= 4(r09 — P/A). 


The extension of the above reasoning to three dimensions is very easy. We 
now take a tetrahedron ABCD, with ¢ = 0 at the vertices and the bounding 
condition (2.2) as before. The functions #; and , are defined as in (2.3), G 
being now the centre of the circumscribed sphere and rp its radius. The in- 
equalities (2.5) are established in the same manner as above, and we get again 
the formula (2.9), wherein dD is now an element of volume and A the volume 
of the tetrahedron. 

We can carry the same type of reasoning on into space of higher dimensionality 
or back to m = 1. In fact, (2.9) holds for m = 1, 2, 3,.... 

The factor & in (2.9) agrees with (1.3), and the special values shown in (1.4) 
are easily obtained from the following equimomental properties [2, pp. 23, 27]: 

(i) A straight segment of mass m is equimomental to particles each of mass 
m/6 at its ends and a particle of mass 2m/3 at its centre. 

(ii) A triangle of mass m is equimomental to three particles each of mass m/3 
at the middle points of its sides. 

(iii) A tetrahedron of mass m is equimomental to four particles each of mass 
m/20 at its vertices and a particle of mass 4m/5 at its centroid. 


3. Proof of the bounding formulae (1.1), (1.2). Consider a domain D in the 
plane and a function F in D with the bounds 
(3.1) |F| < P, |\@F/ds*| < M. 


Divide up D into triangles covering an area 7, and so leaving an area D — T 
between the triangulation and the boundary of D. Then 
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(3.2) J rap = [ rav+&, 
D T 

where 

(3.3) |E,| < F (D —T). 


Now define a function F’ throughout T by the conditions that F’ = F at 
every vertex and F’ is a linear function of the coordinates in each triangle. 
Thus, in the triangle (x1, y:), (x2, ye), (xs, ys), F’ is given by the determinantal 
equation 


| FP .« y 1 
F, Xi V1 1 - 
(3.4) 2 aaa 4 ee 
| Fs X3 ¥3 1 


where F;, F2, F; are the values of F at the vertices. 
Write ¢ = F — F’. Then @ = 0 at the vertices, and 


(3.5) d*o/ds* = dF /ds*, 
since F’ is linear. Hence, by (3.1), 
(3.6) \d*o/ds*| < M. 


Application of (2.9) gives, for integration over a triangle, 
(3.7) J dD | < kMA, 
| A ' 


where 4 is the area of the triangle and & as in (2.9). 

So far we have not assumed the triangles to be congruent with one another. 
Let us now make this assumption, so that k and A are the same for all triangles. 
Then, for integration over the whole triangulation, we have from (3.7) 


(3.8) J eap < RMT. 
T 
Now 
(3.9) [ rap =f rap+ { sap, 
~J/T - T 


and so, by (3.2) and (3.8), 
(3.10) J rap = J F’ dD + E, + E2, 
D T 


where E; is bounded as in (3.3) and 
(3.11) \E.| < RMT. 


To complete the derivation of (1.1) from (3.10), we have to evaluate Sr F’ dD. 
To do this, we integrate (3.4) over the triangle, obtaining 
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| [rap ZA gd A 
PF, 1 y 1 


(3.12) Fr; a » 2 | = 0, 





F; X3 ¥3 a. 
where (Z, 9) is the centroid of the triangle. This gives at once 
(3.13) fr dD = 4A(F, + F. + F;). 


Summing for all the triangles, we get 


(3.14) f rap = 4A>> m, F,, 


where 7 runs over all the vertices and m, is the number of triangles meeting at 
the vertex 7. When we substitute in (3.10), we get 


(3.15) J rap = 440 m, F, + Ei + Es, 
D i 


which is (1.1) for m = 2. 

For n # 2, the initial factor must be changed from 1/3 to 1/(m + 1), because 
in using the analogue of (3.12) in m dimensions the centroid of the cell has 
coordinates Z, 7, ..., where 


(3.16) oe — ee +t... +d. 


But otherwise the reasoning is the same, and so we may regard the bounding 
formulae (1.1), (1.2) as established for all values of n. 


4. Two special triangulations. In dealing with an integral over a plane 
domain, two triangulations are particularly simple: (i) equilateral triangles, 
(ii) isosceles right-angled triangles. For them we have: 


(4.1) Egquilateral triangle of side 2a: A = 34a’, k = 4a’. 


(4.2) Isosceles right-angled triangle with hypotenuse 2a: A= a?, k = 4a’. 
If we use WN triangles (all of one type or the other), (1.1) gives 
(4.3) Equilateral triangles: 
f rap = 340°>° m, Fi + E: + E2, 
D i 


| E.| < $3'a‘ NM. 


(4.4) Isosceles right-angled triangles: 
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f rap = 40°) m, Fi) + EE. + E:, 
D i 
|E.|<4a'NM. 


At each vertex inside the triangulation, m, = 6 in (4.3) and m, = 4 or 8 in (4.4). 
On the boundary these values are reduced. If the domain D is such that it can 
be broken up into equal equilateral triangles or into equal isosceles right-angled 
triangles, then for such triangulations we have E, = 0 


5. The method of von Mises. Professor W. Prager has drawn my attention 
to a paper in which von Mises [1] develops a powerful method of approximating 
an integral over a plane domain by a weighted sum of the values of the integrand 
at selected stations. The error term depends on derivatives d*F/ds* and by 
taking sufficient stations we can make yu as large as we like. Thus the method of 
von Mises is more elastic than that of the present paper, for which only u = 2 
occurs. Moreover his method does not use triangulation. 

Let us compare the results given by the two methods in the simple case where 
the domain of integration is a triangle with sides s;, se, s3. Then (1.1), (1.2), 
(1.4) of the present paper give for the approximation and the error 


(5.1) J rap = 4D(Fi + Fr + Fi) + Es, 


|E:| <M D(s,* + s2’ + s;’). 


To apply the method of von Mises, let us take the stations at the vertices of 
the triangle. By formula (3) of his paper, weighting factors A;, A», A; are to be 
found to satisfy 

A, +A: +A; =D, 
(5.2) A, % + A2X2 + Az x3 = Di, 

A, w+ Ao 2 + As 3 = Dj, 


where (x, yi), ... are the vertices and (Z, 9) the centroid. The solution is 
(5.3) A, = A; = A; = 3D, 


and formulae (1), (II) of his paper give for the approximation and the error 
J rap =3D( + Fit F) +E, 


| Ey | < 3M [Po + $D ri? + r' + rs’)), 


where P, is the polar moment of inertia of the triangle with respect to the origin 
O and ri, f2, 73 the distances of the vertices from O. This bound for the error 
is minimized by choosing O at the centroid, and then we have 


(5.5) | Et | <M D(s;’ + s2” + s;”) + MP,, 


where P, is the polar moment of inertia with respect to the centroid. 


(5.4) 
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Comparison with (5.1) shows that the method of the present paper gives 
closer bounds than the method of von Mises. But of course it must be remembered 
that this particular example favours the method of the present paper, and that 
in applying the method of von Mises the stations need not be taken at the 
vertices, nor need there be only three of them. 
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CONSTANT HOLOMORPHIC CURVATURE 


N.S. HAWLEY 


We shall present in this paper a certain theorem concerning complex manifolds 
provided with an Hermitian metric satisfying the Kaehler restriction. The 
variables 2, Z2,..., 2, denote local complex coordinates in the manifold and 
2, Z2,..., %, their conjugates. The subscripts a, b, c,... run from 1 to m and 
by &@ we mean a + n (mod 2n), e.g., d = a. 

An Hermitian line element is given by 


(1) ds* = gasdzqd2, 


(summation over a, 6 from 1 to m); the coefficients gas shall satisfy the con- 
ditions 


(2) Zan = 2 = 9, Lab = Loa = GJan- 


In matric notation (2) has the following form: 


Zintiy +> + » Zi, 2n 
G= 


Zn.n+1s ++» 2n.2n ; 


If we let G’ be the transposed conjugate of G, then by (2), G’ = G and @, the 
matrix of the fundamental tensor, is as follows: 


OG 

6 - (9°), 
,_(0@ -(9°)-« 
e - (2°) - GO = ©. 


The line element (1) is said to satisfy the Kaehler condition if 


so we see that 


(3) Ofas _ Os 
. 02, O24 


The relations (2) and (3) thus imply 


Oger _ Ife» 
(4) 02. 3° 


We may form an affirie connection and a Riemann-Christoffel tensor from the 
fundamental tensor by the usual formulae. The only components different from 
zero of the Riemann-Christoffel tensor are 


Received August 22, 1951. 
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(S) Rated, Rivza, Ravea, Rasca- 


A two-dimensional element at a point is called an holomorphic section if it is 
tangent to a complex curve through that point. If the sectional curvature at a point 
is the same for all holomorphic sections at that point then the Riemann- 
Christoffel tensor can be written in the form 


(6) Rasea = — FD(ga5 Sea + Sat 0d), 
where b is the curvature on holomorphic sections; moreover b is a constant on 
the manifold if at each point the sectional curvature is the same for all holo- 
morphic sections at that point [1, p. 184]. For a more thorough-going exposition 
of this subject matter, one is referred to the papers of Bochner [1; 2]. 

Bochner has shown that relation (6) implies the existence of an analytic 
coordinate system in which the line element has the form 


2 idz,|" + bf p> lzel* Do \dz,|* — p> 2,d2,|°} 


(7) ; 
ds*s = 
(1+ 4b>0 |z,|*)° 





The line element (7) is, for b > 0, the Fubini-Study line element for complex 
projective space. If b = — 2 then (7) is the invarient line element of the unit 
cell |z;/? + . . . + |z,|2<1 under the group of all linear fractional transformations 
into itself. Let us call the complex projective space P* and the unit cell E* when 
they have these line elements. 

We can now state our theorem. 


THEOREM I. If a complex manifold S has an Hermitian metric satisfying the 
Kaehler restriction and if it has constant holomorphic curvature and if the space is 
complete in this metric then its universal covering space is analytically isometrically 
equivalent with E* forb = — 2 and with P* for b > 0. 


Let us assume from the start that S is simply connected. From Bochner's 
result we know that there is a local coordinate neighbourhood of each point of S 
such that the line element can be put into the form (7). Since S is complete there 
is a 6 > 0 such that each coordinate neighbourhood can be chosen with radius 
> 25. Let K(x, 5) denote the cell of radius 6 and centre x in S and let L(y, 8), 
denote the cell of radius 6 and centre y in E* or P* according to whether b = — 2 
or b > 0. Then K (xo, 5) can be mapped analytically isometrically onto L(yo, 4) 
where x» is any point of S and yo is any point of E*(or P* as the case may be). 
Let x; be a point of K(xo, 5) and let @ be the analytic isometry of K(x», 8) onto 
L(yo, 6) considered, and let y; = ¢(x,). We now wish to show that ¢ can be 
extended to an analytic isometry of K (xo, 5) UK(x:, 6) onto L(yo, 5)UL(j, 4). 

So far we have discussed only the mapping of a K(x, 5) onto an L(y, 4), but 
we could also discuss the analytic isometry of a K(x, 25) onto an L(y, 24) since 
one of these is contained in a suitable coordinate neighbourhood (in which the 
line element has the form (7)) about each point. We now observe that there 
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exists a mapping of K(x», 25) onto L(yo, 26) which agrees with @ on K(x», 8), 
for an isometry of K(xo, 4) is characterized by the directions at yo into which 
the directions at x» are sent. Since ¢ is analytic it sends holomorphic directions 
into holomorphic directions. Now the unimodular unitary group acts transitively 
on the set of all holomorphic directions and we have the unimodular unitary 
group acting on both K(xo, 26) and L(yo, 24) (leaving invariant the metrics in 
each). Thus an analytic isometry y of K (xo, 26) onto L(yo, 24) can be made into 
one which agrees with ¢ on K(xo, 5) by first performing ¥ and then performing 
an unimodular transformation of L(yo, 24) into itself. Therefore ¢ can be extend- 
ed from K(xo, 6) to K(x, 28) so in particular to K(x, 6) \U K(x, 4). Now let 
x's be any point of S then we can continue (analytically and isometrically) the 
mapping ¢ to x’». To do this choose a path C with end points x» and x’; take 
points x;,..., x, on C so that x; € K(xo, 4), x2 € K(x, 6),..., x0 € K (x,, 4), 
then by the method described above ¢ can be extended to each cell in turn and 
so to x’». In this way we extend ¢ to all of S. But ¢ is single valued since S is 
simply connected. Now S satisfies all the conditions of being a covering space 
over part of E* (or P*), however since E* (as well as P*) is also simply 
connected ¢ is actually an homeomorphism of S into E* (or P*). Also ¢ is an 
analytic isometry by the way it was defined, and it must be onto since S is 
complete; thus we have the theorem. 

Since we know that the complex projective space cannot cover, we have the 
following corollary: 


Coro.iary. If b > 0 then S must already be simply connected, hence S = P*. 


On the other hand, E* may well occur as the covering space of a non-simply 
connected S, even of a compact S. We now have 


TuHeoreM II. Jf S is a compact complex manifold with Hermite-Kaehler metric 
which has constant negative holomorphic curvature then the group of analytic 
homeomorphisms of S is finite. 


The proof of Theorem II follows quickly from Theorem I and [3, Theorem 
VII]. For S has Z* as universal covering space and E* is a bounded domain in 
E* hence a fortiori a Picard domain. (In Theorem I we assumed that b = — 2, 
but this was only so that Z* would be the unit cell; this is clearly inessential 
and 6 < 0 is all that matters.) 

We now state: if S satisfies the conditions of Theorem II it uniformizes an 
algebraic variety. For let § be a fundamental domain for S in E*, then § gene- 
rates a discontinuous group IT in E* such that E* (mod lr) = S. Given such a 
group, where S is compact modulo I’, Siegel has proved there exist m analytically 
independent automorphic functions in E* relative to [ [4, pp. 132—136]. To 
obtain an algebraic variety which S uniformizes we note that any m + 1 functions 
on E*, automorphic relative to I’, satisfy an algebraic relation (4, pp. 137 — 145]. 
Making this polynomial homogeneous we obtain a projective model of our 
algebraic variety. 
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Finally one notes that for S compact and b > 0 or b < 0 we have the existence 
of m analytically independent functions meromorphic on S. But for b = 0 we 
have no information as to the existence of functions. Indeed in this case S 
is a multi-torus which may have any number from 0, 1,..., analytically 
independent non-constant meromorphic functions. This phenomenon arises 
because of the need for period relations in the theory of Abelian functions; 
which in turn arises from the fact that in some multi-tori a local sub-variety is 
not part of a sub-variety in the large, not because the local sub-variety cannot 
be extended (it can be continued indefinitely), but because it winds infinitely 
often around the multi-torus and never meets itself (much as the familiar 
everywhere-dense integral curve winds on the two-dimensional torus). In fact 
a multi-torus (in m complex variables, m > 1) may be constructed which has no 
proper sub-varieties other than points. (For an example see [4, pp. 104—106].) 

In concluding I wish to thank Professor S. Bochner for suggesting Theorem I. 
It is also explicitly stated in [2, p. 21]. 
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HARMONIC p-TENSORS ON NORMAL HYPERBOLIC 
RIEMANNIAN SPACES 


G. F. D. DUFF 


Introduction. The subject of this paper is the study of boundary value 
theorems for harmonic p-tensors on a Riemannian space with an indefinite 
metric of the normal hyperbolic signature. The p-tensors or p-vectors ¢;,...«, 
are alternating covariant tensors of rank p, which are closely related to differential 
forms ¢ of degree p,  < m, on an m-dimensional manifold. 

On a Riemannian space with positive definite metric, the harmonic p-tensors 
have been studied by Hodge, de Rham, Kodaira, and others, and a theory 
which generalizes the classical potential theory, and possesses in addition 
certain new features, has been developed. 

The present paper continues a program of extending the boundary value 
theorems of the theory of partial differential equations to p-tensors. In the 
positive definite (elliptic) case, boundary value theorems have been given [3] 
for the Beltrami-Laplace equation 


Ag =0 


and for the harmonic field equations 
d¢=0, 6¢= 0. 


These equations are now investigated under an indefinite metric. 

The principal result is the solution of the Cauchy problem for the Beltrami- 
Laplace equation. After a preliminary section, the auxiliary conditions for the 
Cauchy problem are formulated, and uniqueness under these conditions is 
established. To construct the solution, I use the method of Riesz potentials. 
For the sake of brevity, reference to Riesz’ paper [9] has been made wherever 
possible. The solution of the Beltrami-Laplace equation thus obtained is then 
applied to the construction of solutions of the harmonic field equations. In a 
concluding section, some special cases, in particular the electromagnetic field 
equations, are examined, and the main theorems are compared with the results 
which hold in the elliptic case. 


1. Preliminaries. Let M be an m-dimensional orientable Riemannian space, 
with metric 


ds* = g,, dx‘ dx? 
having the Lorentz signature. That is, on being transformed to a sum of squares, 
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ds* contains one positive and m-1 negative terms. For convenience we shall refer 
to displacements as timelike if ds? > 0 and spacelike if ds? < 0. We shall also 
distinguish the timelike coordinate, say x”, by the letter ¢ in certain cases. 
This metric is also known as a normal hyperbolic metric (Hadamard), the 
corresponding Laplace equation being of the normal hyperbolic type. During 
most of this paper we assume that the metric is analytic, but this restriction 
can be removed. 

Throughout we shall use such terms as continuous, compact, convergent, on 
the understanding that they refer to a suitable positive definite distance, such 
as the Euclidean distance. 

We introduce on M the skew symmetric covariant tensors ¢,,... 4,; correspond- 
ing to these are differential forms of degree p: 

(11) @6=@= DD onde" A... A dx” = b4,...4,) de" A... A de”. 
s<...<t, 

The bracket enclosing the group of indices shall mean that the indices inside it 

are arranged in strictly increasing order. The differentials dx‘ are multiplied 

together by the exterior multiplication indicated by the sign A ; hence these 

differentials anticommute. We now define the differential operator d; 


(1.2) do = (ddxi,...1,)) dx" A... A dx” 
Thus d¢ is a differential form of degree p + 1. Let 


| 
| Bias see Bi. 4s | 
(1.3) D ticscted Pincede 


then 


a | Bins os Livis! F 


Pr... ae 

is just the Kronecker symbol often denoted by 
3, Jaseedp 
ta-eet ° 


Also let 
(1.4) Cir te...te = Ty..01 pee 


a) 14 
I 12...", 12..."| «+ 





We denote covariant derivatives of a p-tensor 
$1,...4, by D, Pi,... tps 
and set D‘ = g“D,. Thus 


e O¢d:.,... . 
(1.5) D194... i, = =, 2 ; fad nn or 


In the notation of [8a], 


(d$);,...44. = Pee "DO. 4,.. je) 
(1.6) 
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The operator d.d is identically zero. If dé = 0, ¢ is said to be closed; if @¢ = dx, 
¢ is said to be derived, and is therefore closed. The converse statement that if 
¢ is closed then ¢ is derived, is true “locally,” but not in the large. 

The dual *¢ of a p-form ¢ is an (m — p)-form whose components are defined 
by the formula 


(ty... ty) 


(1.7) (#D) go... jmme = Cbs. 0- ty) Saeedm—e 
If p is the degree of ¢, then 

eng = (— 1)""9. 
Hence, save for a sign, the * operation is its own inverse. It will be seen that 


¢@ and *@ are “perpendicular” forms at each point. 
The operation of dual derivation is given by 5, where 


(86) s....te—1 = (— 1) (nde) 4... 


a (SaoeeSed ry8 
m <= T'66,... toms "Do 4,...3 


- J») 


(1.8) 


It follows from the preceding remarks that 6.6 = 0. If 6@ = 0, @ is said to be 
coclosed, and if ¢ = dx, ¢ is coderived. A p-tensor ¢ which is both closed and 
coclosed will be known as a harmonic field [6]. 

The Beltrami-Laplace (B.L.) operator A for p-tensors is given by 


(1.9) —A=di + éd, 


and in expanded form may be written 


p 
(1.10) (Ad).,...4, = D'D ii... - > r, —— oR’ 5.5105, a or 


where R",» is the Riemannian curvature tensor and summation over the indices 
i, j, h, k and (j:...j,) is understood. If ¢ is a scalar, then A@ reduces to the 
usual Laplacian. A p-tensor ¢ whch satisfies the B.L. equation 


(1.11) Ag =0 
is said to be a harmonic form. 

From (10) it is clear that the B.L. equation is of the normal hyperbolic type 
under our Lorentzian metric. In fact, the equation (1.11) stands for a system of 


(™) equations, one for each component of ¢. Note that each of these equations 
has the same principal part, namely, 


tj 8° 4%... 4 
eax! ax? * 
Actually all terms containing first or second derivatives are of the same form in 
each component equation. Furthermore the form of the second order (principal) 
terms is independent of ~, so the theory of the characteristics will be carried 
over unchanged from the scalar case. The characteristic surfaces of each com- 
ponent equation being the same, we may speak of a characteristic surface of 
(11). These surfaces are given by 











where 
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(1.12) ae 0. 

Let O be an arbitrary point of M which we select as origin of coordinates. 
The geodesics of zero length through O form a conoid [5; 8] which is a charac- 
teristic surface (1.12). 

Let S be a sufficiently differentiable spacelike “‘initial’’ surface. For conven- 
ience we assume that S is compact. We construct a system of geodesic normal 
coordinates in which S is the hyperplane x" = ¢ = 0. Since the normals to S 
are timelike, we can write 


(1.13) ds* = dt” — gas dx* dx’, 


where the Greek indices range from 1 to m — 1. In this system [10], gin = 5,‘ 
D,, = D". 

Associated with S we define a region R as follows. R shall consist of those 
points P on the positive side (¢ > 0) of S, such that the retrograde characteristic 
cones (null cones) Cp with vertex P, together with that part of S intercepted by 
Cp, bound a simply-connected region of M. Denoting this region by Ds’, 
following Riesz, we see that if P € R, every point of Ds” also belongs to R. 
We may refer to R = Rg as the region of exclusive dependence upon S, in the 
sense that, as will be shown, solutions of the B.L. equation with data assigned on 
S are determined throughout R. 

Thus defined, Rs may contain non-bounding p-cycles. Any p-cycle of R is 
however homologous to a p-cycle of S, since a continuous cylindrical (p + 1)- 
dimensional surface can be constructed on the cycle, lying in R and joining the 
cycle to a (homologous) p-cycle which is its intersection with S. An analogous 
remark holds for the relative p-cycles of R (mod S): these can all be deformed 
into S and are therefore zero. 

A p-form ¢ induces on any surface S (supposed given by x” = 0) a p-form /¢, 
the components of which are precisely those components of @ which are not 
multiplied by dx”. The form ¢@ is known as the tangential boundary component 
of ¢ on S. The residual part of ¢, which contains the factor dx”, is known as the 
normal boundary component and is denoted by n@. Thus 


(1.14) @=teo+ nd = tot o1 A dx”. 


If p= 0, 6 = ¢; if p = m, te = 0. On S, and only there, this division into 
tangential and normal components has an invariant meaning. From the proper- 
ties of the dual operator (7) it follows that the commutation rules 


(1.15) at = Ne, «n = lx 


hold. Also we note that if ¢ is assigned on S (both ¢@ and n@), then td and 
nd are thereby determined. However, nd¢ and té¢ can be assigned independently, 
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since each component of these forms contains a normal derivative of a component 
of . Note especially that nd¢ and té¢ have (">") and (—) components respec- 
tively, making a total of (") components. Finally, we easily verify that 


(1.16) td = ds td, 


where ds is the differential operator (6) in the surface S. 

The surface S will be the carrier of the “initial” data for the Cauchy problem 
for the p-tensor B.L. equation. Since this equation is of the second order, it is 
to be expected that values of the components of ¢ and of their normal deriva- 
tives will be assigned on S. In order to express the data on the surface S in an 
invariant form, we shall use the following lemma. 


LemMA I. Let S be a surface whose equation is x™ = 0 in a sufficiently differ- 
entiable system of coordinates in M, and let o be a p-form defined in a neighbourhood 
of S. Then the specification of any one of the following sets of data is equivalent: 


rs] 
(1.17) (a) o, nd¢, td; (b) , Dud, (c) ¢, Se 


The vanishing of one implies the vanishing of the other two. 


Proof of the lemma will be given in cyclic order. First we show that knowledge 
of (a) enables us to calculate (b). If i;...%, are all less than m, then from (6) 
follows the value of 


Dn di, ot 

in terms of 
(d@) 4... tym 

and covariant derivatives of @ along the surface. Similarly, if i, = m, say 
Dads... tem 


is given by (8) in terms of data included in (a). Next, (b) clearly implies (c), 
since all components of ¢ are known on S. Finally, (c) implies that d@ is known 
on S, by (6). Since +@ and D,,«@ are obtainable from (c), 5¢ is also known on S. 
This proves the equivalence. The last statement of the lemma is evident. 

Let C be a (p + 1)-chain of M with real coefficients, and let its boundary be 
denoted by 6C. Then Stokes’s formula 


(1.18) Jas = Je 


holds. Hence, if @ is closed, its integral over a bounding cycle C vanishes. It 
follows that a closed p-form ¢ has periods f. ¢@ on p-cycles Z, which depend only 
on the homology class of Z. In a manifold with boundary [2; 8b], a closed 
b-form ¢ is derived if and only if its periods (on absolute cycles) vanish. 

If ¢ and y are two forms, the sum of whose degrees does not exceed m, we 
have 
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(1.19) d(o¢A ¥) =coA¥+(-1) GA dy, 
where p is the degree of ¢. 


Consider a subregion of M, which we may take to be the region R. Let a and 
B be two forms of degree p defined in R, and let 


(1.20) ~~ Jia A= Rade 


denote the scalar product over R. We remark that in the elliptic case (ds? > 0), 
the scalar square (a, a), is positive definite; however this no longer holds for our 
metric. However, for any two forms ¢ and y of degrees p and (p + 1) respectively, 
we have the formula of Green: 


(1.21) (d¢, V2 — (68 Ve = Je A ot, 


which follows from (18), (19), and (20). The metric being indefinite, we should 
change the sign of the surface integrals over any timelike surface. However, for 
all applications which we have in view, the surface dR will consist of spacelike 
surfaces and null cones, and the integrand will be made to vanish on the latter; 
so the formula will apply as written. From (21) we obtain by formal transforma- 
tion the extended formula of Green, namely, 


= fo Andy — VA ado +56. oh — 59 A +9), 


where ¢ and y are now of equal degree. Regarding the right-hand side of (22) 
as a linear functional of the form ¢, w being held fixed, we observe that the 
successive terms contain t¢, nd¢, té@, and n@; and that each of these expressions 
can be assigned independently of the other three. All four together are equivalent 
to any one of (17). 


2. A uniqueness property of the B.L. equation. Let R = Rs be the region 
of dependence associated with an initial surface S as described in the preceding 
section. In terms of the geodesic normal coordinate system (13), we may write 
the partial differential equation (11) in the form 


= 84.0.4 ap 9 $4,...4 
(2.1) L(¢...4) = “se to 


R( Js... Ip) CG] (ja--- Jp) 
+ G4,...% or. D6 5s.-049) HOt ty PEt) = 9, 


where t = x", i < a, 8 < m — 1; and 


R(ji.-- dp) (js1-+- Jp) 
G4... % aati ’ ae ' 
are polynomial functions of the components of the metric tensor and their first 


derivatives, which are skew-symmetric with respect to the indices 7, ... 4; 
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ji-++Jp- In this coordinate system the Cauchy problem consists of solving (2.1) 
(in general with a non-homogeneous right-hand term) given the values of the 


0 
%:,...4, and ry) P1,... 4, 
on S. 
Lemma II. The solution of the Cauchy problem for (2.1) is unique in R,. 


The argument which we use is of a standard type [1, II, p. 310]. Let P bea 
point of R with retrograde null cone Cp, and suppose ¢ and 4¢/dt are assigned 
the value zero on the section Sp of S enclosed by Cp. We have to prove that 
¢(P) is zero. Let $, denote a typical component of ¢. For Q € R, we have 


so that by the Schwarz inequality 


$1 (Q) < #(Q) f° ya = tf outer, 


where ¢ = ¢(Q) and the integration is taken along a parametric line of the 
coordinate system. 

The subscripts a and ¢ will denote partial derivatives with respect to x* and 
t, respectively. Let D, be the region enclosed by Cp, and the planes ¢ = 0, 
t = h, and Sp(h) that part of the plane t = h cut off by Cp. We have 


(2.2) a $2(Q) dA <h Ss.ay@4 f $2 dt 


h 
< nf c2dV< nf E,(h) dh 
D, 0 


where 
(2.3) Eu(h) = fs (oui + £%6ube) dA. 


This follows since the metric form gag dx* dx® in the spacelike surface S, is 
positive definite. Integrating (2.2) with respect to ¢ from 0 to h, and noting that 
E,(h) > 0, we have 


rh 
(2.4) f @¢dV< wf E,(t) dt. 
D, 0 


We next consider the identity 
(2.5) bee (b1) = 2herdiee — 2g bras dir + Bil Gy, Gay O41) 
= (ou) + (e% babu): — 2(e btadisds + Bildy sar O51), 


where B, denotes an expression linear in the components of ¢ and their first 
derivatives, multiplied by ¢,; Let a, and ¢t, be the direction numbers ds/dx-, 
ds/dt of the conormal » to a surface. On the surface Cp we have ds? = 0, whence 
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(2.6) t,? = g” aB, on Cp. 


We integrate (2.5) over D, and apply the divergence theorem. Since L(¢@,) is 
zero we have 


' oe ; 
0= Sov. (o00"t + (g "biabis)ty —_ 2(g "bia i1) Br dA + Jip, Bedl : 


On Sp (t = 0), ¢ and its derivatives are to be zero; on S, we have t, = 1, a, = 0. 
Using (2.6) we then find 


1) pe a a3 af, 2 \ 
—_ 1, \8 AaB dit 2t,g B biadvit + g t, Piades f dA 
2 a8 
+ Js4 oe + 8° babu dA 
= ‘ 1 as - pn 
= | ry g (tr Oia aris) (trou Beit) dA + E,(h) 


=- JpBeav = Ry, 


denoting the value of the last integral by R; Here we have denoted by Cp(h) 
that portion of Cp which lies between the planes ¢ = 0, t = h. On Cp, t, is 
positive. Recalling that the metric form gas dx* dx® is positive definite, we 
conclude that 


(2.7) E,(h) < Ry. 
Adding together these relations for all components @¢,;, we find 
(2.8) Zh) = >> Eyhh) <> RR = R. 


Now the quantity R defined by (2.8) is an integral over D, of terms ¢i 
$j: OF $s $12, each of which is less than half the sum of the squares of the two 
factors. The coefficients of these terms in the integral R are bounded. Hence, in 
view of (2.4) we have an inequality 

wh 
(2.9) R< K | Z(h) dh, 
0 
for some positive constant K. From (2.8) and (2.9) follows 


vh ek 
Z(h) <K | S(t) dt <K { S(t) dt 
e/0 0 


for every k > h. Integrating between 0 and h, 


ah h 
(2.10) | Z(h) dh < Kh f =(h) dh. 
0 0 


But (2.10) is manifestly false for h < 1/K, unless =(h) = 0, (0 < A < 1/K). 
Hence (hk) must vanish in this range, so that all derivatives of components of 
¢@ vanish for t < 1/K. Hence ¢ itself is zero in this interval of values of ¢. Since 
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K can be chosen to depend only on S and P and the metric, we can repeat the 
above process starting from the surface S,,x, and so on. Hence finally we obtain 
¢(P) = 0. This completes the uniqueness proof. 


3. The Riesz kernel. A detailed exposition of the solution of the Cauchy 
problem for normal hyperbolic equations in Riemannian spaces has been given 
by M. Riesz [9]. We shall construct the solution for the p-tensor B.L. equation 
using Riesz’s method. For a complete discussion of matters upon which we only 
touch, the reader is referred to Riesz’s work. 

It is required to construct a double p-tensor kernel 


a ra a \ 
(3.1) V"(x, y) = V,(x,¥) = Vice ¥) «. "™ freee dof 
having the properties: 
(a) 4,V**? = V% 
(b) V%(x,y) = 0(s*™”), 
where s is the geodesic distance from x to y, 


(c) V“(x,y) = V*(y, x), 


(d) ‘. V"(x,y) A «V*(z, y) = V***(x, y), 


where Dj, is the double conoid enclosed by the retrograde null cone of one 
argument point and the direct cone of the other, 


~ 1 
(e) V"(x, y) lyr = ye pli dx“ A... A dx"*dy" A... A dy”. 
i ee 
We proceed to construct V* by means of a series expansion in powers of s. 
Properties (c) and (d) will be established later. 
Certain facts from the theory of geodesics, which we now set forth, will be 
required. Riesz gives a detailed account of these matters. Choosing the point x» 


as origin, we denote by x'...x” a coordinate system in M. On a family of 
curves C, let ~:... Pm be the conjugate variables defined by 
dx* 
= sk gk OX 
(3.2) Pi = fur, & da’ 


where ¢ is a parameter along C. We introduce the Hamiltonian function H (xp, x), 
where 


(3.3) 2H (xo, x) = g"D: Pe, 

and note that the differential equations of the geodesics can be written 
9H , _ _ oH 

(3.4) + = ap,’ py _ ax' 


The geodesics passing through x» are determined by the quantities #{. As we 
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shall need to consider only timelike geodesics (ds* > 0), we may take s as 
proportional to the parameter ¢ on the curve. If we set o(x) = 1, then since H 
is constant along a geodesic we have 


1 - 1 — 
s = s(xo,x) = f V gat’ do = f V 2 bide da 
1 ——— —— 
-f V/ 2H de = ~/ 2H; 
0 


letting P = s* we find P = 2H. 
We now select normal coordinates &* which determine [9, p. 173] the point x: 
(3.6) f= 0,%5 = s £. 


Taking S as parameter, we have & = & along the geodesic arc. 
A number of useful identities may now be derived. Since H is constant along a 
geodesic, we have 


(3.7) 2” (0)ni m = ge” (x)ni m 


where the 9; are variables conjugate to the #*. From (3.7) follows in view of 
(3.2) the relation. 


(3.8) gu(0)E‘e* = ga(0)s* (0, E)£0 
= gu(t)s*(0, te = ga(tye't’; 


(3.5) 


hence from (3.5) 


(3.9) P(O, ) = s°(0, &) = gu(O)e'e® = galeee’. 

Thus, choosing Riemannian coordinates at x» such that gnm = 1, gx = — 1, 
k # m, g* = 0, i # k, we find for P the Lorentz distance 

(3.10) P = s* = (¢)’ — (")* -—... — )”. 


Next, denoting dP/dt‘ by P,, we have 
P, = 2g,¢’, P* = 28°, 
whence 
(3.11) PP, = 4g,4't' = 4P. 
This partial differential equation for P shows, that, as has been mentioned, the 
geodesics of zero length through x, form a characteristic surface for the B.L. 


equation. In fact, the geodesic lines are just the bicharacteristics. 
Again, let V be a differentiable function of position; then 





QU _ 91 8U _ 4 OU _ 4 dt'aU _ , dU 
(3.12) P ae" = 2é ae = 2s€ ae = 2s ds ae = 2s ds ° 
Finally, we note that 
a°P a 
ij ~~ ij tjpk OL pe 
(3.13) o aE ar 2g "S13 + 2g °E ae! 


2m + 0(s°). 
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These formulae will now be applied to the construction of the kernel. We 
shall need to calculate the Laplacian of a p-tensor ¢ multiplied by a power of 
the geodesic distance. Let F(P) be a scalar function of P. Then we have 


A(F(P)¢) = D'D,(F(P)¢) + BF(P)¢ 

= F(P)A¢ + (D'D,F(P) + 2D‘F(P)D,)¢, 
where B denotes a suitable multiplying matrix. We find 

D'D,F(P) = D‘(F'(P)P,) 
(3.15) = g°[F'(P)PPit F(P)Pa — {eh F(P)Po) 
4PF"(P) + 2mF'(P) + F’(P)sA,, 


where A, is a suitable bounded multiplying matrix. For the last term in (3.14) 
we find 


D‘F(P)D1$\),.... 49) 


te] i< 
r)| Pf DU 4s... 49) - P Dy { 0) Oy... dosh — 


F'(P)| 2s FSi! + sA — 


where A, is a second matrix with bounded elements near xp. 

In the preceding equations, ¢ is to be understood as a vector with (") com- 
ponents, while A,, A; and B are square matrices of order ("). Continuing with 
this notation, we have 


A(F(P)¢) = F(P)A¢ + [4PF"(P) + 2mF'(P)]¢ 


(3.14) 


(3.16) 


ll 


(3.17) de 
+ sF(P)A¢ + 4F(P)s =, 


where A = A; + Az is also bounded, the bound depending solely upon the 
bounds of components of the metric tensor and their first derivatives. 

Following Riesz, we set 

a—m-+ 2k a—m o~ 2k 
ra < 5 V, 5 Ss V,(0, §) 

L V 0, = wae. > 
(3.18) , 8) = 2 lak) ~ Kula) alate + 2k) 
and determine the double (matrix) forms V; successively so that (a) shall hold. 
Since P = s?, we have 








a+m+2k+2 
ae As V, 


L Kula + 2) Lm (a + 2k) 


(3.19) © _ ge (+ ’ ? . 
2, Hala + 2k)\° AV; + (@ — m + 2k + 2)(a + 2k) Ve 





AV**? = 





+ 2(a — m + 2k + 2)sAV, + 2(a - m+ 2k + 2)s 44s) 
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——— > a s 
ds 


s af AVes a—= + m+3| 2.4%, 
Ln(a + 2k)  Lm(a + 2k + 2) 


+ (a + 2k)V; + 2sA v.]). 


We now choose 


(3.20) Lnla + 2) = (a+ 2 — m) L, (a), 

and determine the V, by the recurrent system of differential equations 
7 

(3.21) 2s os + (2k + 2sA)V; + AVi-1 = 0, 

with 

(3.22) V_.,=0, V.(0,0) = J, 


where J is the unit matrix of order ("). The expression (3.19) now becomes 


a—m+ 2k 
a . Ss V;, 


a+2 —— ioe 
AV" = Kula+2) $s Lala + 2k) 


which is equal to V* if we choose 
(3.23) K,, (a + 2) = aK,,(a). 
From (3.20) and (3.23) we have 
K,(a) = K2*r (:); Ln(a) = L2* 7 (2+3-=), 
and we choose K, L so that (e) is satisfied: 


(3.24) Hy (a, k) = gi gett (2) r (2 + 2k + 2—- m). 





In these details we follow the scalar treatment. 
We have now to prove the convergence of the series for V*. The relation 
(3.21) can be written 


(3.25) ois Vz) +As* Vz = — 4415" ". 
To solve (3.25) consider the adjoint system 
(3.26) fyea'y 

e ds —_— ’ 


and let Y be a (nonsingular) matrix solution of (3.26) with | Y(0)| = 1. It 
follows from the Lyapunov relation 


Y(s) = Y(O) exp | fe A as| 


that Y(s) is nonsingular, hence Y-'(s) exists and is bounded for |s| < 1. Since 
s* V, vanishes for s = 0, we find as solution of (3.25) 
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(3.27) *V,=—-—} ys) f V(c)AV,-3(c)o" “deo. 
0 
With the vector and matrix norm 
XI] = V zx, 


(sum over all components), we have 


3.28) IMAI < BLY CDI J 1YGDI| [l4Ve-a(@)lo* ae. 


Let us assume that for suitable constants r and k and for |s| < «, 


(3.29) [|Ve-s|| < x(1 - e)-m, p= 2 |e; 


Then [5; 6], 
—2z-—2 
[|AVe-a|| < KK,(1 _ e) : 


and from (3.28) we have 
s k—1 rk —2z—1 
TT TAT te ( _¢ do » MK Ki -¢) 
(3.30) s*IIVil] < MEK, fey < st MEE | 


where A and M are independent of k. Since an estimate (3.29) holds for V4 it 
holds for V; and hence for all V,. It follows from (3.30) that the series 2s" V, 
converges geometrically for s sufficiently small. From (3.24) it follows that the 
series for V* converges for all values of a if s # 0, is sufficiently small; and 
uniformly in any right half plane of a. Hence V* (0, £) is defined as an analytic 
function of a, and regular for R(a) > m and & € R [ef 9; 5]. 

Exactly as in [9, p. 186] it follows that an estimate 





—a| 


(3.31) Iv", 911 < cs*-*|(2) | jas 


holds, where R(a) > m, and C, q are independent of a. 

We remark that the present construction follows that of Riesz very closely; 
the only additional features being that it was necessary to show that the matrix 
A was bounded and to solve a system of recursive equations for the V,. For the 
non-analytic case Riesz has given a method of construction for the kernel, and 
the same formal considerations apply to the present p-tensor kernel. Since 
the details are lengthy we omit them. 


4. Construction of the solution. We continue the adaptation of Riesz’ 
method for the p-tensor B.L. equation. We now define the Riesz fractional 
potential which plays the central role in the solution of the Cauchy problem. Let 


(4.1) I“o(x) = (¢(y), V(x, ¥)) ne = / o(y) A «V"(x, y), R(a) > m — 2. 


Here D*% is the conical region bounded by C, and S, and ¢ is a p-form defined in 
R and assumed sufficiently differentiable for the prupose in hand. Riesz has 
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proved certain relations in the scalar case, which as we now show, can be extend- 
ed to the present p-tensor treatment. If @ is of the class C[4(m — 1)], then 


(4.2) Pr o= ¢, 
and in general, 

(4.3) AI*o = I*~*4, 
and 

(4.4) I*I°¢ = I*** 9. 


Riesz has also given a method of analytic continuation of J*@ into the left 
half-plane R(a) < m — 2. 
First we shall prove (4.2). A typical component of J*¢ is 


(I$) innte = Ss > ar | ee | A 
(4 5) D, Ri<...<Bp 


. : 
= . = i‘. oe....te FF Ce, 0) Vix... tp; Ra--typ(%) ¥) CVn, 


where v*(x, y) is a double p-form regular on the cone D% and such that 


(4.6) v"(x,x) = Vo(x,x) = J,. 
The recurrence relation (3.25) for V» can be written 
ot + A Ve = 0, 
ds 


whence it follows that the off-diagonal elements of Vo, hence also those of x 
(x, y), are O(s). The typical component (4.5) then breaks up into a sum of 
integrals of the type studied by Riesz. We may assume normal coordinates 
lending to the Lorentz distance (3.10). Using this coordinate system, it has 
been shown [4, 9] that an expression of the form 


re So: Hata.) 90) 2Vm 
depends only on the value of ¢ at the point x, and is in fact equal to ¢(x). This 
result implies that as a — 0, all terms of (4.5) save the diagonal term tend to 
zero due to the factor O(s), and that the remaining term is precisely @,,... :,(x). 
This establishes (4.2). 
We now introduce the formula of Green (1.22) for the region D%, and set 
yy = V***(x, y). We obtain 


(Ad, yet? — (¢, av 


vy he (V*** A add — OA ad VV? + 5V"" A 2g — 56 A # V™"). 


Supposing that R(a) > m, the surface integrals over the cone vanish because of 
the factor s*+*-™ contained in V*+*, and we are left with the surface integrals 
over S*. Supposing further that A¢@ = p and replacing AV**+* by V*, we find 
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(4.7) I*¢@ = r)— f cy" A «do — oA edV*** 4. gy°** A #6—5dA oV***), 
To facilitate the study of (4.7), we write 


Iie) = — f VA ado 


ai Jul, a) + f orsav 


J3(¢,a) = — f 5V" A #¢, 


Juda) = + f 360 oV*; 
then (4.7) may be written 


4 
(4.9) I*o@ = I***p + » J4(, a + 2). 
We now formulate our theorem. 


THeorREM I. Let S be a spacelike surface of the class C[4(m + 3)], and let 
R be its region of exclusive dependence. Let ¢ be a p-form of the class C{[4(m +- 3)} 
defined in a neighbourhood of S (or let tt, nt, ndt, t8t be given on S); and let p be a 
p-form of class C[4(m + 4)] in R. Then there exists in R a unique p-form @ such 
that 


(4.10) Ad=p 
in R, and 
(4.11) t@¢@=t, no=nt, ndd=ndt, tid = tet 


on S. This form is given sufficiently close to S by 


: 4 
(4.12) @= I’p+ Du J,(é, 2). 


The main steps of the proof are as follows. It must first be shown that (4.9) 
can be continued analytically to a = 0, so that (4.12) has a meaning. Then the 
conditions (4.10) and (4.11) are to be verified. The uniqueness has been estab- 
lished in Lemma II, so the proof will then be complete. 

The analytical continuation has been carried through by Riesz [9] and by 
Fremberg [4, pp. 19-53]; since the essentials of the matter are in no way altered 
in the present case we shall not repeat the details. We remark that the continua- 
tion is possible if the forms p, — and the surface S have the differentiability 
properties stated in the theorem. It also follows that (4.12) is then continuous 
together with its first and second derivatives. 

At a = 0, the function H,,(a, 0) has a simple pole. Fremberg has also proved 
that the surface potentials can be continued toa = 0, leading to finite expressions 


apart from the factor 1/H,,(a, 0). It follows that the expressions J, vanish for 
a= 0: 
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(4.13) J,(é, 0) = 0. 


We must now establish a number of relations involving the B.L. operator 
and the integral /*. First we make the important observation that in calculating 
derivatives with respect to x of J* and the J,, we may disregard the variation at 
the limits of integration [9, pp. 67-70]. The justification of this is that for R(a) 
> m + k, the kernel V* is zero, together with its derivatives up to order k, 
on the cone D,. It follows by analytic continuation that for all a we may calculate 
derivatives by formal differentiation under the integral sign. In view of condition 
(a) of §3, it therefore follows that (4.3) holds whenever it has meaning. In par- 


ticular, if the continuation is valid to a = — 2, we have 

(4.14) Ag = I*¢, 

in view of (4.2) and (4.3). Similarly, the relations 

(4.15) AJ,(¢,a) = J,(¢,a — 2), » = 1,2,3,4 
hold. From (4.13) we see that 

(4.16) AJ,(é, 2) = 0, p = 1,2, 3,4 


The relation (4.4) follows as in the scalar case from the formula (d) of §3 
which may be written 


(4.17) (V", V*) = Jye¥*2) A «V*(y, 2) = V°*(x, y), 
u u 


and which we prove in essentially the same way. We apply Green’s theorem to 
D; for R(a) > m, R(8) > m, and, noting that the surface terms vanish, we 
have 


a+2n 78 . os ya+2n = An 7B+2n . 
(V (x, 2), J 0 )) pe (J (x, 2), At 2) ps 


= (ATV***"(x, z),V°r*(y, 2) pF = (V“(x, 2), very, z) )p* 
n=0,1,2,..., and, letting x, 8 0 we find [9, p. 196] 
v(x, y) = V'"(y, x). 
It follows that the coefficient forms V; are symmetric, hence 
(4.18) V"(x,¥) = V“(y, x), 
so that condition (c) of §3 is satisfied. 
Consider the expression 
(v*, v*) 


2- ve, R(a) > m, R(B) > m. 
D, 


A typical component g(a) of this double form is an analytic function of a which 
satisfies the estimate 


(4.19) lg(a)| < (2) “er 
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for R(a) > m. Also g(a) is zero for a = 2n (nm = 0, 1, 2...) by the preceding 
formulae. We may now apply the following theorem [7, Abs. III, problem 298, 
pp. 142, 327]: If two analytic functions g(z), h(z) are regular for R(z) > a, 
and h(z) # 0 for R(z) > a, and 


lg(z)| < |h(z)|, 


then, denoting by z, the zeros of g(z) with R(z,) > a, |z,| > 1, the divergence of 
the series 


> R(1/2,) 


implies that g(z) is identically zero. Taking h(a) equal to the expression in the 
absolute bars on the right-hand side of (4.19), and noting that the series for the 
zeros of g(a) is the divergent harmonic series, we conclude that g(a) = 0. 
Hence (4.17) holds for R(a) > m, R(8) > m, and therefore in general. 

Using (4.17) we find by inverting a certain double integral that 


(4.20) I*J,(¢, 8B) = J,(o,a + B), yp = 1,2, 3,4; 


for R(a), R(8) sufficiently large, and therefore in general. 

The kernel V* has certain additional properties which follow from the fact 
that if an operator T commutes with A it commutes with /*. We first establish 
this assertion. Let T be a (linear) operator such that TA = AT. We wish to prove 
that 


(4.21) TI*@ = I°T¢. 


Now (4.21) is certainly true for a = 0, by (4.2). Using (4.3), we can prove 
that it holds also for a = 2n (nm = 0, 1, 2,...). For, 


TI’¢ = PATI’¢ = PTA ¢ = IT 4, 


since 7 and A commute, and Ais inverse to J?. By iteration of this, we find that 
(4.20) holds for a = 2n. Consider now any typical component of 


(4.22) TI*o — I°T¢, 

this component is an analytic function of a which admits the estimate (4.19) 
and vanishes for a = 2n. By the theorem we have just quoted above, we must 
have (4.22) identically zero. This proves that T commutes with /*. 


Now A commutes with », d, and 6. Applying our remark to the + operator, we 
have 


*(O(y), Vole, ¥))oe = (#609), Vin-v(, ¥))oe = (—) (60), * Vin») o» 
for all ¢. Taking the dual with respect to x, we find 
(4.23) V3(x,y¥) = ++ Vin—v(%, y). 
Similarly, from dJ* = J*d we find that 


(4.24) d Vi(x, y) = 6 Vius (x, y). 
z uv 
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The relation 6J* = J*6 yields nothing additional. Note that (4.23) and (4.24) 
are essentially the same as relations satisfied by de Rham’s Green’s form [8a] 
and for the same reason. 

From these relations it follows that the “single layer’ potentials J, and J, 
are dual in the sense that 


(4.25) Ji(*@, a) a «J 4(¢, a). 
Similarly, the “double layer’’ potentials J; and J; are related by the equation 
(4.26) J2(*o, a) = «J3(¢, a). 


Since derivatives may be calculated formally for the potentials J,, it follows 
from (4.8) and (4.24) that 


(4.27) 5J2(¢,a) = 0, dJ3(¢,a) = 0. 


In order to show that our solution (4.12) satisfies the conditions (4.11), we 
must examine the potentials J* and J,(£, 2) as x tends to the surface S. First let 
us consider the volume potential J*. Just as in the scalar case it follows that 
I*p + 0 since the conoid D% becomes infinitesimal as x tends to S. The first 
derivatives of J*p also tend to zero. 

As for the surface potentials, it will be sufficient, in view of (4.25) and (4.26), 
to consider J, and J. We shall reduce these expressions to integrals of the type 
studied by Fremberg, and apply his results. We may assume that S has the 
Lorentz form (3.10); it follows that as x tends to S, s becomes small of the same 
order as &*. We may discard all terms of V* except the first, noting that these 
terms, together with their first derivatives, become negligible in comparison 
with the first term. 

Considering first J,, we have 


” ean 
~—- | ——I, A eft, 
Jule, @) hia pm 
whence, on S, 


(Si) a... = — Jim (T4,...% A «dt)1. m1 dx" A... Aa 
In the integrand we have 


| 7 A #&) 1... (m—1) = ( ies 1)? (dE) 4,... tym 


Hence the components of ¢/, are simple layer potentials of Fremberg’s type, so 
that [4, p. 44], 


d p 
tJ, = 0, dn tlt) ts. te = ( — 1) (dE) 4,... tom: 


Here n denotes the inward normal to the surface of the conoid. The integrand is 
zero if, say, 1p = m, so that 


d 
nJ, = 0, in ty) = (), 
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~I 


or 


We now calculate 


(mT) to... tom = Tagen. tyme OT 1h 5...059) 
mg tg On J3,...49) Since nJ, = 0 
= (-1)’ Pe J iG... ty) 
= (dE) ¢,... tom: 
Finally, we have 
(t6J)4,...4-. = — Tus, ited | 


If the summation index i is less than m, then none of the j, are equal to m, 
hence we get a tangential derivative of a component of t/J;, which yields zero. 
If i = m, we get the normal derivative of a component of nJ;, which is again 
zero. Altogether, we have 


(4.28) tJ; = 0, nJ = 0, nd J, = ndé, tb; = (), 
A similar analysis holds for J,. We find 


AY 
Jx(—,a) ~ Juz A +d F(a) L, 


"ts = 
- ft: A oly fo) mn ai. 


Taking a typical normal component of J2, we find that the integrand is of the 
form 
eer 


err ‘Ox Hala)’ 


6, Cx occ heeg XM 


which shows that components of mJ; are tangential derivatives of single layer 
potentials. As we have seen above, these are zero. Hence nJ2 = 0. For tJ, we 
find 


r = ° aE ieee 
(« A #Ylotinnt A tf })..we = §4...4 In Hela)’ 


which leads to a double layer potential of Fremberg’s type for (¢/2);,...«,. 
From his results [4, p. 48] we conclude that 


d 
tJ: = t, tJ, = 0. 
We note from (4.27) that /é/2 is zero. Again, 


(dJ2) ,...tgm — Ps... ant 0 Jats... 


For 7 = m, we obtain a normal derivative of a component of ¢/2, which is zero. 
For j7 # m we obtain a tangential derivative of a component of nJ2, which also 
yields zero. Collecting our results, we have 
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(4.29) tJ = té, nJs = 0, nd J» = 0, tJ = (. 


The corresponding results for J; and J, follow from (4.25) and (4.26). 
The solution (4.12) may now be verified. From (4.14) and (4.16) we have 


Ag = A{I’p +X J, (&,2)} =, 


so that (4.10) holds. To verify (4.11) we note that the volume potential and its 
first derivatives contribute nothing on S. Applying (4.28), (4.29) and their 
duals, the result (4.11) follows. 


The proof can be used to show that the kernel V* (x, y) is uniquely determined 
independently of the coordinate system in which it was constructed. The 
volume potential J*p is uniquely determined, hence so is V*. It follows from 
(4.17) that V™ is unique for m = 1, 2, 3,...; and the uniqueness theorem for 
analytic functions quoted above shows that V* is unique for all values of a. 
The kernel V* is defined in a small region only, but repetition of the integration 
process shows that the solution can be extended throughout R. 


5. The harmonic field equations. We consider the system of first order 
equations for the components of the p-tensor @,,...;,; 
(S.1) do=p, i¢=¢0 
and the problem of solving this system under conditions sufficient to determine 
the solution uniquely. When p and o are replaced by zero in (5.1) we have the 
equations of harmonic fields. The solution of our problem will utilize Theorem I. 
Let S be a spacelike surface and R its region of dependence, as before. Obvious 
necessary conditions that (1) should possess a solution are that p be a derived 
form, and o the dual of a derived form. Hence we impose the conditions 


dp = 0, f p = 0, 
(5.2) Zp+. 


bo = 0, J xo = 0 


for all cycles 2,41, Z,—p+1 Of R. 

Since (5.1) is a first order system, it is to be expected that the assignment of 
t@ and n@ on S will determine the solution. That this is indeed true follows from 
Lemma II. For, the homogeneous system 


(5.3) d@ = 0, 66 = O0inR, 6 =O0o0nS 


has only @ = 0 as solution. To show this, we observe that Ag = 0, ¢#¢ = 0, 
nd = 0, ndd = 0 and téd = 0. Hence ¢ = 0 in R. 

Returning to the system (5.1) we observe that if ¢ = — on S, then ~ must 
satisfy the conditions 
(5.4) tdt = tp, nit = no, 


since the corresponding component equations contain only derivatives in 


| 
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directions lying in S. A p-form ¢ defined on S which satisfies (5.4) will be called 
admissible with respect to (5.1), or, more briefly, admissible. 

As a simpler example of the method which we shall use to solve (5.1), we 
consider the homogeneous equations 


(5.5) dé = 0, 86 =0 
with the boundary condition ¢ = & € C [4(m-+ 1)] on S, where ¢ is admissible 


with respect to (5.5). Let @ be the solution of A¢ = 0, such that t@ = éf, no = 
nt, nd@ = 0, nid = 0 on S. Thus 


(5.6) @ = J2(&,2) + Js(€, 2). 


It remains to be shown that ¢ is closed and coclosed. But td@ = td§t = 0, ndé = 
0, A(d@) = d Ad = 0, and nd(dd) = 0 since d.d@¢ is zero. Finally, we have 


ti(d@) = — tdid = — dtid = 0, 


since 6 = 0 on S. Hence, by Lemma II, dé = 0. Reasoning exactly dual to the 
preceding shows that 5¢ = 0 in R. Hence ¢ is the unique harmonic field which 
assumes the admissible value ~ on S. 

Since the equations (5.1) are non-homogeneous, we shall need to treat them 
in a less direct fashion. We divide the problem into two parts. Consider first the 
problem of soiving the system 


(5.7) dod = p, 56 = 0 

where 

(5.8) to = tt, nd =0 

on S, and ¢é is admissible with respect to (5.7). That is, 

(5.9) tdt = tp, 

on S. 

Assume for the moment that a solution ¢ exists. Then we have from (5.7), 

dso = 0, 

and 


f o= f «od = 0 


from (5.8), where z,_, is any (absolute) cycle of R, and 2*,_, is a homologous 
cycle of S. The cycle z’,_, homologous to 2,_, always exists (see §1). It follows 
that +¢ is a derived form in R: 


(5.10) @ = 6X. 
We now obtain 
diX = p, 
which suggests that we attempt to solve (5.7) by solving 


(5.11) AX = —p 
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under conditions which ensure that dX = 0. This is a kind of “gauge’’ condition. 
It is not difficult to write down the auxiliary conditions which X must satisfy 
on S. They are 


(5.12) iX = 0, nX = 0, ndX = 0, HX = to = Kk. 


These conditions ensure that (5.8) are satisfied. To prove that X is closed, note 
that 


A(dX) = d(AX) = — dp = 0, 
that dX = 0, ndX = 0, and that nd . dX is zero. Finally, we have 
tédX = — idiX + tp = — 1dbX + tdt 
td(é — 6X) = d tg — 6X) =0 


in view of (5.9) and (5.12). Hence dX = 0 in R. 
From Theorem I we have 


X= -Te+f EA «V3ti 
Ss? 


’ 
a=(0 





whence from (5.10) 


. 
ll 





—sFo+s fen «VSti 


— Tip + Ja(é, 2), 
in view of (4.8), (4.24) and the permissibility of differentiating under the 
integral sign. 
The dual problem, in which 
(5.14) d¢@ = 0, 6¢ = a, 


and the boundary conditions are easily written down, can be solved in the 
same way, or by use of (4.23) and (4.26). The solution is 


(5.15) @= — Ido + J3(E, 2). 


Summing up these results, we have 


a=0 


(5.13) 


THEOREM I]. Let S be a spacelike surface of class C[}(m + 3)], and p, ¢ 
forms of class C|4(m + 3)], of degrees (pb + 1) and (p — 1) respectively, which 
satisfy (5.2). Let & C[}(m + 1)] be a p-form admissible with respect to (5.1). 
Then there exists a unique solution of (5.1) defined in R, given by 


(5.16) = — (8p + do) + Ja(t, 2) + Ja(E, 2). 


We note that if p and o are harmonic fields, the first term of (5.16) vanishes, 
and the solution is formally the same as (5.6). However ~ must satisfy the 
conditions of admissibility (5.4). 


6. Concluding remarks. We consider certain special cases of the two theo- 
rems which result when further restrictions are imposed on the manifold M or 
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the metric. A comparison of the present results with the boundary value 
theorems for p-tensors in the elliptic case will also be made. 

Under certain circumstances the solutions discussed in Theorems I and I] 
may be valid in a wider region than has been described. To begin with, we have 
restricted ourselves to one “‘side”’ of the initial surface and it is clear that the 
solution is valid both for ¢ > 0 and ¢ < 0 in the appropriate regions R+, R-. 
The solution is also determined in the region R; associated with any surface S, 
which is spacelike and lies in R. Suppose for example that M is the product of 
an (m — 1)-dimensional closed manifold and a line (— © <t < + ©); with 
metric 
ds* = dt” — gag(t) dx" dx’ 


where gas(t) gives rise to a positive definite line element depending on ¢ as a 
parameter. Let S be the surface ¢ = 0: the solution is then determined throughout 
the entire manifold. If the g.s(t) are independent of ¢, the B.L. equation is 
separable: 

ou = Aso 
where A, is the B.L. operator for S. Setting ¢ = e®'u, we would have 


Au +r*u =0 
just as in the scalar case. 

If the metric is of certain restricted types, the B.L. equation will separate 
into (") independent component equations. This is the case for the Lorentz 
metric 

ds* = dt’ — (dx')’ —... — (d&x”™")’. 


Each component of ¢ then satisfies a wave equation 
ou = Zz 92,0. 


The Riesz p-tensor kernel reduces to the first term of the series: 


ya in , aa t; tp * f, ty 
V;(x, y) Hala) ae, ™ Ass. Oh. 1 AD’, 


for which the relations (4.23) and (4.24) are easily verified. This case can 
therefore be treated by scalar methods. 
Of particular interest is the case m = 4, p = 2, since the equations 


do = p, 5¢ = 0 


are then, in the Lorentz metric, equivalent to Maxwell’s equations, with u = « = 
1. In a general metric it would be natural to retain this invariant formulation. 
The two-tensor ¢,,;, is the electromagnetic field tensor F;,;, and the dual of the 
three-form p is the charge-current vector. Theorem II may therefore be inter- 
preted as a characterization of electromagnetic fields in a general relativity 
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metric. The conservation of charge current is expressed by the fact that p is a 
derived form, in other words that 

fone 
for all cycles 23. 


It is interesting to compare our results with the boundary value theorems for 
p-tensors in the elliptic case, and with the classical theory of differential equations 
(p = O or p = m). We consider first the B.L. equation. The result in the elliptic 
case is that the solution exists and is unique if #@ and m@ are given on the boun- 
dary. This result has, however, only been proved under certain topological 
restrictions [3]. In the hyperbolic case the initial surface does not bound the 
domain of determination of the solution, but the additional data md@ and té¢ 
are needed to determine the solution uniquely. If p = 0, n@ is zero automatic a'ly, 
and so is té¢. Hence the data reduce to ‘@ = ¢ and 


in the hyperbolic case, and to ¢ in the elliptic case. A similar remark holds if 
p = m. 

For the harmonic field equations the elliptic boundary value problem is 
uniquely solvable given either the normal boundary component and the absolute 
periods, or the tangential boundary components together with the relative 
periods. In the hyperbolic case both #¢@ and n@ are to be prescribed on the 
(smaller) initial surface. The region R in which the hyperbolic solutions are 
defined has no relative cycles (modulo the initial surface) and the periods on 
all absolute cycles are fixed by the given data t@ and n@ on S since all absolute 
cycles of R are homologous to cycles of S. In the hyperbolic problem, therefore, 
period conditions are not needed explicitly. 
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PARTITIONING AN ARITHMETIC INTERVAL 


WILLIAM GUSTIN 


The purpose of this paper is to characterize all ways in which an initial 
interval of natural numbers can be partitioned into a unique arithmetic sum of 
certain of its subsets. 


Preliminaries. We set forth below certain explanations, conventions, and 
definitions pertinent to our subject. 

Numbers and sets. By number we mean a natural (that is to say, an arithmetic 
or finite ordinal) number 0, 1, 2, . . .. Every set of numbers is required to contain 
0. A set containing a number other than 0 will be called proper. The improper set 
comprising solely 0 will be written 0. The set consisting of all numbers dx as x 
ranges over a set X is denoted by dX. 

Intervals. An (initial arithmetic) interval is a set of numbers containing 
every predecessor of each number in it. A finite interval containing d numbers 
thus consists of the first d numbers r < d where d > 1. An infinite interval 
contains all natural numbers, that is, all numbers r < w where w is the first 
transfinite ordinal. Denote by J(d) the interval consisting of all numbers r < d 
where 1 < d < w. J(d) is finite or infinite according as d < w or d = w, proper 
or improper according as d > 1 or d = 1. For this reason we call an ordinal d 
proper if 1 <d <a. 

Proper sequence. A finite (terminating) or infinite (non-terminating) numeri- 
cally indexed! sequence d,, dz, d3, . . . of ordinals will be called proper if not 
only is each ordina! in it proper but its ordinal product is also. Thus 1 < d,, 
< w for each index m and 1 < d,; dz d;. . . < w. Consequently, if there is an 
index m such that d, = w, the sequence terminates at this index m. 

Partitions. We shall say that a set X is partitioned into (a unique arithmetic 
sum of) finitely or infinitely many numerically indexed? sets X* if every number 
x in X can be uniquely expressed in the form x = }>x* where the unique co- 
ordinate x* of x belongs to X* for each index yu; and if furthermore every number 
x of this form belongs to X. To indicate this we write X = }/X*. It is easily 
seen that order and grouping of the terms X* in such a sum are not significant 
and that any two differently indexed terms X* have 0, and only 0, in common. 

Division algorithm. As an apposite example of partitioning an interval 
consider the division algorithm. Let d,, d. be a proper sequence of two ordinals. 
According to the division algorithm every number r < d;d, can be uniquely 
expressed in the form 


Received April 14, 1951; in revised form August 14, 1952. This work was done under a 
grant from ONR. 

1Subscript indices indicate that the order of indexing is significant. 

*Superscript indices indicate that the order of indexing is not significant. 
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r=n+dro, 


where r; < d,; and rz < dz; and furthermore every number r of this form is 
< d,; ds. That is, 
I(d, d2) = I(d;) + d; I(d:). 


Statement of results. Our first result concerning partitions of an interval 
merely extends the division algorithm to any proper sequence and provides 
a recipe for constructing interval partitions. Our second result asserts that every 
interval partition can be constructed, and in a certain sense uniquely, by 
following this recipe. 


THEOREM I. Let dj, de, ds, ... be a proper sequence of ordinals. Then 
I(d; dz ds...) = I(d;) + dy I(d2) + di dz (ds) +... 


Coarser repartitions of I (d, dz d;...) can be formed from this partition by grouping 
the above terms together in arbitrary fashion. 


TuHeoreM II. Let I = >\X* be a partition of a proper interval I by indexed 
sets X*. There then exists a unique proper sequence of ordinals d,, d2, ds, . . . gene- 
rating the following finer repartition of I: 


I= I(d; dz dz...) - I(d;) + dy I(dz) + di dz I(d3) +..., 


the sets X* of the given partition of I being formed by appropriately grouping 
together the terms of this finer repartition of I, with consecutive terms allocated to 
different X* sets. 


Proofs. We shall establish the results stated above somewhat formally 
by using three lemmas. The first lemma is used in proving the first theorem, all 
three lemmas in proving the second. 


LemMA 1. Let there be given a proper sequence of ordinals d,, and sequences 
of sets Ay, and Xm, all sequences of the same length, such that: 


fle (m terminal), 
An = 
Am + dm Xm+1 (m non-terminal). 
Then 
Xi =A,+d,Az2+d,d2A3+.... 


Proof. For each x; in X; the sequences x, in Xp, dm in A» are uniquely 
determined by recursion as follows: 


Gn (m terminal), 
tn = 
Om + dm Xm+1 (m non-terminal), 
so for each m, 


Xi = 0, + dj aeo+... +d; de... dmi1 Xm. 


— tr 


bam | 





) 
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If m is terminal x, = ad». If, on the other hand, the sequence d,, does not ter- 
minate, then by choosing m so large that d; dz. . . dy, > x; we find x, = 0 and 
a, = 0 for all » > m. Therefore, whether the sequences terminate or not, 


X= 4 +d,4.+d,d2a3;+..., 
this representation being unique, as was to be shown. 


Proof of Theorem 1. Let d,, be a proper sequence of ordinals, terminating or 
not. According to the division algorithm 


T (dm Ami +» -) = I(dm) + dm I(dm4i - « «) 


for all nonterminal m; so the hypotheses of Lemma 1 are satisfied by taking 
Am = I(dm) and Xm = I(dm dnyi . . .), Whereupon 


I(d; dy dz...) = I(di) + dy I(d2) + d, dz I(ds) +..., 
as was to be shown. 


LemMA 2. Let 1=A+B bea partition of a proper interval I and d a proper 
ordinal such that A contains I(d) but not d. There then exist: an interval I with I 
= 0 if d = w, and sets A and B such that 


(1) A = I(d) + dA, 
(2) B= dB, 
(3) I =A+8B. 


Proof. lif d = w it is evident that A = J = I(w) and B = 0, so the lemma is 
verified by taking I, A, B all equal to 0. Furthermore, propositions (1) and (2) 
imply the remainder of the lemma. For since A + B is a unique sum, A + B is 
also a unique sum. Define J = A + B. Thus I satisfies (3) and in addition 


I = I(d) + di. 


From this it follows that J is an interval. For if u < v with v in J, then du < dv 
with dv in J, whence du is in J and hence x is in J. 

It therefore remains to prove (1) and (2) for d < w. Now (1) and (2) together 
are equivalent to affirming for each g = 0, 1, 2, . . . the following proposition. 

$,: af r < dand r + dq belongs to A or B, then p + dq is in A for all p < din 
case r + dq is in A, and r = 0 in case r + dq is in B. 

We establish these propositions $, by induction. Obviously > is true. 
Assume then that §, holds for all » < g with g > 1: to prove $,. We do this by 
first establishing from the induction hypothesis the following weaker prop- 
osition p, obtained from $, by putting r = 0: 

Dy: p + dqis in A for all p < d in case dq is in A. 
Then from p, and the induction hypothesis we prove §,. 

To prove ),, let dg be in A: we are to show that p + dg is in A for all p < d. 
Because g > 1, d < dq and dq is in J, sod is in J. Since A contains J(d) but 
not d, the B-coordinate of d is > 0 and hence > d, so must equal d. Therefore 
d is in B. This, together with dg in A, shows that dg + d belongs to A + B = J, 
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Consequently p + dq, being < dg + d, is also in J. Let p + dg have coordinates 
a in A, } in B. We wish to show that p + dg is in A, or, what amounts to the 
same, that b = 0. Suppose, to the contrary, that b > 0. Write 65 in the form 
b = r+ dB with r < d. Hence 8 < g. If 8 = g, then r > 0, for otherwise we 
would have dg = 6 in B as well as dg in A. The number dg + d with coordinates 
dq in A, d in B thus has alternative coordinates d — r in A, r + dq = b in B; 
which is impossible. On the other hand, if 8 < g, then r = 0 by the induction 
hypothesis, so ) = d8. Thus a has the form a = p + da where a + 8 = g and 
0 < a < q since 0 < 8 < g. Therefore da is in A by the induction hypothesis: 
’ so dg in A has positive coordinates da in A, d8 in B; which is impossible. This 
dilemma proves our contention that p + dg is in A. Thus ), is proved. 

To prove $,, let r < d and r + dg belong to A or B: we are to show that 
p + dq is in A for all p < d in case r + dg is in A, and that r = 0 in case r + dg 
is in B. Now r + dq is in I, so dq is in I also. Let dg have coordinates a in A, 
in B. Then 6 has the form 6 = d8 with 8 < gq: obviously in case 6 = dg, by the 
induction hypothesis in case b < dg. Hence a has the form a = da with a < g. 
Consequently r + da is in A: by ), in case a = g, by the induction hypothesis 
in case a < g. The number r + dg then has the coordinates r + da in A, d@ in 
B. Therefore if r + dg is in A, d8 = 0 so dq = da is in A; whence p + dg is in 
A for all p < d by p,. And if r + dg is in B, r + da = 0; whence r = 0. This 
completes the proof of $, and hence, by induction, the proof of the lemma. 

LemMa 3. Let I = >(X* be a partition of a proper interval I. There then 
exist: a proper ordinal d, an interval I with I = 0 if d = w, an index a with d not 
in X*, and sets X*, all of these unique, such that 
(4) X* = HI (d) + dX*, 

(5) I =X, 
where 6 is Kronecker’s delta: * = 1 or 0 according as u = a or not. 

Proof. Since I ¥ 0, 1 is in J. Therefore 1 = }°x* with x in X*. Clearly a 
unique index a exists such that x* = 1 and x* = 0 for all remaining indices 
8 ~ a. Define A = X* and B = }-X*. Then J = A + B. Let d be the smallest 
ordinal not in A = X*. Since 1 is in A, d is a proper ordinal. Now A contains 
I(d) but not d, so the hypothesis of Lemma 2 is satisfied, and thus the conclusion, 
whereby the unique sets J, A, B are furnished. Therefore J = 0 if d = w. Define 
X« = A; then by (1), 

X = I(d) + dX«. 
By (2), d divides every number in B and hence every number in each X*, so 


X® is of the form he 
X® = dX8, 
with B = > X*. Evidently (4) summarizes the two formulae above. As for (5), 
T=A+B= X*+ 5 X=> 
by virtue of (3). This completes proof of the lemma. 
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Proof of Theorem 11. Let I = > X* be a partition of a proper interval J. 
By recursive use of Lemma 3 we construct: a proper sequence of ordinals d,,, 
a sequence of indices a», consecutive ones being different, and sequences of sets 
X* with Xt = X*, all these sequences unique and of the same length, such that 


a" I(dn) (m terminal), 
xXi= 

8°" (dm) + dn Xii41 (m non-terminal). 
Therefore, by Lemma 1, 


t = 8°" I (d;) oa 5”""d; I (dz) a 5”""d; d, I (ds) + eee, 
whereupon 


IT = > X* = I(d:) + ds (dz) + ds da T(ds) +... , 


as was to be shown. 

This finer repartition of the given partition of J has evidently been constructed 
in a unique fashion; it will be called the resolution of the given partition. Any two 
consecutive terms of this resolution lie in different X* sets, or, as we shall say, 
are separated. 


Number of partitions. Consider for given positive numbers m and n a 
partition of the interval I(m) into m sets: 


I(n) = Xi+...+ Xm. 


Call the ordered sequence of sets X;, . . . , Xm an ordered m-partition of I(m), 
and let ,,(m) be the number of such partitions. 

The values of this partition counting function’ can be obtained recursively as 
follows. Let m > 1. The (terminating) resolution of an ordered m-partition of 
I(n) may be formed from the resolution of an ordered m-partition of J(d), 
where d < n is a divisor of n, by adding the term T = dI(n/d) as separated 
terminal term to one of the partitioning sets of J(d). If d = 1, all m partitioning 
sets of I(d) are 0, so T can be added to any one of these m sets. If d > 1, T 
can be added to any one of the m partitioning sets of J(d) except that one 
which contains the terminal term of the resolution of the partition of J(d). This 
procedure uniquely delivers all ordered m-partitions of I(m), so 


Pm(n) = 1+ (m — 1) D pnd), 
the summation extending over all divisors d < » of n. Though derived for 


n > 1, this formula is also valid for » = 1, since obviously p,,(1) = 1. 


Institute for Advanced Study 
Princeton, N.J. 


*Two explicit formulae for this function are derived from the recursion relation developed 
here in Proc. Amer. Math. Soc., 3 (1952), 31-35. 











PROOF OF A CONJECTURE OF SCHOENBERG ON THE 
GENERATING FUNCTION OF A TOTALLY POSITIVE 
SEQUENCE 


ALBERT EDREI 


Let 
(1) Qo, 21, @2,... 
be a sequence of real terms with which we associate the generating power series 
(2) @+ai2+a.27+... = f(z). 
We consider the following definition due to Schoenberg [7, p. 362]: 
DEFINITION 1. The sequence (1) is said to be totally positive if the infinite matrix 


a 0 0 0 
a, ado 0 0 
G2 a, a O 
a3 G2 GQ 


(3) 


has only non-negative minors (of all orders, with any choice of rows and columns). 


Schoenberg [7] proved that the coefficients of the series (2) form a totally 
positive sequence if f(z) is a function of the form 


(4) f(e) = CHe*T] + yf IT  - 6.2), 


where A is a non-negative integer and C > 0, y > 0, a, > 0, 8, > 0, }- (a, + 8,) 
<+ o-, 

From now on, we shall assume a) ~ 0 (i.e., A = 0); obviously this is only an 
apparent restriction. 

Schoenberg [7, p. 367] also conjectured that functions of the form (4) are 
the only ones which generate totally positive sequences. 

At the International Congress of Mathematicians [Cambridge, Massachusetts, 
1950], Aissen, Schoenberg, and Whitney announced that they had proved the 
following proposition : 


A. If the coefficients of the series (2) form a totally positive sequence, then the 
series (2) is necessarily the expansion of a function of the form 


f(e) = 2°] (1 + a) / IT (1 - 8s), 
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where (2) is an integral function and a, > 0, 8, > 0, > (a, + 8) < + ©. 


In trying to prove this proposition, I was able to complete the solution of 
Schoenberg’s problem by showing: 


B. Under the assumptions of A, $(2) necessarily reduces to a polynomial of the 
form yz + 6 (y > 0). 


A detailed account of the results of Aissen, Schoenberg, Whitney, and the 
author will be given in the Journal d’Analyse Mathématique [1; 3; 8]. These 
results yield a proof of Schoenberg’s conjecture which differs substantially 
from my first approach to the subject. My excuse for presenting here this early 
arrangement is that my method essentially reduces to the construction of the 
Padé table of the generating function of a totally positive sequence. It discloses 
a number of surprisingly simple properties of such a table. Their investigation 
has been continued by R. J. Arms and the author and has led to theorems on 
the convergence of the Padé table and continued fractions associated with a 
totally positive sequence. 


1. Definition and immediate consequences of the property (P). 
DEFINITION 2. Consider the determinants 


| Om Am-1 Am—2 cee Am—n+1 | 
| 
| Im+1 Om Om-1 «++ Om—n+2) 
(n) (0) 
i> = | m2 Om+1 am ++ Am—n+3 hea @1 (m > 0," > 0) 
Omin—1 Omin—2 Omin—-3 «+--+ Om 


where a, = 0 for u < 0. We say that the power series (2) has the property (P) if 
A® > 0, for all m > 0 and all n > 0. 


The inequalities 


an >0, A®>0 (m = 0,1, 2,...), 
imply 

1 _ pe Gmt1 U1 

R = = < ay 


so that, if (2) has the property (P), its radius of convergence R, is positive.' 
It has been proved by Schoenberg [6, p. 558] that, if f(z) has the property (P), 
(1) is totally positive. On the other hand, it is obvious from the definitions that 
the converse cannot be true.” 


1By a similar argument [7, p. 362[, it is easily verified that the generating series of a totally 
positive sequence has a positive radius of convergence. 

*The connection between the property (P) and total positivity is completely clarified by an 
unpublished result of R. J. Arms. Mr. Arms has shown that, if f(z) denotes the generating 
function of a totally positive sequence, and if f(z) does not have the property (P), then f(z) 
reduces to a rational function. 








88 ALBERT EDREI 
Two observations, of purely formal character play an important part in this 
investigation. 
I. If f(s) has the property (P), then 
a 
f(— 2) 
also has the property (P). 


= by + byt + bos’ +... 


This is an immediate consequence of the elementary identities 
ay"An = bo"B,” (m > 0, > 0), 


where B® is the determinant obtained by replacing, in the expression for A, 
the a’s by the b’s. 


II. Let f(z) and g(z) denote generating functions of totally positive sequences. 
Then the sequence {c,} defined by 


f(z) g(2) = cot cet ce? +... (f(0) g(0) ¥ 0) 


is totally positive. Moreover f(z)g(z) has the property (P) if at least one of the two 
factors f and g has this property. 


It is easily verified, by direct computation, that e* has the property (P). 
Hence the function e**f(z) has the property (P) as soon as (1) is totally positive 
and « > 0 (a) ¥ 0). 


2. Poles and zeros of f(z). Let f(z) have the property (P). Consider the 
polynomials Q“”(z) defined by 


1 2 2° vial 
Omn+1 Om Gat «+. Gu-ati 
OG) = gex|Pmrt Ome Om Mmet (> 0, > 0), 
— a ee 
QO" (s) = 1 (a, = 0 for u < 0), 
S™") (2) = OO™™( — 2) (m > 0, > 0). 


We prove the following seven assertions: 


(i) the coefficients of S™™” are all positive; 


(ii) lim QO” (z) = Q@” (z) exists for all n > 0; 
(a+1) (n+1) 
(iii) lim aa An = Basi > 0 exists for all n( = 0); 
™M-sc0 m+1 m 


(iv) the polynomials Q‘@-” (z) are connected by the recurrence relations 
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Q"tP(z) = (1 — Buyiz) Q(z) (n > 0); 
(v) 2 Bs converges so that 
Q(e) = lim Q"(s) = T] (1 — Be) 


is either a polynomial or an integral function of genus zero; 


(vi) f(z) F(z) = pe + pis + pis? +... 
converges for |z| < 1/By+1, the coefficients p™ are all positive; 
(vii) F(z) Q(z) = pot pit + pw? +... 


is an integral function, the coefficients p, are non-negative. 
Proof. All the above assertions are easily deduced from the identities? 


AMY IA (n+1) 


(m,n+1) i (m-+1,m) (m,n) 44 m+1 = 
(n+2) (n+1) 
(6) grr as smn"t (2) ———s 25" (6) 40H Lan — << 2Y mn” (2). 


Assertion (i) is a consequence of a theorem of Schoenberg [6, p. 558]; it is 
also easy to deduce directly from (5). We now prove (ii). 
From (6), we obtain 


(7) Sterner (gs) = S™*t(s) — 80 eS” (8) (Yun > 0). 


—m™m 
As the cvefficients of the S’s are positive, it follows from (7) that all the 
coefficients (except the constant term, which is 1) decrease if we increase / and 
keep m fixed. This proves (ii) and also the convergence of 


-. te 


y=0 


Observing that, for a positive z, 


lim 2S” (z) > 0, 


we deduce from (5) the existence of the limit in (iii). This proves (iii) and (iv) 
simultaneously. 
Consider (7) for 1 = ©, m = 0, and compare the coefficients of 2; this yields 


Bi + Bot... + Bai < @1/do, 
which proves (v). 
To prove (vi), we consider the expansion 





*These formulae are straightforward consequences of Jacobi's theorem on the minors of the 
adjoint of a given determinant. 
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(8) f(2)Q™” (z) = a pe” 2" 
and observe that the well-known identity [5, p. 430] 
| dy Gy—1 Ay—2 ee Ayn | 
i< Gm+i1 Om Am-1 eee On—o+1| 
(9) fo ”(z) = A™ ,* | Om+2 Om+1 am «+ + Om—n+2 | ’ 
nn 
Rs Omin—1 Omtn—2 «+--+ Om 


gives us explicitly the coefficients p,°"”. Now the sequence of functions 


(10) f(z) \o"™"(sz (z) (m=0,..., @) 
is regular in the circle |z| < R, where R denotes the radius of convergence of (2). 
For m fixed and m — , (10) converges uniformly to 
(11) f(z) OP” (z) = po + piM2? +..., 
in any circle |z| < R* < R. 
From (7) we deduce 


(12) f(z) Ow") = fe) Oz) + :> Yeni f(z) O**-?(2)}. 


=m 


Using (8), (11), (12) and Weierstrass’s double-series theorem, we obtain 


(13) pe = py” + 2 ta per” Sete e.. 3 
The total positivity of (1) implies p&“” > 0, when k < m, so that, (13) yields 
(14) > pp” (k < m). 
Observing that 

Pmii = 0, 


and using (14), we also obtain 


(15) oS. < a? + Ym+1i.n—1 + Ym+2.n—1 + ** _ (m > 0, n > 1). 
Putting 


im sup|(p,, ) = 0, coi } & Se: 
(16) li \(pe’)**| ( 2 


mM +c 


we deduce, from (15), 


(17) oe £ ex 1 lim sup (x . .)™ 


mid 


Moreover, as all series of the form 


D Yun 


“=m 


are convergent, (17) trivially implies 
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(18) On <S On-1 (mn = 1, 2, 3, ) 
We now verify the inequalities 
(19) Cx < Best (k => 0, 1, _ ee a, 


by an induction over k. 
We note that p® = a,, implies o) = 8,, and assume that (19) is true for 
k=0,1,2,...,2”— 1. We distinguish two cases 


(a) B62... 8, = 0, 
(b) BiBe eee B, ~ 0. 


The case (a) is immediately settled by observing that our induction assumption, 
and (18), then yield o, = 0, for k > n. 
In the case (b), we observe that 


es) l/m | 
(20) lim sup ($ ra0-1) | < Bass, 
is trivial if Basi > Ba. 


If Bai < 8,, we return to the formulae (6) which define the quantities ma 
as ratios of determinants. 
Observing that (iii) implies* 


(n+1) 
lim Aart = 6:B:... BeBor, 
(21) lim (A)'"| = Bib2... Be 


we obtain 
| = tn Aett_4e) 
im raas)™ = tim (44a) = Be 
We may then compare 
} Yun-1 
=m 
to the remainder of a geometric series. This comparison leads again to (20). 
Combining (17) and (20), we complete the induction which proves (19). As 
(13) implies p” > 0, assertion (vi) is also proved. Assertion (vii) follows from 
(vi) and Weierstrass’s double-series theorem. Although the fact is not used in 
this paper, we note that (19) can be replaced by 


(22) On = Bers (k = (), > er * 





‘It should be noted that the validity of (21) is independent of the arguments following (14). 
Interrupting our proof after (14), and using (21), a direct application of Hadamard’s theorem 
on polar singularities (see, for instance [2, p. 333])would have enabled us to verify that f(z) 
is a meromorphic function. However, as the wording of Hadamard’s theorem does not meet 
all our requirements, and as we do not need its more delicate parts, I have replaced its use by 
simple and elegant arguments which I owe to R. J. Arms. 
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To verify this, we assume that, for a definite value of k( = nm), a, < Basi. The 
point 1/8,,: would then be interior to the circle of convergence of the 
expansion of f(z)Q@."(z). The function f(z)Q@"*»(z) would be regular and 
would vanish at z = 1/8,4:. This is impossible in view of the positivity of the 
coefficients p**. We finally observe that (18) and (22) imply 


Aii>be>8s> 


3. Proof of Proposition A. The function e*f(z) certainly possesses the pro- 
perty (P); it is therefore sufficient to prove the proposition for functions which 
have this property. 

By I and the seven assertions, there exist two integral functions P(z) and 


Q(z) such that 


P(z) = Tl (i — as) (a, >0,>.a,< + @), 


v=l 


Q@)=-Ta-6s) >0D8<+=) 


and such that both P(z)/f( — z) and Q(z)f(z) are integral functions. Moreover, 
when z > 0, the function P(z)/f( — z) does not vanish because, by assertion 
(vii) the coefficients of its expansion are real, non-negative, and ap~' = P(0)/f(0) 
~ 0. Hence f(z)/P( — 2) is regular for negative values of z and 
pote) = FOSe) ~ ois) 
is obviously an integral function. Similarly, 1/v(z) is an integral function so that 
v(z) has no zeros. 
Hence 
v(z) = &™, 


which proves the theorem of Aissen, Schoenberg, and Whitney. 


4. Proof of Schoenberg’s conjecture. As (1) is totally positive, the coeffic- 
ients of the function 


f*(z) = do +awtag2?+...+ayn2+.... 


also form a totally positive sequence. Applying Proposition A to both f(z) and 
f*(z), we can then assert the existence of 


1. two integral functions ¢(z) and ¢*(z), 


2. four convergent series of non-negative terms, say 


- mm at we. 
such that 
(23) jie) = vol] G+ a) 





I] @ — 6,2)’ 


Eee 


ei 





p = FF CF 


at 


eee 
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e *ol] fit <2). 
(24) f(s) = pg 'v 
Now 
f(z) thr). = f*(z*) 


using (23) and (24), this becomes 
(25) sol] + a,Z) + eo oft — a,z) = 2 o * 1] (1 + a* ry ‘) 


Tl @ — 32) Tla+é2)  [la-é<.) 
Put o(2) = yo + yiz + yor? +... and 


e+rne+... Ty (L+a4,2)(1 + 6,2) 
%6 3) = ett 1 + Ys2* + ) (a+. ee. 
( ) u( ) vel (1 inal a,z)(1 _ B,2) 
equation (25) is then equivalent to 


on g(a) -o(-) FT (1 + a*,2" )(1 + 8,2) 
(27) writes LL (r— 6%") (1 — at)” 


This relation, together with (26), gives us complete information on the 
distribution of the values 0, ©, and — 1 of the meromorphic function u(z). 

The value — 1 has an exponent of convergence not exceeding 2, whereas both 
the zeros and the poles have an exponent of convergence not exceeding 1. 

By a well-known theorem of Nevanlinna [4, p. 72], the order of u(z) does not 
exceed 2, so that 

O= 72 = ¥s = ¥7=...-- 

Hence f(z) is of the form 


ad Vittyat*+yett+ . . JTa+ a,zZ) 
(28) f(z) = Ce Tla — 2s) (C> 0). 


Now if ¢ > 0, both e** f(z) and e**/f( — z) have the property (P), so that, by 
assertion (vii), the power series expansions of 


u;(z) —_ Cele trvetres trast cee I] (1 4. a,z), 


v=l 


and 


ua(2) = Cotetrrererer—-- TT (1 + 8.2), 


ven! 
have all their coefficients non-negative. 
Hence, putting 


M,(r) = max \e1(z)|, M,(r) = max \ue(z)|, 
we have 


Mi(r) = u(r), M,(r) = u2(r), 
and 
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(29) Myi(r)M2(r) = eS" Il (1 + a,r)(1 + 8,r). 
v= 
As M:(r) > M:(0) = C-, 


we see that the increase of M,(r) is dominated by the increase of an integral 
function of order 1. Hence 


O=72= 11 = ¥e=..-.-.- 


We complete the proof of Schoenberg’s conjecture by proving y = y: > 0. 
If y were negative, we could choose « > 0 so that 


(30) y¥+e< 0. 


Now the increase of 
I] (1 + ar) (1 + 6,7), 
is at most of the minimal type of order 1. Hence (29) and (30) would imply 


lim M,(r)M;(r) = 0, 


which is contrary to the maximum-modulus theorem. 
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THE BOOLEAN ALGEBRA OF REGULAR OPEN SETS 
R. S. PIERCE 


1. Introduction. Let S be a completely regular topological space. Let 
C(S) denote the set of bounded, real-valued, continuous functions on S. It is 
well known that C(S) forms a distributive lattice under the ordinary pointwise 
joins and meets. For any distributive lattice L and any ideal JC L, a quasi- 
ordering of L can be defined as follows: f > g if, for allh € L, f (\h € I implies 
g(\h € I. If equivalent elements under this quasi-ordering are identified, a 
homomorphic image of L is obtained. In this paper, the particular case where L 
is C(S) and J is any principal ideal will be studied. It will be shown that the 
homomorphic image obtained in this way is isomorphic to a sub-lattice of B(S), 
the Boolean algebra’ of regular open sets of S. This homomorphic image will be 
denoted by %(5S); it does not depend on the choice of the principal ideal J. For 
the case where S is a normal space, a topological characterization of (5S) will be 
obtained. Finally it will be proved that 8(S) is isomorphic to the normal com- 
pletion (see [1}) of 2(S). Incidentally these results prove, without using trans- 
finite methods, that for any completely regular space S, 8(S) is determined to 
within isomorphism by C(S), a fact which could also be inferred from the 
theorem of Kaplansky [3] that a compact Hausdorff space S is determined to 
within homeomorphism by C(S). 


2. A congruence relation. In everything that follows, S will denote a 
completely regular topological space. Denote by <, (\, the ordinary point- 
wise inequality, meet and join operations in C(,S). 

Definition 2.1. f D gif, for allh € C(S),f Oh < 0 implies? gh < 0. Write 
f~egiff2gandg of. 


LEMMA 2.1. On C(S), > is a quasi-ordering and ~ is a congruence relation. 


Proof. lf f>g and g Dh, then f > h is an immediate consequence of the 
definition. Also f > f is clear. Thus > is a quasi-ordering and consequently ~ is 
an equivalence relation (Birkhoff [{1)}). 





Received September 17, 1951. The author wishes to thank Professor R. P. Dilworth for 
his generous help and inspiration. 

1An open subset of a topological space is called regular open if it coincides with the interior 
of its closure. It is well known (see Birkhoff [1, p. 177]) that the regular open sets form a 
complete Boolean algebra. 

*This definition puts the zero function in a distinguished position—a position which the 
lattice structure of C(.S) does not support. It is easy to see however that if 0 is replaced by 
an arbitrary (fixed) fo€ C(5S), all the results of the paper go through with only slight modi- 
fication. 
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Let f >g and f; be any continuous function; then f/\f;: > g\f; and 
fUfi2eVUfi For suppose (fOf:1) 14 <0. Then fA (fi Oh) <0, so 
since f > g, g(\ (fil \h) < 0. Thus (g fi) 4 < 0 and, since & was arbi- 
trarily chosen, f (\f; > g(\fi. Next suppose (fU f:) (1\ h < 0. This means 
fi-Q\h <0, fO.\hk < 0 and, because f D> g, g(\h < 0. Hence (gU fi) Nh = 
(g\h) U (fi Oh) < 0 and again by definition f Uf; D> g U fi. 

It follows that if f~g, then fM\fi~gO\fi and fUfi~g Uf; for any 
fi € C(S). Thus ~ is a congruence relation and the proof of the lemma is 
complete. 


LemMa 2.2. f Dg if and only if * {x| f(x) > 0}- D {x| g(x) > O}-. 


Proof. First suppose {x| f(x) > 0}- D {x g(x) > O}-. Let h € C(S) be 
such that ff\h < 0. It will be shown that this implies g/\h < 0, so (by 
definition) f > g. 

If xo € {x| g(x) > O}-, then by hypothesis, x» € {x | f(x) > 0}-. This means 
that if N is a neighbourhood of xo, there is an x € N with f(x) > 0. Since 
fC\h < 0, it follows that h(x) < 0. Thus x» € {x | h(x) < O}- = {x | h(x) < 0} 
and h(xo) < 0. Consequently (gM h) (xo) < h(xo) < 0. If xo ¢ {x| g(x) > 0}-, 
then g(xo) < 0 and (g/\h)(x0) < g(xo) < 0. Hence for all points x» of S, 
(g (\ h) (xo) < 0; that is, gh < 0. 

Now suppose f D g. Let M = {x| f(x) > 0} ~*\ {x| g(x) > 0}. Clearly, M 
is open. Also if y € M, then f(y) < 0 and g(y) > 0. To show that {x | f(x) > 0}- 
> {x| g(x) > 0}-, it suffices to prove that M is empty, since then {x | f(x) > 0}- 
> {x| g(x) > 0}, and taking closures on both sides gives the desired result. 

Suppose M is not empty. Let x» © M. By complete regularity, 4 € C(S) 
exists so that A(x9) = 1 and A(y) = 0 for all y¢ M. Then (g/\h)(xo) = min 
{1, g(xo)} > 0, so g/\h non < 0. On the other hand, (f (\h)(y) < f(y) < 0 if 
y € M, while (f (1 h)(y) < Aly) = Oif y ¢ M. Thus f 1h < 0. This contradicts 
f > g and shows that M is empty. 


Denote by &(5S) the set of equivalence classes under ~. Then by Lemma 2.1, 
2(S) is partially ordered under > and, with this ordering, it is a distributive 
lattice. If f € C(S), let A, denote the congruence class containing f. 


LEMMA 2.3. &(S) is isomorphic to the sub-lattice of the Boolean algebra of reg- 
ular open sets of S which consists of all sets of the form I({x | f(x) > 0}-) where 
f € C(S). 


Proof. The mapping A, — I({x| f(x) > 0}-) is one-to-one by Lemma 2.2. 
Also* 


*For any set NC S, N~ will denote the closure of N, N° the complement of N in S, while 
the interior of N, N°~°, will be represented by I( NV). 
‘The lattice operations in the Boolean algebra of regular open sets are given by PAQ = 


I(UPN Or) = PO\O, PVO = KP-U0-) and P’ = P~. 
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A sng > T( fx | (F Og) (x) > O}-) = I([{x| f(x) > 0} A {x| g(x) > 0} 
I({x | f(x) > 0}-) A I({x| g(x) > 0}>, 
I({x | f(x) > 0}- U {x | g(x) > 0}-) 
I({x | f@) > 0}-) Vv I({x| g(x) > 0}>>). 


Ay, > I({x| (f U g)(x) > 0}-) 


The following lemma gives information on the image of 2%(S) which will be 


needed later. 


LEMMA 2.4. Let P bea regular open set of S. 
(a) If P # 0 (the null set), f © C(S) exists so that 0 # I({x | f(x) > 0}-)CP. 
(b) If P # S (the unit set), f € C(S) exists so that S # I({x | f(x) > 0}-)DP. 


Proof. (a) Let x» € P. By complete regularity, choose f € C(S) so that 
f(xe) = 1 and f(y) = 0 if y¢P. Then xo € {x| f(x) > 0} CP so x € I({x| 
f(x) > 0}-) SC I(P-) = P. 

(b) P- # S since otherwise P = J(P-) = 1(S) = S. Let x» € P~*. Choose 
f € C(S) so that f(xe) = — 1 and f(y) = 1 for all y € P-. Then 


PC {x| f(x) > 0} C I({x| f(x) > 0}-) C {x| f(x) > 0} ¥S. 


It may happen that the image of &(S) under the mapping of Lemma 2.3 
coincides with the whole Boolean algebra of regular open sets. For this to be the 
case, it is clearly necessary and sufficient that every regular open set P of S be 
of the form P = I({x | f(x) > 0}-) for some f € C(S). This condition prevails 
in two cases of special importance. The first is that in which S is a metric space. 
For then, if P is a regular open set and f € C(S) is defined by f(x) = p(x, P*), 
p being the metric of S, P = {x| f(x) > 0} = I({x| f(x) > 0}-) and the 
criterion is fulfilled. The second important case is that of an extremal space, that 
is, a space in which the closure of every open set is open. In such a space, the 
regular open sets are evidently just the open and closed sets. Hence the charac- 
teristic function ¢p of a regular open set P is continuous and P = {x} dp(x) > 0} 
= I({x | p(x) > 0}-). An important consequence is: 


THEOREM 2.1. The Boolean algebra of regular open sets of a completely regular 
topological space S is isomorphic to the complete Boolean algebra obtained by re- 
ducing the normal completion of C(S) modulo the congruence of Definition 2.1. 


Proof. By the theorem of Dilworth [2], the normal completion of C(S) is 
isomorphic to C(S) where © is the Boolean space associated with G(S), the 
Boolean algebra of regular open sets of S. Then (see Stone [4; 5]) S is extremal 
and its Boolean algebra of open and closed sets is isomorphic to @(S). This 
means that S and S have the same Boolean algebra of regular open sets. The 
theorem then follows from the above remark. 


In spite of the above results, it is not true that, for all completely regular 
spaces S, the image of &(S) is the whole Boolean algebra of regular open sets. 
This is shown by the following example. 
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Let S be the set W, consisting of all ordinal numbers < Q (the least ordinal 
of the third class), together with the real interval J = [0, 1], where Q is identified 
with 0. Simply order S by defining a < x for all a € W and x € J, and by re- 
taining the usual order in W and J. Impose the interval topology (see [1]) on S. 
It is easily seen that S is completely regular. Now the interval P = {a € S|a < 
Q} is a regular open set of S. Suppose f € C(S) exists with P = I({a| f(a) >0}>). 
Then f(a) < 0 holds for all a > 0 = Q. Let a, = sup {a f(a) > 1/n}. For all 
N, & <Qsince {ae S | f(a) <1/n} is open and contains J. Thus a, is an ordinal 
of the second class. a = sup a,, as the limit of a sequence of ordinals of the 
second class, is itself an ordinal of the second class. Hence a < 2. Because 


{a | f(a) > 0} = U {a | f(a) > 1/n} C U {a|a <a} C fa |a <a}, 


it follows that J({a | f(a) > 0}-) C {a| a < a} # P. This is a contradiction, 
showing that no such f € C(S) can exist. 

The preceding example indicates the need for a topological characterization 
of £(S). For the important case where S is a normal space, the following theorem 
answers this need. 


THEOREM 2.2. Let S bea normal topological space. In the mapping M — I(M—) 
of the lattice of open sets onto the Boolean algebra 8(S), L(S) is the image of the 
sub-lattice of all open sets which are the countable unions of closed sets (that is, 
open F, sets). 


Proof. %(S) consists of all sets of the form P = I({x| f(x) > 0}-) where 
f € C(S). But ; 
{x | f(x) > 0} = U {x | f(x) > 1/n} 


is an open F, set. Thus every P € &(S) is the image of an open F, set. 
On the other hand, suppose P = J/(M~—) where M is open and 


M= UF, 
n=l 


is the countable union of the closed sets F,. By Urysohn’s lemma, for each 2, it 
is possible to select f, € C(S) so that f,(x) = 1 for all x € F,, f,(x) = 0 for all 
x € M‘and 0 < f, < 1. Define 


f = > (2" fn: 


Then f € C(S) and f(x) = 0 if x € M*. If x € M, then x € F, for some n, 
whence 


f(x) > 2" fax) = 2" > 0. 


Consequently M = {x | f(x) > 0} and P = I(M-) = I ({x| f(x) >O0}-) € 2S) 
Thus the image of any open F, is in &(S) and the proof is complete. 
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3. Completion. Having constructed the homomorphic image of C(S), it will 
now be shown that the normal completion of £(S) is isomorphic to the Boolean 
algebra of regular open sets. 

Let it be recalled that the normal completion of a lattice & consists of the 
normal subsets (or closed ideals) of & ordered by inclusion, that is, those subsets 
which contain all lower bounds to the set of their upper bounds. The mapping 
A—{Bé 2|BCA} imbeds & as a sub-lattice of its normal completion (for 
proofs, see Birkhoff [1)}). 


LEMMA 3.1.° Let & be a sub-lattice of the complete Boolean algebra B. Suppose, 
moreover, that the following dual conditions are satisfied: If X € Band X #Z 
(zero of B), then A € & exists with Z# ASX. If X € Band X # I (unit of 
B), then A € & exists with I # A D X. Then the normal completion of & is 
isomorphic to B. 


Proof. For X € , let Ny = {A € &| A CX}. The first task is to show 
that Ny is a normal subset of &. 

First note that, if A > B holds forsome A € Qandall B € Ny, then A D X. 
For otherwise A’ (\ X # Z and, by hypothesis, B € & exists with A’ (\ X D 
B # Z. This means A non 2 B while BC X which contradicts the assumption 
on A. Thus the upper bounds in % of all B € Nx are precisely those A satisfying 
A D> X. As a result, the preceding argument can be dualized to show that a 
lower bound, A, of the set of all upper bounds of Nx satisfies A C X. Consequent- 
ly Nx is normal. 

Next, every normal subset of & is of the form Ny for some X € B. In fact, if 
M is normal, let X = UM (sup in B). Then if ACM, ACX, soMC 
Nx. On the other hand, by the definition of X, if C D> B holds for all B € M, 
CD> X. Thus A € Ny (that is, A C X) implies A C C for all upper bounds 
of M. Since Mt is normal, this means A € Mt, so Ny € Me. 

The above argument shows X — yx is a mapping from $8 onto the normal 
completion of %. The mapping obviously preserves order: X, C X, implies 
Nx, S Rx,. To complete the proof, it is sufficient to show that this mapping 
is one-to-one. 

Suppose X, non C X». Then Z # X,f)\ X2', so A € L exists with Z # A C 
X,0\ X2/. It follows that A € Nyx, and A ¢ Ny,, that is, Ny, non C Ny,. 


CorROLLarRyY. If S; and S: are two completely regular topological spaces such that 
£(S,) is isomorphic to &(S2), then S, and S, have isomorphic Boolean algebras 
of regular open sets. 


According to Lemma 2.4 the above principle is applicable to the image of 
£(S) in the complete Boolean algebra of regular open sets. The result of this 
application is: 


‘This lemma, in essence, was presented by Professor Dilworth during his 1950-51 seminar 


"git hy 
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THEOREM 3.1. The Boolean algebra of regular open sets of a completely regular 
topological space S is isomorphic to the normal completion of the lattice obtained by 
reducing C(S) modulo the congruence of Definition 2.1. 
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CONTINUOUS SOLUTIONS OF THE FUNCTIONAL 
EQUATION /*(x) = f(x) 


G. M. EWING AND W. R. UTZ 


In this note the authors find all continuous real functions defined on the 
real axis and such that for an integer n> 2, and for each x, 


(1) f"(x) = f(x). 


The symbol /”" denotes f iterated n times. 
The following two classes of functions occur as solutions. 


Cxass I. 
(a:) The function f(x) 1s continuous for all real x, 
(bi) f(x) = x on a connected subset S of the x-axis, and 
(cx) g < f(x) < G, in which g and G denote, respectively, the infinum and 
supremum of f(x) on S. 


The set S must be a point, a closed interval, a closed ray, or the entire 
x-axis. Thus Class | includes all constants and the function x. If S is a closed 
interval (a, b] then f(x) is arbitrary outside of [a, 5] except for continuity and 
the condition a < f(x) < 6. If S is a ray, f(x) is similarly described. 


Cxass II. 
(a2) The function f(x) is continuous for all real x and either, 
(be) f(x) = x, or 
(co) f?(x) = x on a non-degenerate closed interval |a, b], f(a) = 6, f(b) = a, 
anda < f(x) < b. 


A function satisfies (bz) if and only if y = f(x) implies x = f(y). Its graph is, 
accordingly, symmetric with respect to the line y = x. If f(x) is a solution of 
(be), then the inverse of the transformation x -+ f(x) is clearly single valued 
and continuous. Hence the transformation x — f(x) defines a homeomorphism 
of the x-axis onto itself. 

One can easily see that Classes | and II have only the function x in common. 


LEMMA 1. f(x) is a continuous solution of f?(x) = f(x) if, and only if, f(x) 
is of Class 1. 


Proof. That every Class I function is a solution is easily verified. Con- 
versely, if f(x) is a continuous solution then x = f(r) satisfies f(x) = x for 
every real r. In case f(x) = x has only one solution, f(x) is constant and hence 
of Class I. If f(x) = x has two solutions a and 6, a < b, then f(a) = a and 


f(b) = b; and given c between a and 6 there exists, from the continuity of f(x), 
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a number m, a < m < b, such that f(m) = c. It follows that f(c) = ff(m) = 
f(m) = c, and hence if f(x) = x at the ends of an interval the relation holds 
identically on the interval. The maximal set S on which f(x) = x is thus con- 
nected. The continuity of f(x) implies that S is closed. Finally f(x) has property 
(c:), for if there exists an 7 not in S such that f(r) does not satisfy (ci), the 
relation f(x) = ff(x) = f*(x) with x = r, contradicts the fact that S is maximal. 


LemMMA 2. If f"(x) = x on a non-degenerate closed connected subset S of 
the real axis and if f(x) maps S continuously into S, then 

(i) f(x) is a homeomorphism of S onto itself, 

(ii) tf S is an interval [a, b), f(x) = x on S or f?(x) =x on S and f(x) is 
equivalent to a reflection of |a, b| about the single fixed point p, 

(iii) af S is a ray, f(x) =x on S, and 

(iv) if S is the entire axis, f(x) = x on S or f(x) = x on S and f(x) ts equiva- 
lent to a reflection of S about the single fixed point p. 


Proof. Conclusions (i) for the case of an interval and (ii) are special cases 
of results in Whyburn [2, pp. 240, 264). 

If Sis a ray, the mapping [2, p. 240] is (1-1) and onto. Thus f(x) is monotone 
on the ray and the end point is fixed under f(x). If there were an interior 
point 5 of the ray such that f(b) + 5, the monotonicity of f(x) would imply 
that f"(b) + b. Hence (iii), which implies (i) for the ray. 

If S is the real axis, the mapping is again (1-1) and onto and f(x) is mono- 
tone. If f(x) increases with x, we see that f(x) = x by the argument employed 
for the ray. If f(x) is monotone decreasing its graph cuts y = x in exactly one 
point and f(x) is topologically equivalent to a reflection of the x-axis about 
the abscissa of this point. 


COROLLARY 1. Jf mn is odd, the functional equation f"(x) = x has only the 
function x as a continuous solution. If n is even, the continuous solutions of 


f"(x) = x are those of f?(x) = x. 


Proof. By conclusion (iv) of Lemma 2, there are two possibilities. If 
n= 2m+1 and if f(x) =x then P(x) =f... f(x) =x, and hence 


f(x) = f(x) = x. If m is even, the stated result is immediate from Lemma 2 
since f(x) = x is a solution of f?(x) = x. 


THEOREM 1. The continuous real solutions of f"(x) = f(x), nm > 2, are the 
functions of Class | if n is even and the functions of Classes | and II if n ts odd. 


Proof. lf f(x) is of Class I then f?(x) = f(x). Whence 


f(x) = f(x) = f(x),.... 
If f(x) is of Class II we verify that f*(x) = f(x). Then 
f*(x) = f(x) = f(x), f(x) = f(x) = f(x), .... 


Conversely, let f(x) be a continuous solution of f"(x) = f(x). Then 


ft fle) = f-2 PM(x) = f"-2 f(x) = f(x) 





A FUNCTIONAL EQUATION 103 


so that f"~'(x) if of Class I by Lemma 1. Let S be the maximal subset of the 
x-axis on which f*~(x) = x. 

If S is a point, then f*""(x) = c and f(x) = ff"""(x) = f(c), so that f(x) is 
of Class I. 

If S is the closed interval [a, 6] then f*~"(x) = x on [a, b] anda < f*"(x) <b 
by (bi) and (c,). Moreover a < f(x) < 6, as a consequence of the relations 


(2) f(x) = f(x) =f" f(x). 


Thus f(x) maps the real axis into [a, 5]. In particular [a, 6] goes into [a, 6] and 
Lemma 2 is applicable. If f(x) =x on [a, 6], f(x) is of Class 1. The other 
possibility is that f*(x) =x on [a, d], and that f(a) = b, f(b) = a, in which 
event f(x) is of Class II. 

If f"""(x) = x on aray, we see that f(x) maps the reals into [a, ©] or [— ~, 5]. 
Hence f(x) = x on the ray by Lemma 2 and is of Class I. 

Finally, if S is the x-axis, the desired conclusion is given by Corollary 1. 


The functional equation 


(3) f"(x) = f(x), 

m and n integers, 1 < m < n, has among its continuous solutions the functions 
of Classes I and II if m + nm is even and those of Class I if m + n is odd. 
In either case, if f(x) is a solution of (3) there exists an integer k such that 
f*(x) is of Class I. However, each equation (3) has continuous solutions in 
neither of our classes. For example f*(x) = f?(x) has the solution 


sin x x>I1 

— 0 x| <1 
(=) = —x-—l —-2<qx< -l 
l x <-—2. 


Let y = f(x) denote a transformation of a topological space X into itself. 
A necessary and sufficient condition that f*(x) = f(x) have f(x) as a continuous 
solution is that the set S of fixed points under f(x) be non-vacuous and that 
f(x) map X continuously into S. That f(x) be a continuous solution of f*(x) 
= f(x) isa restriction on both S and f(x). If X is the real axis these restrictions 
are given by our Lemma 1. In general the solutions of f*(x) = f(x) are the 
“retractions’” of the space X. These have been studied by Borsuk [1] but 
results for higher-dimensional or abstract cases comparable to those of 
Lemma 1 do not seem to be available and appear difficult to achieve. 
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GENERATING FUNCTIONS FOR BESSEL AND 
RELATED POLYNOMIALS 


E. D. RAINVILLE 


1. Introduction. Krall and Frink [4] aroused interest in what they term 
Bessel polynomials. They studied in some detail what may, in hypergeometric 
form, be written as 


(1) Ya(x, a, 6) = oFo( — n,a — 1+; — ; — x/d). 


The simple Bessel polynomial is the special case a = b = 2. In essence the 
polynomials considered by Krall and Frink, and by others [1; 3], are the 
terminating 2F>o’s with numerator parameters ( — m) and (c +m), m a non- 
negative integer. We shall therefore work with 


(2) oFy( — 2, ¢ +"; — ; — x) 
for brevity. From (2), the polynomial (1) is easily obtained. 

Burchnall [1] produced a convergent generating function for y,(x, a, 0). 
It is the purpose of this note to obtain relations which yield in simple special 


instances two generating functions, one of them Burchnall’s, for Bessel 
polynomials. 


2. Generating functions of Burchnall type. Consider the polynomials 
(3) ¥a(x) 


—2,c+n,1— Bi —2,1— B—2,...,1— Be — 2; 
= oF, (— 1)?** "x 


1—a,—n,1—az:—n,...,1—a,—n; 
in which the c and all the a's and §’s are independent of nm, and in which m is a 
non-negative integer. The a’s and #’s are not to be non-positive integers, but 
are otherwise unrestricted. 
If in (3) no a’s and no #’s are used, the resulting polynomial is that of (2), 
the generalized Bessel polynomial. 
We now prove that 


2 e-1 Oi, 2, ~. ++ » Ap; 2 
a oe — a —— 
9 M25+++s Mes 


_ > Walx) (@1)n(@2)n - - - (Opn 
- 2, n! (B1)n(B2)n eee (By)n . 


in which (1 — 4xt)!-+1 as t-—>0. Burchnall’s [1] generating function is a 
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special case of (4) in which there are no a's and no §’s, x is to be replaced by 
(x/b) and ¢ by (4dt). In Burchnall’s paper the misprint, {} — 4(1 — 2xt)*}2-« 
where {} + 4$(1 — 2xt)!}?* was intended, should be corrected. 
For our proof of (4) we need a simple result from the theory of the »F;. 
LEMMA. 
| 2 P. 
F Be De = on BE celine 
oFi(y,7 + 4; 27; 42) (1 42) 1+ v (1 — 42) ° 
To obtain the result in the lemma, use Gauss’s formula 
oF; (a, b; 2b; 4x(1 + x)~*) = (1 + x)™ oFi(a,a + § — 536 + 3; x’), 
for which see Magnus and Oberhettinger [5], put b = y, s = y + 4, x(1 + x)" 
= z, and use the fact that a ,F> is a binomial. 
Since (1 — a — n), = ( — 1)*(a@),/(a),-x, it follows that 
Wn (x) (a1) n (2) n sss (a) - . (C) m+n (01) n—e (2 ne “+ . (ay \n—2 x* 
n! (B1)n(Be)n “+. (Bo)n k=0 (C)n(B1)n—e(B2)n—z **. (Ba)n a(n — k)! ki 


Therefore, 


= n(x) (a1) (2) n ese (cry) nt” 
2, n! (B1)n(B2)n eee (By)n 


cs) 











: (©) nsx (1) n—z (2 ne ** « (ay) nx xf 





nad bmd (C)n(B1)n—e(B2)n—e - -- (By)n—e( — k)! k! 





a * (C)ns-20(c%1)n(ai2)n - - » (tp) xe'0"** 
n=0 k=O (C)n+e(B1)n(B2)n see (By)n n! k! 
= FF C+ M)nlxt)’ (ar)n(or)n -- - (p)nl” 





mao eo (C+) R! © (Bi)n (Bon. -- (By)n 1! 


D aFi(de + 4m, bo + dm + Bsc + m; det) Sr dalaadn --- (@p)af 


n=0 (B1)n(Be)n eee (Bo)n n! 
= = ——— - | enaten + (@%y)n f 
a 2d, (i txt) F + V (1 _ 4xt) (B1)n(B2)n “ee (Be)n n! 


™ y 2, «+ + y My, 
(1 — 4xt)™ a Oak a F, ‘ 2 ; 
‘ 1+ (1 — 4xt) "Bs. Bs g.. 1 + V{l — 4xt) 
ly p++ Bes 


which completes the derivation of equation (4). The method used is nothing like 
that in Burchnall’s work. 

If p < q + 1, as is true for the generalized Bessel polynomials, the generating 
function series, the ,F,, even has a region of convergence near ¢ = 0. For gene- 
rating function purposes.the luxury of convergence is, fortunately, not at all 
necessary. 


3. Another generating function for Bessel polynomials. Consider next 
the polynomials 
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= — M,C PM, Oy, Oa, 2 2 Uys | 
(S) a(x) raPf ee tae 


of which the Bessel polynomials are again a special instance, namely that in 
which no a’s and no #’s are used. This extension of the generalized Bessel poly- 
nomials is, of course, distinct from that in equation (3). 

We shall obtain the generating function relation 


$c, 3c + 4, cs, ae, . . . , Op; dxt = g(x) (c)qf" 

atn|- Ee“ 

Bis Bares Ba; 

which contains a divergent generating function, a Fo, for the Bessel poly- 

nomials, namely the choice of no a’s and no §’s. The series in (6) converge near 
t=O0ifp<q-1. 


(6) (1 — t) “42s 








We find that 
oale)(clat” _ Sd yh ( — m)alc + m)s(c)a(ar)s(as)s--- (@p)e( — 2)'F 
2d, n! Dy 2 (B:)e(B2)e.. ~ (By) R! mn! 
a > “  (C)npe(ar)e(2)e - « « (tp) xt" 
n=0 t= (B1)e(B2)e. . . (By) R! (mn — k)! 
- Ss (C)npze (1) x (a2) - - - (ay), x't"™* 





n=0 k=O (B1)e(Bo)e.. - (By)y R! n! 


mF Ft CA 2k )n (Chan (1) 2 (as) e - - - (pe (xt) 
2, 2, m! (B:)z(B2)e.- - (Bade R! 
= Fi — anlar )e(ar2)s - - - (ap) (xt) 
£0 (B1)e(Bo)e... (Bye RE (A — 8)" 
a bc, 3c + 4,01, a2, ... Gy; At 
= (1 — ty p+2Fe @—o:|' 
Bas Bays cios Ba; 


which completes the derivation of equation (6). 
For c = 1, the equivalent of equation (6) has already appeared in Sister 
Celine’s work [2]. 
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ON TWO-DIMENSIONAL BERNSTEIN POLYNOMIALS 
P. L. BUTZER 


1. Introduction. Let the function of two real variables f(x, y) be given over 
the square 


S:0<¢x<1, 0< y<1. 


Then the Bernstein polynomial in two variables x and y, corresponding to the 
function f(x, y) is defined as 


(1) Bin 9) =D DA, *)p, ws (#) Prans(9). 


vy, =) v,=0 


where 
Prn(u) = ("\wrc —u)"”. 


D Prn(u) = 1. 


v==0 


Obviously 


The main purpose of this paper will be to prove the following result. 
THEOREM. If the function f(x, y) is bounded in the closed square S, then 
" . a” 


a” , 
es on “ap 

(2) = ax4 dye * Bin, (Xo, Yo) ax} ays * f (xo, Yo) 

at every point (Xo, yo) belonging to the open square 0 < x < 1,0 < y < 1 where 

the pth total differential of f(x, y) exists, provided ny, nz tend to infinity in such a 

way that 


(3) 0<reBti<s<te. 

If, moreover, the partial derivatives of f(x, y) of the first » orders exist and 
are continuous in S, then the relation (2) holds uniformly in S as m1, m2 approach 
infinity in any manner whatsoever. This latter result, which is used in estab- 
lishing our theorem, was recently proven by Kingsley [1], but we shall deduce 
it in a more direct manner. 


2. Preliminary results. If 5 > 2,;~*, 0 < a < 4 then for every k > 0 there 
is a constant C (depending only on a and k) such that 


Received March 5, 1952. This paper forms a portion of the author's thesis, On Bernstein 
Polynomials, written under the direction of Professor G. Lorentz and accepted for a Ph.D. 
degree at the University of Toronto in November, 1951. The author gratefully acknowledges 
the assistance received from the National Research Council, Ottawa. 
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(4) > Prn(u) < or, 


|¥ wis 


Furthermore, for a given y > 0 and any 0 < u <1, 


(S) u ly — nul” pr.n(u) < Cn” 


= 


where C; is a constant depending only on y. 
For the proofs of these inequalities see [2]. 


Ocul. 


We now recall several facts on differentials. By definition, the pth total 


differential of f(x, y) exists [3] at the point (xo, yo) if 
(6) f(x, y) = P(x, y) + g(x, y) 


where P(x, y) represents a polynomial of degree p, and 


? 


(7) g(x,y) = > a(x, y)(x — x0)” ‘(y — yo) 


i=0 


for the sake of brevity, where the a;(x, y) have zero as a unique double limit as 


x—> Xo, ¥ — Yo. 
We now define 


A’ fe, y) =f, 9), AS y) = fe + hy) — I ») 


and in general 


(8) AS, y) = 2 ( - '(? jc + ih, y), 
and similarly , 
(9) ASG y) = > (= (re, y + jk). 


From (8) and (9) it follows that 
(10) A” A*f(x,y) = AS A’ f(x, y) 


(zh) (y.&) (y.k) (2.a) 


ix=0 j=O 


p @ p 
> >d(- yrre(?)(2) pce + ih, y + jk). 


If the partial derivatives of f(x, y) of order < p + q all exist at (xo, yo) and 


are continuous in S, it can be shown that 


a" ¢ 


(11) A” AS f (x0, yo) = h” k* = f (xo + pO: h, y + qOek) 
) 


a 4,.P-a 
(zo.h) (yo.k Oxo IYO 


where 0; = 6;(xo, yo; A, k) (0 < 0, <1; i =1, 2). 


Applying (10), by induction we can show for g = 0, 1, 2,..., p, 


Fi 
(12) ax® ay Bins (X,Y) 


= m(m, — 1)... (m1 —g+1)mo(me — 1)... (me -—-pP+q4+1) 


ni—@ Na—P+@ 1 y 
: > > A’ A?-? f — = 
) 


¥5=0 va=0 (2.m,~*) (y.m,-? 


2 


2 


Ypnua-a(t) P», na~ p+e(¥)- 
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3. The main theorem. We state Bernstein’s Theorem for functions of two 
variables. 


Lema 1. If f(x, y) is bounded in the square S, then at every point of continuity 
(x,y), 
lim Bh, n(x, ¥) = f(x, 9); 


the result holding uniformly in x and y if f(x, y) is continuous in S. 
The proof follows simply by extending the proof of the one-dimensional 


Bernstein’s Theorem to two dimensions. 


LemMA 2. [If all the partial derivatives of f(x, y) of order < p exist and are 
continuous in S, then 


a” , a” / 
(13) axa By, n(x, 9) > axay? aS (x,y) 


uniformly in S as m1, M2 approach infinity in any manner whatever. 
Proof. Applying the relations (12) and (11) we obtain 
a” . 
(14) dx" dy’ * Bi, ns (X, y) 


= (my —1)...(my — g + 1)mo(me —1)...(m2 —- p+ q4+1) 


1 ‘(2 y" Ri-@ e—P+¢ o” (2: 7 vy, tus) 
(1) Nn» 2d, 2, ax" dy” if ny +4 nm,’ Ns + (- 49) No 
Po..ac~q(*) Poon p+a(d)s 
where 
0 < 6..., <1, 0 < dus, < 1. 


Now, in the relation (14), the product of the factors outside the double sum 
will approach 1 as m, m2 approach infinity. We denote 0?f/dx*dy’~* by h(x, y) 
and observe that the difference between the double sum in (14) and 


h 
Ba.—¢.0.-p+e*: y) 


approaches zero uniformly in x and y as m,, mz — ©. But by the previous lemma, 
the last expression approaches 


uniformly in x and y. Hence the limit relation (13) follows, proving the lemma. 
LEMMA 3. 


Tal [u"(1 — u)"~"] = O(u)u”"(1 — u)””' 


where 
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Q(u) = > n‘(v — nu)’hi, ,(u), i,j >0,2%+j<l 
J 
and the hi ,(u) are polynomials in u independent of v and n. 


The proof of this lemma follows by induction, see [2]. 
We are now in a position to consider the theorem stated in the introduction. 


Proof of Theorem. Since the pth total differential of f(x, y) exists at (xo, yo), 
we have the representaion (6) with P(x, y) a polynomial of degree p and g(x, y) 
as defined by (7). We therefore obtain 


a” SU. a” P a ¢ 

axtayn * Ban, (Xo, Yo) — axiaye * Bans (Xo, yo) + axtays © Ban, (Xo, Yo) 
and by Lemma 2 we only need to show that the second term on the right-hand 
side approaches zero. Using the definition of g(x, y) in (7), it is sufficient to show 
that 


a” . 
(15) Sehayg* Bases Xo, Yo) — 0, k=0,1,2,...,p. 
By Lemma 3, we find that 
a” Be 
axhays * Bain, (Xo, yo) 
ns ns p—k 
= 2 Miso) Do WS.) DU ou( 24 n)(2 . x.) 
2iitsi<e 2ist ja<p—@ vi ==0 y,=0 nN, Ne nN; 
A 
(x - y») ni' (v1 — myxo)”*n3' (v2 — ne)"(™) a3" * xynre(*) 
No Vi V2 
° yo “(1 aa gapenrers 


and we can rewrite this in the form 


a : 
(16) axtayh* Brin. (X0, Yo) = a. hi..s(%0) D> i. 40) 


2i:+3:<¢ 2ist+ ja<p— 


ms Ns p—k k 
Vy ve V1 Ve 4, js 
> ag Peng — a a Se ee (vy = N Xo) 
¥,=0 v=0 nN, Neo ny Ne 


: 3" (ve — mayo) ”*P»,.n,(X0)Prs.ns (Yo) 


where the h*;,, ;, (xo) are independent of v; and m, and the h*,, ;, (yo) are indepen- 
dent of v2 and m2. Hence it is sufficient to show that the inner double sum in 
(16), which we denote by ¢, approaches zero under the restriction (3). 

Now for any ¢ > 0 there is a 6 such that 


(17) jory (1, y1)| << €, lx = Xo| < 6, ya —_ yo! < 6; k = 0, # ® coe »P 


(the choice of 6 is made more precise later on), and there also exists an M such 
that lay (x1, ¥1)| < M for0 <x, < 1,0 < y1 < 1. 
In absolute value the part of ¢ in (16) corresponding to 
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does not exceed 


Mm, Ms 
1—Ptk_ is—k |p—k+ fj, |R+Js 
ni No «>. ys . je MXo\" , \Vo — Neo} : Pr,.n, (Xo) Py, n, (Yo) 


¥,=0 y,=0 
and by (5), this last expression does not exceed 


1,—p+k+) (p—k+ jf.) tak +4(h+ Js) 
Coen} p+k+i(p ds n;" j ja 





< Cre( fa) tno o/te) **e-* 


As 24; + ji < q, 212 + j2 < p — Gg, we have 2%; + 2i2 + j1 + je — p < 0, and 
so this last expression does not exceed C;« for all sufficiently large m,, m2 which 
satisfy (3). 
In absolute value the part of ¢ in (16) corresponding to 
| 
| — x0| < 6, |=2 — yo > i 
| Ne 


| 1 | 


does not exceed 


ms 
at nM > |v a 83X0|” "Dyn, (Xo) Ze, lve _ nayo|’* "Dy, .n,(Yo)- 


{ v,=0 |» 


| ae —vo\>s 


Using (3) and the relations |» — nu| < n and 


) Do Prams (0) = 1, 


we conclude that the above expression does not exceed 


(18) ny nye Mint ns DS Pans (Yo)- 
= —yol>s 
Provided , m2 are sufficiently large and 0 < a < 3, the number 6 = m,;-* + 


n2~* satisfies the condition (17). Using the inequality (4) with this 4, we see that 
the expression (18) does not exceed 


t:—p+k is—k — + k+ tit tet JitJjetl 
n> P+ ns" M, n> ds Ns ds C./ni atjitds 


| which approaches zero as mz tends to infinity. 


) A similar conclusion holds for the part of corresponding to 
' | | | | 

| Vy | Ve 

“1 — x9|>8, | — yo| <8. 

| nN; = > | Me yo | 
Finally, the part of o corresponding to 

jm _ ,.| 2 _ | 
) lee xo| > 6, ~ yo| > 6 


is in absolute value not greater than 
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t,.—pt+k —k |\p—k : k . 
nin *M Don — miro!” “ty, n(x), D> lve = mayol** py, (Yo) 


\> 


|z* ze [as x* —yo|>8 
and by (4) this does not exceed 


n p—k+ J, k+Js 


—ptk__is—k 1 
n;' ?''n2" MC; - 
ny 
which approaches zero as m1, m2 — ©. Hence 


— a = BP . 

axtay, * mins (0 Yo) 

can be made arbitrarily small by taking m;, m2 large enough, satisfying (3). 
The theorem is now established. 

We shall now show that the limit relation (2) also holds at the four corners 
(0, 0), (0, 1), (1, 0) and (1, 1) of the square S. It will be sufficient to discuss the 
corner (0, 0). 

Writing 
d v, e—Oa te vs Ne—"s 
D,,»,(x, ¥) = gxt * (lL — x) aye ag 


we have D,,,,(0, 0) = 0 for 1» >¢g+1 or »>p—q+1. For » <q and 
ve < p — g by Leibniz’s theorem, we have 


D.,»,(0, 0) = («2 * a al — x) v. jdy” y 


(19) =~ Cnt "ns: *”, Ny, Nz @, 


Again we have to show that the limit relation (15) holds, but now at the point 
(0, 0) instead of (xo, yo). 

We have 

3” 


— . 
ax" dy’ * Buin, (X; y) 


r=y=0 


r n— p—k k 
" 3 3252 ... TE 1) ( %2 : 
on 2 E a( n)(n r) (2: y) () ( ") Daan y) 


and by (19), this in absolute value is 


1 V1 Ve hae Ay 
2 5 .¥s O—? —q-» 
~C:} >} Dal -=)h\— ~} ni'ny'nt "nt" 
y,=0 »,—0 nN, Ne ny Ne 


r=y=( 











—, 
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and since m, and mz, satisfy the inequality (3), the last expression does not 


exceed 
aia = 
k n,' Ns 


which approaches zero as m, m2—> © because a,(x, y) ~0 as x 0, y— 0. 
Hence we have the relation 


a a 
axa Bin. (x, ¥) | a axtayp tl ¥)| , 


. 


Cd 


vy =O r,=0 








’ 


=—y=0 


if the pth total differential of f(x, y) exists at (0, 0), provided m, and m; satisfy (3). 

We now consider an example to show that the relation (2) does not necessarily 
hold at a boundary point of S different from one of the four corners. Consider the 
point (xo, 0), 0 < x» < 1 and the function 


f(x, y) = e(y) |x — xol, O0O<x<1 


where «(y) is defined for 0 < y <1 and has the properties e(y) 0 and 
y te(y) ~ + © as y— 0. Take m, = n; = n. Obviously 


7 for v. = 1, 


d 
dy Po, .n(y) _ = 0 for Vo 


lI 
a 
ba 
x 


Therefore 





j S ve\d 
n a | raat 2) 2 Prva) | (2-0) 


n | | 
= >» 2 nXo| Poateade()m 


I 
iM: 


0 
By Baa, ¥) | — 


which approaches + © as m—> © because it is known [2] that forevery 0 <x <1 
there exists a constant C, > 0 such that 


D |v — nx| Prn(x) > CoV 2, n—- @, 


vy=0 
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A CONVOLUTION TRANSFORM ADMITTING 
AN INVERSION FORMULA OF 
INTEGRO-DIFFERENTIAL TYPE 


D. B. SUMNER 


1. Introduction. Following Wiener’s fundamental work on the Operational 
Calculus [9, pp. 557-584], Widder suggested [5, p. 219] that the inversion 
operator for the convolution transform 


(1) f(x) = Jee — t)¢(t) dt 


should be a suitable interpretation of E(D), where D = d/dx and the function 
E(s) is defined by 


(2) E(s) = frerew dt. 


Two methods of interpreting E(D) have been used. The first, which appeals 
only to real-variable methods, has been used by Widder and his collaborators 
[8, p. 659; 5, p. 217; 4] in cases where E(s) is entire with real zeros. The substance 
of the method is to express E(s) as 

lim E, (s), 


n+ 


where E,(s) is a product of m factors of the form 


(1-3) 


and where, under suitable convergence conditions, 


lim E,(D) - f(x) 


can be computed. 

The second method, which uses the complex variable, has been applied by 
Widder and Hirschmann [8, p. 692], and by the present author [2, p. 174] in 
cases where the entire function E(s) can be expressed as a Fourier integral. 
It is naturally of less general applicability. 

At the conclusion of their article Widder and Hirschmann conjectured the 
existence of convolution transforms, whose inversion operators are of integro- 
differential form, and not purely differential, as in all the then known cases. 
Widder has followed this by using such integro-differential operators to invert 
the Lambert transform [6, p. 171] and the Fourier sine transform [7, p. 119]. 


Received October 17, 1951. This investigation was carried out while the author held a 
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In the last-named work, he showed that the methods of [8] were still applicable, 
even though the function defined by the integral 


(2') Ae = foerew a 


“2 


was meromorphic. In both applications (to the Lambert and to the Fourier 
sine transform), a distinction was made between the factors F(D) and 
|[E(D)}"' of the inversion operator. 

The objects of this note are: 

(i) to show that transforms requiring an integro-differential invertor occur in 
quite simple cases, the example treated being one in which the kernel is a function 
little more complicated in form than the classical Stieltjes case [1, p. 473], 
where the kernel is (1 + e~*)'; 

(ii) to prove that the order in which the factors of the invertor are applied is 
material. 

We call F(D) the differentiating factor, and [E(D)}~' the integrating factor. 
We use the methods of the complex variable, and interpret the factors of the 
invertor by expressing them as integrals. 


2. Preliminary results. We consider the transform 
(3) f(x) = [He - Ho a, 


where H(x) = [e~* + 2e-* cos x8 + 1]-', and 0 < 8 < 1. By adapting classical 
methods due to Widder [3, pp. 10-11], after an exponential change of variable, 
the following theorem is easily seen to be true: 


THEOREM 1. Let f(x) be defined by (3), and o(t) be such that the integral 
converges for at least one value of x in the strip |¥x|\ < (8 + 1), then it converges 
for all such x, and converges uniformly in any compact subset of the strip, so that 
f(x) ts analytic in the strip. 

It is well known that 

* a ra sin rB(1 — s) 
(4) f er ED nv nde “qr pee (0 < Rs < 2). 


ri sin ms sin +8 


We take the reciprocal of this function, with s = D, as the invertor for equation 
(3), and we write 


F(D) = 222? = tim 5+. fe'*? — "| = lim F,(D), 
647 649 
sin rB sin 78 (° e "e? dv 


- —l — a<_om =F SS _ 
[EO = sisal —D)~ «8 J.i1+e" 


and shall as usual interpret e*?(f)x as f(x + k). 


3. The inversion theorem. Our main result is the following theorem: 


THEOREM 2. Let f(x) be defined by (3), with 0 < 8B <1, and o(t) be such 
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that the integral converges for at least one x in the strip |3x| < x(8 + 1). Then 


F(D) - (E(D)I"f(x) = Hot) + o@-)], 
whenever the right-hand side has a meaning. 


We have 


sin 8B ('“e 'f(x + v) dv 
xB - 1 + e 4 ’ 


° ) ae = 
op EB [ime+—omne 


= fixe ~ 1) dt, 


[E(D)I'f(=) 








where the interchange of the integrations is justified by the uniform convergence 
established in Theorem 1, and 


K(r) = =—s . ce aera 





On considering the integral 


1 e” dw 
eur &i1i- ev" 
taken round the rectangle with vertices R + iS, — R+ wif, and making 
R— o, we see that 





K(r) = [1+ 6°)". 


Thus 
all _ £° ot) dt 
(5) EOrye) = f° POS. 


This is the classical Stieltjes transform (after an exponential cliange of variable); 
and if we denote it by g(x), it is well known [4, p. 340] that 


lim F,(D) - g(x) = 3[o(x+) + o(x—)]. 


Theorem 2 is thus seen to be true. 

We do not here allow the value 8 = 0, since the integral (4) then has the 
value (1 — s)/sin xs, whose reciprocal has no poles. It is easily seen that in 
this case, after suitable changes of variable, the equation (3) takes the form 


_ Fv) dt 
h(x) = J. (x +0? 
which has been discussed by the present author [2]. 





4. Commutativity of the factors in the operator. That the order in which 
the factors F(D) and [E(D)}|“' were applied to the function defined by (3) is 
material, is brought out by the following example. 








ce 


iS 
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e', t>0, 
o() = “4 t<0. 


It is then easily seen that the corresponding 








f(x) = f H(x — t)-e‘ dt, 
& | sin +8 |, 
~ sin <B- ete e* + cos #81’ 
while 
F,(D)f(x) = = : | en arc tan( _ =B ) 
sin +8 e* "+ cos #8 





— “arc tan( ein =f ) —0 
eo + cos rB 


as 6—>+ x. The differentiating factor, if applied first, is thus seen to annihilate 
the function, a fact which makes pointless the subsequent application of the 
integrating factor. 
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SYSTEMS OF PARTIAL DIFFERENTIAL EQUATIONS 
WHICH GENERALIZE THE HEAT EQUATION 


I. I. HIRSCHMAN, Jr. 
1. Introduction. The distribution of heat in an infinite rod is closely bound 
up with the theory of the Weierstrass transform. This connection is exhibited 
most clearly in the theory of Widder [3]. Consider the heat equation 


1.1 Ss = —. 


Widder has shown that if u(x, h) satisfies 1.1 and if u(x,h) > 0 for — ©- <x < @, 
a <h < b then we have 


13 u(x, h”) = f k(x — y, hh” — h’)u(y, h’) dy 
ifa <h’,h” <b, h" > h’. Here 

1.3 k(x, t) = (4xt)* exp ( — x’ /4t) 

1.4 k(x, t) = (2x)~* j exp [ — ty’ — ixy] dy. 


Let &,, and §, be real m- and n-dimensional vector spaces respectively. We 
shall write the elements X, Y, Z of &,, and H, H’, H” of $, as column matrices, 


Xy1 y1 21 
X= (:" Y « (" Z = (* } 
Sas Vm Sa, 


hy hi \ ui 
H=|")) we=l|™) gral” I 
hy hn hy’ 
In place of u(x, h), we consider u(X, H) depending upon the m + n real variables 
X1,.. + 5Xmy My,...,h,, and instead of the heat equation 1.1, we consider the 


system of partial differential equations 


ou _ a°u ' 
: — = —— a; fc a | 
15 Oh, f5h1 ax, ax," . 


Here the a’s are real constants and 


Received November 12, 1951. Research supported by the Office of Naval Research. 
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i i 
Qizy = A5z. 


In the present paper we shall develop for such systems a theory analogous to that 
of Widder. 
Let Xj (m+1) be the space of real symmetric m X m matrices, 


tia tig «2. bam Gia ia «ee Ge 
T= (" too... tom T’ = (‘ fe... Oe ) 
tm tm see tmm tnt bin? see a 


We may define a partial ordering in Xjmim41) by writing 7’ > T if T’ — T is 
positive definite. We introduce corresponding to the system 1.5 a mapping 
H — H, of , into Timms, the (7, 7) entry of H, being defined as 
> ai, hi eo] Sar m). 
l=1 
We shall show that if uw(X, H) satisfies 1.5 and if u(X, H) > 0 for X € &,, 
H € N where MN is an open convex subset of H,, then 
1.6 u(X,H") = fi, R(X — ¥,H",— H')u(¥, H’) d¥ 
if H’, H” € R, H”, > H’,. The function k(X, T) is defined by the formula 
1.7 R(X, T) = (24) ™ ‘ exp [— Y*7Y — iY*X]dY; 


here T must be positive definite. An explicit formula for k(X, 7) (due to Cziiber) 
is given in §3. 
If A is any matrix, A* is its transpose; thus 


Y*X = _ Vi Xy, 


i=l 


y°7Y = dy bey Vj 


i, j= 


2. A special case. The equation of heat transfer in m-dimensional space is 


au OU du 
2.1 ah = u ax? or ah = Au. 


In this section we shall establish our theory for this special case. If m = 1 in 
2.1 we have the equation considered by Widder. It is to be noted that only 
relatively minor adjustments are needed to make the extension from 1 to m 
dimensions. We set 


2.2 R(X, t) = (4nt)"™ exp | = i x*x | 


k(X, t) = (2")™ 3 exp [ — tY*Y — iY*X] dV. 








120 I. I. HIRSCHMAN, JR. 


Let us write |X| for 
(x*X)! = ( xi)*. 
We begin with a discussion of the convolution transform 
2.3 V(X, h) = Pc: — Y,h)¢(Y)dY 


under radial symmetry. A function ¢ defined on &,, is said to be radically 
symmetric if ¢(X’) = ¢(X”) whenever |X’| = |X”|. If in formula 2.3 ¢(Y) is 
radially symmetric, then for each h, ¥(X, h) is also radially symmetric and the 
m-fold integration in our formula may be collapsed to a single integration with 
respect to a new kernel, which we need later and which we now compute. 
Since ¢(Y) depends only upon | Y| = p we write ¢(p) instead of ¢( Y). Similarly 
since ¥(X, h) depends only upon |X| =r and h, we write ¥(r, h) instead of 
V(X, h). 


Introducing spherical coordinates, 


¥1 =pcos di 
Ye = psin $1 Cos go 
Ys = /psin ¢2Sin ¢2 Cos ¢; 


Ym-1 = P sin o1 sin Pieces sin dm—2 COS Om—1 
Ym psin @) Sin do... SIN dm—2 SIN Pm—1, 


and setting x) = 7, X%. = x3;=.. = 0, we have 


¥(r, h) = yrds J. ¢: f (pen 7* 278 contatoed 144 


- SIN Dp 2d¢ +. - 1bm— 1 dp. 


From this we obtain 
$(m—3) 
v(r, t) i (4 ae an o(p)e~” +p? at, m— 17 fe rp COS, (24, ° in” *¢, d¢y. 
By Watson [2, p. 79] we have 
‘ Tp cosd, /2t m—2 ai Tm _ns b) Vr (2) 
f e sn dg; = (rp/4t) = Tym-1 ° 
We finally get 


2.4 v(r,t) = (2t)"* J eiroenniginy imei, (12) o(9) dp. 


The function u(X, h) is said to belong to S, u(X, h) € S, at (Xo, ho) if there 
exists an open set in E,, X H containing (Xo, ho), throughout which the partial 
derivatives 

Ou 0°u — 
ah’ Ox,dx, (i,j =1,...,m) 
are continuous and if 
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ot = Au for X = Xoth = he. 
We denote the set of points 
IX| < R, 0O<hec 
by De, and the set of points 
IX| < R, t= 0, 
by Bae. 
LEMMA 2a. If 
i. u(X,h) € S for (X,h) € Dre} 
2. lim u(X,h)>O0 for (Xo, ho) € Bre 


XX o Raho 
then u(X, h) > 0 in Dp. 
Suppose that there exists a point (X,,4,) € Dg, at which we have u(X,, hi) = 
— 1< 0. Form the function 
v(X, h) = u(X, h) + k(h — hy) 
where k > 0 is chosen so small that 
lim v(X, h) a> 4/ (Xo, ho) € Bre. 
X4Xo,h+he 


For 6 > 0 suitably chosen »(X, hk) > — 1 within a distance 6 of Bz,. Conse- 
quently the minimum of »(X, h) in Dg, is attained at some point (X2, hz) of 
Dg,. The function v(X, h) satisfies the partial differential equation 


ov 
2.5 Av = ah — k. 
At (Xo, 42) we must have 
2.6 Av > 0, 
Ov 
2.61 ah < 0. 


(If he # c the equality holds in 2.61.) Equations 2.5 and 2.6, 2.61 are in contra- 
diction. 


LemMA 2b. Jf 


1. u(X,h) € S, X€G@,,0<hk gc; 
2. lim wu(X,h) =0, Xo € En; 
XaXe, h 0+ 
3. M(r) = Lub. |u(X,h)|; 
|X|=r, 0<a<c 
4. M(r) = O(¢"), r++ 


for some a > 0; then u(X, H) = 0 for X € Gy, O<h Ce. 
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Take \ > 1 and form the auxiliary function 


|X|*+R* m | y\—}m+ |X| 
Ua(X, b) = AM(R) (xe) at eR RIE xm (RIX), 
Using the asymptotic formula 


2.7 I,~e/ V/ 2x (x + + @~), 
we find that [2] 


Un(X, h)|\xin2n — A(c/h)* M(R), R-@, 
uniformly for 0 < h < c. Thus when R is large 
Up(X, h) > M(R), IX|=R,O0<he. 


The function Upg (X, h) — u(X, h) belongs to S for X € &,, 0 < kh < c. When 
R is large 





lim Upg(X,h) — u(X,h) > 0, (Xo, ho) € Bre. 


XaXe, Raho | 

By Lemma 2a we have | 

u(X,h) < Ur(X, h), (X,h) € Dre. 

Fix X and hk and let R-— + @. Making use of 2.7 we find that u(X, h) = 0 

for 0 <h < 1/4a, X € &,. If 4a < 1/c, our proof is complete. Otherwise 

repeat the above argument with u(X, h) replaced by u(X, h + 1/4a), and so 
forth. 

L(€,) is the class of functions ¢(X) defined for X € ©,, and such that 


Jog. \#eo ax 
exists. 


Exactly as in [3] we may establish 


Se 


LEMMA 2c. If 
1. o(X)e**"" € L(G,) for some a > 0, 


2. F(X, h) = ih. k(X — Y,h)o(Y) dy, 
then F(X, h) is defined and belongs to S in the strip 0 < h < 1/4a. 
LemMMA 2d. Under the assumptions of Lemma 2c, we have 
lim ¢(X)< lim F(X,h); lim F(X,h) < lim ¢(X). 
X+X oh 0+ XX 


XX XX oh 0+ 


LemMA 2e. If o(Y) is integrable for |Y| <A then for any c> 0, we have 


bad 
LemMA 2f. If 
1. u(X,h) € S, O0<hkh<c,X € &; 
2. u(X, hk) > 0, 0<hk<c,X € &; 


| 
lim (ut, R(X — Y,h)o(Y) dY 


0<a<c | ly|<a 





|X |+00 


then we have 
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< R(X — Y, 8’ )u(Y, 6) dY < u(X, & + 8), 
0<6,0<8,5+8 <c¢,X € &,. 
Consider the function 
v4(X,h) = u(X,h + 6) — f , 
ly 
where A is a positive constant. If 8’ < c’ < c — 6 then 2,4(X, h) belongs to S 
for (X, h) € Dee for any R. Let / > 0 be given. Using Lemmas 2d and 2e we 
see that, if R is sufficiently large, 
lim wv,(X,h) > —/, (Xo, ho) € Bre. 


XiX oe haho 


k(X — Y, h)u(Y, 8) dY 
<A 


Lemma 2a implies that »,(X, h) > —/ for (X, h) € Dg. Letting 1-0 + 
and R-— + o, we find that »,(X, h) > 0 for X € &,, 0 <A < c’. Setting 
h = & we have - 


J ) R(X — Y, &’)u(Y, 6) dY < u(X, & + 8). 
Y|<A 
Letting A increase without limit, we obtain our desired result. 


LEMMA 2g. If 


1. u(X,h) € S, O<h<c,X E€ Ca; 
2. u(X, h) > 0, O<ch<c,X € &; 
3. u(X, 0) = 0, X € &; 
then u(X,h) = 0 Oc h<c,X € Ee. 


h 
w(X,h) = f u(X, t) dt. 
0 


It is easily verified that w(X,h) € Sfor0 Ch <c,X € Ep. We have w(X, hy) 
> w(X, he) if c > hy > Ae > O. It follows that 

2 w(X, h) > 0 

oh , 


and hence that Aw(X, h) > 0 for 0 < kh < c; equivalently w(X, h) is a subhar- 
monic function of X for each value of h,0 < h < c. 
Let 6 be an arbitrary number 0 < 6 < cand let 6’ be such that 0 < 8 <c¢ — 6. 
By Lemma 2f we have 
w(0, 8’ + 5) > (445) fe eV ly, a) dY. 


Again we have 


i) e848 ory, 8’) dY > M(X) f w(Y, 8’) dY 
G,, ly—x|<1 


where 
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M(X) =g.b. fe !”'*"*). 
ly—x|<1 
Since w(Y, 6’) is subharmonic 


f w(Y,8’)dY > wuiX, 3’). 
|y—x|<1 
It follows that 
w(X, 8’) = O(exp a|X|’) 
for a > 1/48. Since 


] 
ah w (X,h) >0 
we find that 
lu.b. |w(X, 8)| = O(e*’). 


0<a<é’,|x|=R 


Applying Lemma 2b we see that w(X, hk) = 0 for 0 ch < 8’, X € &,,. Since 
i’ < c is arbitrary we have w(X, h) = 0 for0 <A <c, X € &,. This in turn 
implies that u(X, h) = 0 for X € ©, Och <e. 


THEOREM 2h. If 
1. u(X,h) € S forX €Gn,a<h <b; 
2. u(X,h) > 0 forX €@,,a<h<b; 
2648, 8 <0, B > 8; 
then 
u(X, kh’) = i. R(X — Y,h"” — h’)u(Y, h’) dy. 


By Lemma 2f we have 
Seg tx — Y,h— h’)u(Y, h’) < u(X, h) 
for h’ < h < 6. This together with Lemma 2c implies that 
Seg — VYjh—W’)u(Y,h’)dYes 
for h’ <h <b, X € &,. Thus 
v(X, h) = u(X, h) - Jog ex — Y,h—h’)u(Y, h’) dV 


belongs to S and is non-negative for h’ <h <b, X € G,». Moreover, by 
Lemma 2d, 
lm v(X,h) = 0. 
hah’, YX. 
Lemma 2g implies that v(X, h) = 0. Setting h = h”, we obtain the desired 
result. 


3. The main theorem. Let JT € Zjnim41) be positive definite and let X € &,. 
We set 
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3.1 k(X, T) = (2%) ™ exp [— Y*TY — iY*X] dY. 


Cn 


Since T is positive definite, there exists a constant « > 0 such that Y*7Y > 
e|¥|*. This insures the validity of our definition. Let ® = [4);,s01,.... n bea 
real unimodular matrix. We assert that 


3.2 k(®"X, T) = k(X, 6T@"*). 
We have 


k(@"'X, T) = (24)™ Se exp [— Y*TY — i¥*@"X] dV. 
Making the change of variable Y = ©*Z we obtain 
k(@"X, T) = (2")™ GP [— Z*@T*Z — iZX]dZ = k(X, T4*). 


This formula may be used to compute k(X, T) explicitly, see [1, p. 185]. By 
this method, Cziiber has shown that 


3.3 b(X, T) _ exe( |2%"| /airi) 
(4x)™|7|' 





In particular if T is a diagonal matrix with equal entries, 


then k(X, T) is equal to k(X, r). 
Consider the system of partial differential equations 
ou = d*u ' I - 
dh, ijn Ox, ax, ** (ai, = ay; 1=1,...,m), 


3.4 





where the a’s are real constants. The function u(X, H) is said to belong to S(a) 
at (Xo, Ho) if there exists an open set in E,, X H, containing (Xo, Ho) throughout 
which the partial derivatives 


2 





ou u =F 
ah, = les +++)» Ox, ax, ( I = 11--- +m) 


are continuous and if equation 3.4 holds for X = Xo, H = Hp. 


LemMA 3a. [/f: 
1. u(X, H) € S(a), where N is an open subset of Hy, X € En, HEN; 


2. © = [Pisle. sunt, .... m is a real unimodular matrix; 
3. w(X, h) = u(@"X, h); 
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3.5 a. FS 2-k, 


for X © En, H € RN, where 
bis = > aij dai op j- 


i, j=l 
(Briefly we have w(X, H) € S(b) for X € ©, H € NR.) 
We have w(X, H) = u(Y, H) where X = @Y. Now 
du _ dw du Sw 
oh, 3 oh,’ OY; OY; -_ = OXa OX 
Substituting in 3.4 we obtain 3.5. 
If H € §,, then H — H, is a mapping of , into Tymim+1), the (7, 7) entry of 
H, being 





ai pj. 


D hi bi. 5. 
l=1 
It is easily verified that 
3.6 oH, &* = H,. 


LEMMA 3b. [/f: 
1. u(X, H) is continuous for X € En, H € N where N is an 
open subset of D,, 


2. H’,H” €R, HH". > H’., 
3. u(X, H”) = 1. k(Y — X, H”, — H’.)u(Y, H’) dy, 


4. w(X, H) = u(@'X, Hi), 


then w(X, H) is defined and continuous for X € En, H € N, H", — H’,> 0, 
and 


w(X,H") = fe k(Y — X, H"”, — H’,)w(Y, H,) dy. 
In 3, replace X by ®-X and make the change of variable Y = ®~'Z to obtain 


u(@'X, H"”) = 4 k(@"Z — ®°X, H", — H',)u(®'Z, H’) dZ. 


w(X, H") = & k(@"Z — @'X, HH", — H’,)w(Z, H’) dZ. 


By 3.2 and 3.6 we have 
k(@'Z — @°X,H", — H',) = k(Z — X, 6H", &* — oH’, *), 


k(Z 5 ow X, H", — H’,). 


Our lemma follows. 





—S =. —_—---—_—_——— ee 
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THEOREM 3c. [Jf: 
1. u(X, H) € S(a), X € En, H € R, where N is an open 
convex subset of Dz, 
2. u(X, H) > 0, X€G@, HER 
3. H’,H” ER, a".> Fe 
then 


u(X, HH”) = fe k(Y — X, H”, — H',)u(Y, H’) dy. 


Let ® be a real unimodular m X m matrix. Because of Lemmas 3a and 3b, 
it is sufficient to establish the corresponding relation for w(X,H) = u(@-\X, H), 
w(X,H") = ‘ k(Y — X, H", — H’,)w(Y, A’) dy. 

By 3.6 we have 
H", —_ H', = O(E “n = H,)%*. 
It is possible to choose @ so that H”, — H’, is a diagonal matrix with equal 
entries, i.e., 
-_ ee 
aan Oy ’ 
00 T 


Having chosen # in this manner, we consider 


w |x. ( a ‘) H’ + "| = » (X,?). 
T T 


Since ® is convex and open, v(X, #) is defined for a < t < b where a < 0, 


b> +r. We have 
Ov — dw|1.,, 
mo ella, w') | 





t= T 
Since 
ow — a*w i “ av i 
— = — = —----—§ 
Oh, Fa; Ax, ax, b4s », Ox, Ax, * 
we find that 
w_1 S[y, ( » ) | 
at b> sie iit x, Ox, 
ov 
_ Av. 


Thus 0(X,?#) € SforX € ©,,a < t < bwherea < 0,5 > +r. Applying Theorem 
2h we have 


»(X, r) = fe k(X — Y, r)v(¥, 0) dY. 


Now 
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v(X, r) = w(X, H”), o(X, O) = w(X, A’). 


w(X, HH”) = Se R(X — Y, r)w(Y, H’) dY, 


w(X, H") = 4 k(X — Y,H”", — H’,)w(Y, A’) dy, 


a relation which we have seen to be quivalent to our theorem. 

There are some systems of equations 3.4 for which Theorem 3c gives no infor- 
mation. This is because there do not exist real vectors H’ and H” such that 
H”, > H’,. The system consisting of the single equation 


is of this type. 
Making use of Theorem 3c and the concept of weak compactness, we may 
demonstrate the following result. 


THEOREM 3d. Let N be an open subset of , such that H', H” © N imply 
NH” + WH’ €N for 0< rd’, O< dA", XY +A” <1, and such that HEN 
implies that H, is positive definite. If: 

1. u(X, H) € S (@), X€G@, HER 


2. u(X, H) > 0, X€G@G, HER 
then 


u(X, H) = Jig kX — Y, H,) dm(Y) X €@,, HER, 


where m(Y) is a non-negative measure defined in the a-field of Borel sets of En. 


REFERENCES 
1. S. Bochner, Fouriersche Integrale (Leipzig, 1932). 
2. G. N. Watson, A Treatise on the theory of Bessel functions (Cambridge, 1922). 
3. D. 


V. Widder, Positive temperatures in an infinite rod, Trans. Amer. Math. Soc., 55 (1944), 
85-95. 


Washington University, 
Saint Louis 








